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1.  INTRODUCTION:  We  had  proposed  to  evaluate  lipid  droplets  (LDs)  as  a  disease  marker  and 
therapeutically  target  them  in  ovarian  cancer.  Our  preliminary  results  showed  that  stable 
downregulation  of  HSulf-1  in  OV202  (Shl/2)  and  OV2223  lead  to  increased  LDs  as  evidenced  by 
BODIPY  staining  and  transmission  electron  microscopy  (TEM)  analysis  compared  to  non-targeted 
control  (NTC)  cells.  Conversely,  enhanced  expression  of  HSulf-1  in  TOV21G  and  SKOV3  cells  also 
resulted  in  reducing  the  number  of  LDs  in  these  cells  compared  to  vector  transfected  controls. 
Additionally,  our  microarray  analysis  identified  cPLA2a,  a  protein  involved  in  LD  biogenesis  as  one  of 
the  genes  that  was  significantly  upregulated  OV202  Shi  and  Sh2  cells  compared  to  NTC  cells  .Based 
on  this  finding  and  work  done  by  others  we  hypothesized  that  the  activation  of  extracellular  signal- 
regulated  kinase  1/2  (ERK1/2)  by  increased  growth  factor  signaling  leads  to  the  activation  of  cytosolic 
phospholipase  A2  alpha  (cPLA2a),  which  culminated  in  a  higher  accumulation  of  LDs  in  ovarian 
cancer  cells.  We  further  hypothesized  that  tumors  with  activated  cPLA2a  can  most  likely  be  treated 
with  targeted  therapy  against  LDs  in  addition  to  standard  chemotherapy.  We  also  proposed  to 
elucidate  the  role  of  cPLA2a  in  LD  formation.  The  role  of  LDs  in  chemoresistance  and  metastasis  will 
be  examined  using  in  vitro  and  in  vivo  models.  Lastly,  we  proposed  to  determine  the  expression  of 
phospho-cPLA2a  (p-cPLA2a)  in  primary  ovarian  tumors  and  correlate  it  with  clinical  outcome. 

2.  KEYWORDS:  Ovarian  Cancer,  TEM-Transmission  Electron  Micrographs,  Lipid  droplets,  C-PLA2 
(PLA2G4A)-cytosolic  phospholipase  A2,  iPLA2-  Ca++-independent  phospholipase  A2,  Arachidonic 
Acid,  AACCOF3-  arachidonyl  trifluoromethyl  ketone,  MAFP-  methyl  arachidonyl  fluorophosphonate, 
BODIPY-  Boron-dipyrromethene,  BEL,  Biomol-  Bromoenol  lactone,  DEDA-  7,7-  dimethyl-5, 8- 
eicosadienoic  acid,  FFA-Free  Fatty  Acids,  TMA-  Tissue  microarray,  TTR-Time  to  Recurrence,  PFS- 
Progression  Free  Survival,  OS-Overall  Survival 

3.  OVERALL  PROJECT  SUMMARY: 

Milestone  ^  Animal  Use  Approval-Prepare  IACUC  protocol  for  animal  related  studies. 

The  IACUC  protocol  #A69212-Therapeutic  targeting  of  lipid  droplets  as  a  disease  marker  in  ovarian  cancer  is 
approved  by  the  Animal  Care  Committee  at  Mayo  Clinic. 

Aiml :  Elucidate  the  role  of  cPLA2a  in  LD  formation. 

We  previously  reported  that  stable  genetic  downregulation  of  HSulf-1  in  OV202  cells  achieved  by  two  different 
shRNAs  targeting  HSulf-1  (OV202Sh1  and  Sh2  cells)  displayed  increased  LDs  compared  to  non-targeted 
control  transduced  cells  (NTC)(7).  We  further  confirmed  these  findings  using  OV202Sh1  to  OV202Sh1  Clone 
Cl 7  cells  (generated  by  re-expression  of  CMV  driven  HSulf-1  construct  in  Shi  cells.)  HSulf-1  levels  were 
verified  in  these  clones  by  western  blot  analysis  as  shown  in  Fig  1A.  Subsequently,  lipid  droplets  were 
monitored  by  Bodipy  staining  in  these  cells.  Rescue  of  HSulf-1  in  OV202Sh1CI7  cells  reduced  the  number  of 
LDs  compared  to  OV202Sh1  and  Sh2  cells  (Fig. IB,  Insets).  Mean  fluorescence  intensity  of  Bodipy  staining  in 
these  cells  is  shown  in  Fig. 1C.  We  further  performed  quantification  of  the  changes  in  the  LDs  by  Transmission 
Electron  Microscopy  (TEM).  Data  showed  significantly  higher  numbers  of  LDs  in  OV202Sh1  and  Sh2  cells 
compared  to  NTC  cells  (Figs.  ID  and  E).  Please  note  that  in  our  previous  (Year  2  short  report)  we  had 
indicated  that  increased  lipid  droplets  in  cells  with  loss  of  HSulf-1  are  associated  with  inhibition  of  autophagy. 
We  pursued  this  angle  since  loss  of  LDs  by  lipophagy  is  a  very  novel  concept  not  tested  before  in  ovarian 
cancer  and  therefore  we  had  indicated  that  we  will  focus  on  the  mechanistic  aspects  of  the  connection 
between  lipid  droplets  and  autophagy  in  ovarian  cancer. 

The  increased  lipid  droplets  indicate  the  nutrient  rich  or  energy  proficient  state  of  the  cells  and  therefore  we 
hypothesized  that  cells  will  not  activate  catabolic  process  such  as  autophagy.  Therefore,  to  verify,  we 
determined  the  extent  of  autophagy  via  TEM  analysis  in  these  cells  and  observed  that  HSulf-1  depleted 
OV202Sh1  and  Sh2  cells  exhibit  less  autophagic  vesicles  (AV)  when  compared  to  OV202NTC  cells.  Also, 
rescue  of  HSulf-1  in  Sh1CI7  cells  resulted  in  increased  AV  compared  OV202Sh1  cells.  To  rule  out  the 
possibility  of  cell  specific  effects,  downregulation  of  HSulf-1  in  TOV2223  cells  with  ShRNA  (Fig. IF)  also 
showed  increased  numbers  of  LDs  in  TOV2223Sh1  cells  compared  to  TOV2223NTC  cells  by  both  Bodipy 
staining  (Fig.lG)  and  TEM  analysis  (Fig.  1H)  indicating  that  these  metabolic  alterations  were  not  unique  to  the 
OV202  cell  line.  Similarly,  a  high  number  of  AVs  was  observed  in  both  OV202NTC  cells  and  TOV2223NTC 
cells  as  compared  to  corresponding  HSulf-1  depleted  cells.  The  quantification  of  LDs  and  AVs  in  TOV2223 
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Figurel:  HSulf-1  loss  promotes  lipid  droplet  biogenesis  and  defective  autophagy  in  OV202,  TOV2223  and  in 
HSulf-1  knockout  mouse  embryonic  fibroblast  cells:  Protein  expression  of  HSulf-1  was  assessed  by  western  blot  in 
OV202  clones  (NTC,  Shi,  Sh2  and  CI7  cells).  Actin  was  probed  as  a  loading  control.  B.  OV202-NTC,  Shi,  Sh2  and  CI7 
cells  were  labeled  with  Bodipy  493/503  (green)  and  DAPI  (Blue)  to  determine  the  cytoplasmic  lipid  droplets  (LD)  and 
nuclei  respectively.  C.  Mean  fluorescence  of  Bodipy  in  NTC,  Shi ,  Sh2  and  CI7  is  measured  using  Carl  Zeiss  Zen 
(2009)  software  and  the  data  is  plotted  as  a  bar  diagram.  D.  Representative  trans-electron  microscopic  (TEM)  images 
of  NTC,  Shi,  Sh2  and  CI7  cells  are  shown;  cytoplasmic  LDs  are  marked  with  green  arrows,  whereas  red  arrows 
indicate  autophagic  vesicles  (AV).  E.  LDs  and  AVs  quantified  in  NTC,  Shi  and  Sh2  cells  (total  25  cells  counted  in  each 
category)  demonstrated  in  bar  graph.  F.  Protein  expression  of  HSulf-1  in  TOV2223  NTC  and  Shi  cells  is  detected  by 
immunoblot  analysis;  b-actin  is  probed  as  a  loading  control.  G.  Bodipy  staining  performed  in  TOV2223  NTC  and  Shi 
cells  to  detect  LDs.  H.  Representative  TEM  images  of  TOV2223  NTC  and  Shi  cells  were  shown.  I.  AVs  and  LDs  in 
TOV2223  cells  are  quantified  J.  Protein  expression  of  HSulf-1  in  Mouse  Embryonic  Fibroblast  (MEF)  wild  type  (WT) 
and  HSulf-1  knock  out  (Sulf-1  _/")  cells  is  detected  by  western  blotting.  b-Tubulin  has  been  used  as  a  loading  control.  K. 
MEF  cells  were  labeled  with  Bodipy  493/503  (green)  to  determine  the  cytoplasmic  lipid  droplets  (LD). 


cells  is  shown  in  Fig.  II.  Furthermore,  HSulf-1  knockout  (KO)  mouse  embryonic  fibroblasts  (MEFs)  also 
displayed  increased  numbers  of  LDs  compared  to  wild-type  cells. 


Similarly  upon  ectopic  expression  of  HSulf-1  in  SKOV3  cells,  there  was  a  significant  decrease  in  the  numbers 
of  LDs  compared  to  vector  transfected  controls  (Fig. SI).  In  contrast,  vector  transfected  SKOV3  cells  exhibited 
higher  degree  of  AVs  as  quantified  through  TEM  analysis,  compared  to  HSulf-1  transfected  SKOV3  cells  (Figs 
SIB  and  D  respectively).  Quantification  of  LDs  and  AVs  in  25  cells  is  shown  in  Fig. SIC. 


Fig  SI:  A.  Protein  expression  of  HSulf-1  was  determined  by 
western  blot  in  SKV03  vector  and  cl7  cells.  GAPDH  was  probed 
as  loading  control  B.  Representative  electron  micrograph 
images  of  SKOV3  vector  and  cl7  cells  are  shown;  LDs  and  AVs 
are  indicated  by  green  and  red  arrows  respectively.  C. 
Quantification  of  LDs  and  AVs  from  electron  micrograph  images 
are  determined  in  SKOV3  vector  and  cl7  cells  and  are 
demonstrated  in  a  bar  diagram.  D.  SKOV3  vector  and  cl7  cells 
were  labeled  with  Bodipy  493/503  (green)  and  DAPI  (blue)  to 
determine  the  cytoplasmic  LDs  and  nuclei  respectively 
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Milestone  2: 

Task  1:  Determine  if  increased  ERK  activity  in  OVCA  cells  results  in  increased  cPLA2a  activity 

Our  data  indicate  that  downregulation  of  HSulf-1  in  OV202  cells  results  in  increased  ERK  activity  (He  et  al  Int 
J  Cancer).  Conversely,  re-expression  of  HSulf-1  in  HSulf-1  deficient  cells  results  in  attenuation  of  growth 
factor  mediated  Erk  activity  (2-4).  To  determine  if  increased  ERK  activity  is  required  for  the  activation  of 
cPLA2a  as  previously  reported  (5,  6)  we  had  proposed  to  treat  OV202Sh1  and  HSulf-1  deficient  parental 
SKOV3  cells  (with  higher  levels  of  p-ERK)  with  MAPK  inhibitor  PD98059  in  a  dose  and  time  dependent 
manner  and  determine  if  this  will  result  in  the  inhibition  of  cPLA2a  activity  by  western  blot  analysis.  We  had 
proposed  to  determine  under  the  optimal  condition  of  ERK  inhibition  by  BODIPY  staining  and  TEM  analysis. 
Additionally,  to  determine  if  growth  factor  mediated  upregulation  of  ERK  activity  is  critical  for  increased 
cPLA2a  activity,  we  had  proposed  to  treat  OV202  NTC  and  Shi  cells  and  HSulf-1  deficient  parental  SKOV3 
vec  and  CI7  cells  with  HSulf-1  expression  with  HB-EGF  (50ng/ml)  for  12  to  24  h  and  examine  the  LD 
accumulation.  Shorter  treatment  (10-60min)  with  HB-EGF  will  be  used  to  examine  cPLA2a  activity  using  p- 
cPLA2a  (Ser505)  antibody  by  immunoblot.  Thus  the  rationale  for  this  experiment  is  to  inhibit  the  activity  of 
ERK  which  has  been  shown  to  activate  cPLA2a  at  serine  505  critical  for  LD  biogenesis.  Our  western  blot 
analysis  showed  that  SKOV3  cells  express  very  low  levels  of  total  cPLA2  and  no  p-cPLA2. 

_cPLA2  is  activated  by  phosphorylation 
at  Ser  residues  505,  515,  and  727(7). 
Within  these  sites,  Ser-505 
phosphorylation  leads  to  increased 
catalytic  activity,  (8)  which  depends  on 
MAPK,  ERK  and  p38  .(8-10)  We 
previously  showed  that  downregulation 
of  HSulf-1  in  OV202  cells  results  in 
increased  ERK  activity.(  7  7)  Conversely, 
re-expression  of  HSulf-1  in  HSulf-1 
deficient  cells  results  in  attenuation  of 
growth  factor  mediated  ERK  activity. (3, 
12,  13)  To  determine  if  inhibition  of  ERK 
activity  will  alter  p-cPLA2  levels, 
OV202NTC,  Shi  and  Sh2  cells  were 
treated  with  lOpM  U0126  (MEK  inhibitor)  for  24  hrs.  As  shown  in  figure  2A  (Panels  2  and  4),  treatment  with 
U0126  resulted  in  the  attenuation  of  p-ERK  and  p-cPLA2ser505  levels  in  OV202Sh1  and  Sh2  cells  compared  to 
NTC  cells.  Bodipy  staining  of  OV202  Shl-and  -Sh2  cells  show  figure  3B  shows  reduced  LD  accumulation  as 
a  result  of  U0126  mediated  inactivation  of  ERK  and  p-cPLA2  (F  ig.2B). 

Task  2:  Identification  of  role  of  cPLA2  in  regulating  LD  biogenesis  in  OVCA  cells. 

Milestone  3:  Pharmacological  approach:  We  will  examine  if  cPLA2a  plays  an  important  role  in  LD  biogenesis 
in  OV202NTC  and  Sh  clones,  and  in  SKOV3  vector  and  CI7.  Using  the  pharmacological  inhibitor  approach,  we 
had  proposed  to  use  cPLA2  specific  inhibitors  including  methyl  arachidonyl  fluorophosphonate  (MAFP), 
arachidonyl  trifluoromethyl  ketone  (AACOCF3)  from  Calbiochem.  To  examine  if  these  inhibitors  are  inhibiting 
cPLA2a  activity,  OV202  NTC  Sh  clones  and  SKOV3  -vec  and  CI7  cells  will  be  labeled  with  0.1  pci  of 
radioactive  arachidonic  acid  (AA)  ([3H]  AA,  ARC)  overnight.  Next  day,  cells  will  be  washed  with  PBS  three 
times  and  then  cells  will  be  treated  with  10pM  of  MAFP  and  AACOCF3.  Release  of  radioactive  AA  will  be 
examined  in  cell  supernatant  at  various  time  periods  from  2-24h.  After  standardization,  we  will  examine  the 
number  of  LD  in  treated  and  untreated  Sh  clones  with  cPLA2  inhibitors  after  24h  of  treatment  to  determine  if 
inhibition  of  cPLA2  leads  to  abrogation  of  LD  biogenesis.  As  a  control  for  these  experiments,  we  will  determine 
AA  release  under  similar  experimental  conditions  in  the  HSulf-1  rescued  Sh1C7  cells. 

Our  western  blot  analysis  showed  that  SKOV3  cells  express  very  low  levels  of  total  cPLA2  and  no  p- 
cPLA2.  For  this  reason,  we  scanned  other  ovarian  cancer  cell  lines  with  loss  of  HSulf-1  expression  and 
determined  that  OV2008  express  high  levels  of  p-cPLA2  and  decided  to  use  this  cell  line  in  addition  to 
OV202Sh  cells. 


Untreated  20uM  U0126 


U0126  (IOuM) 
p-ERK 


NTC 


Shi 


Sh2 


Figure  2 


OV202Sh1 


OV202Sh2 


Figure  2:  Inhibition  of  growth  factor  mediated  signaling  reduced 
p-cPLA2  and  LD  biogenesis:  A.  OV202  NTC,  Shi  and  Sh2  cells 
were  treated  with  20mm  of  U0126,  an  MEK  inhibitor  for  24  hrs  followed 
by  western  blot  analysis  to  detect  the  protein  expressions  of  phospho- 
ERK,  total-ERK,  phospho-cPLA2  and  total-cPLA2.  B.  Bodipy  staining 
for  LDs  in  OV202Sh1  and  Sh2  cells  after  treatment  with  20  mM  U0126. 
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cPLA2  is  activated/phosphorylated  in  HSulf-1  depleted  ovarian  cancer  cells  Since  our  initial  observation 
indicated  increased  lipid  droplet  accumulation  upon  HSulf-1  depletion  in  ovarian  cancer  cells,  we  next  wanted 
to  determine  the  underlying  mechanisms.  LD  often  accumulates  from  anabolic  process  known  as  lipid 
biosynthesis.  Previously  we  showed  that  lipid  droplet  associated  proteins  were  upregulated  in  HSulf-1  depleted 
cells  ( 1 ).  Among  them,  one  of  the  key  proteins  involved  in  lipid  droplet  formation  is  cPLA2a  previously  shown 
to  be  essential  for  LD  biogenesis. ( 14)  cPLA2a  is  activated  and  phosphorylation  at  Ser505  by  p-ERK  has  been 
shown  to  be  critical  for  cPLA2  mediated  LD  biogenesis. (15)  Therefore,  we  first  determined  the  activated/ 
phosphorylated  levels  of  cPLA2a  by  Immunoblot  analysis.  Our  data  show  that  p-cPLA2ser505  was  clearly 
increased  in  OV202Sh1  and  to  a  lesser  extent  in  OV202Sh2  cells  (Fig.  3A).  Re-expression  of  HSulf-1  also 
reduced  the  cPLA2  phosphorylation,  indicating  that  phosphorylation  was  reversed  by  increasing  levels  of 
HSulf-1. 

To  further  investigate  the  role  of  cPLA2  activation  and  lipid  droplet  accumulation/biosynthesis,  we 
treated  OV202Sh1  and  OV2008  cells  with  a  cPLA2  specific  inhibitor,  AACOCF3  (10  and  20pM).  Western  blot 
analysis  show  decreased  levels  of  p-  cPLA2  in  inhibitor  treated  cells  (Figs.3B  and  C).  Bodipy  staining  upon 
inhibition  of  cPLA2  activity  with  10pM  of  AACOCF3  and  MAFP  (another  cPLA2  inhibitor)  in  OV202Sh1  and  -2 

cells  (Fig.  3D)  showed  almost 
complete  inhibition  of  LD 
biogenesis  in  these  cells  compared 
to  untreated  control  cells. 
Consistent  with  this  data,  TEM 
analysis  of  OV202Sh1  cells  with 
AACCOF3  and  MAFP  also 
revealed  significantly  lower 
number  of  LDs  (Fig.3E,  with 
quantitation  of  LDs  in  25  cells 
shown  next  to  the  representative 
TEM  images).  Furthermore, 
transient  downregulation  of  cPLA2 
expression  with  two  different 
ShRNA  (ShcP-1  and  ShcP-2) 
against  cPLA2  (Fig.  3F)  in 
OV202Sh1  cells  resulted  in  the 
decrease  of  LDs  as  revealed  by 
Bodipy  staining  (Fig.  3G). 
Similarly,  stable  downregulation  of 
cPLA2  in  OV2008  cells  (Fig.3H) 
also  resulted  in  reduced  numbers 
of  LDs  (Fig. 31).  To  determine 
whether  cPLA2  inhibitors  were 
able  to  inhibit  cPLA2  activity,  we 


Figure  3:  Inhibition  of  cPLA2  attenuates  LD  biogenesis  in  OV202Sh1  cells.  A.  Immunoblot  analysis  shows 
phospho-cPLA2  and  total-cPLA2  levels  in  OV202  NTC,  Shi,  Sh2  and  CI7  cells.  B.  &  C.  Immunoblot  analysis  shows 
phospho-cPLA2  and  total-cPLA2  levels  in  OV202  Shi  and  OV2008  cells  after  treatment  with  10mM  and  20mM  of 
cPLA2  inhibitor,  AACOCF3.  D.  Shi  and  SH2  cells  were  treated  with  either  10mM  AACOCF3  or  MAFP  for  24  hrs 
followed  by  Bodipy  (493/503)  staining  to  show  the  LDs.  E.  Representative  TEM  images  of  OV202  Shi  cells  untreated 
and  treated  with  lOmM  AACOCF3  and  MAFP  for  24  hrs.  LDs  are  marked  with  green  arrows,  whereas  red  arrows 
indicate  autophagic  vesicles  (AV).  LDs  are  quantified  from  25  cells  from  untreated  and  treated  cells  from  TEM  images 
and  demonstrated  in  a  bar  diagram.  F.  &  H.  Total  cPLA2  expression  was  transiently  downregulated  with  two  different 
ShRNAs  targeting  the  open  reading  frame  of  cPLA2  in  OV202Sh1  and  OV2008  cells.  Immunoblot  analysis  shows  total 
cPLA2  expression  in  non-targeted  control  transduced  (NTC),  ShcP-1  and  ShcP-2  cells.  G.  &  I.  Bodipy  staining  of  LDs 
in  OV202  Shi  and  OV2008  cells  after  transiently  transfecting  them  with  either  empty  vector  or  sh-cPLA2.  J. 
Arachidonic  acid  (AA)  release  is  evaluated  in  OV202  NTC,  Shi,  Sh2  and  CI7  cells  treated  or  untreated  with  10mM 
AACOCF3.  Cells  were  incubated  with  H3-AA  under  serum  starved  condition  for  24  hrs.  Fresh  medium  was  added  to  the 
cells  after  washing  and  aliquots  of  growth  medium  were  measured  for  radioactivity  after  24  hrs.  Radioactivity  is  shown 
as  counts  per  minute  (cpm)/ml.  K  and  L.  AA  release  is  evaluated  in  OV202  Shi  and  CI7  cells  and  in  OV2008  NTC, 
ShcP-1  and  -2. 
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performed  arachidonic  acid  (AA)  release  assay  in  OV202Sh1  and  Sh2  cells.  Cells  treated  with  10pM  of 
AACOCF3  showed  significant  reduction  of  AA  release  indicating  that  activity  of  cPLA2  was  inhibited  (Figs.  3J, 
K  and  L). 

Milestone  4:  Molecular  approach.  4a)  Generate  stable  clones  using  inactive  cPLA2a  expression  vector 
mutated  at  Ser  505  at  active  site  (cPLA2S505)  in  OV202Sh  clones. 

Our  several  attempts  in  generating  the  mutant  (cPLA2S505A)  clones  in  OV2Q2Sh1/Sh2  were  not  successful. 

Milestone  5:  Task  3:  Mitochondrial  beta  oxidation  to  assess  the  role  of  LD  to  supply  energy  source  in  in 
OVCA  cells. 

Rationale:  Once  we  established  the 
relationship  between  cPLA2a  and  LD  in  Sh 
clones,  we  next  examined  the  role  of  LD  in  Sh 
clones.  Lipid  droplets  consist  of  a  core  of 
neutral  lipids,  mainly  triacylglycerol  (TAG)  and 
sterol  esters,  that  are  surrounded  by  a 
monolayer  of  polar  lipids  (16).  TAG  is  a 
mechanism  for  storing  unused  calories  and 
free  fatty  acids  (FFA)  can  be  released  from 
Sn-2  position  of  TAG  by  action  of  cPLA2a. 
Resulting  FFA,  depending  upon  the  size,  can  be  metabolized  by  mitochondrial  beta  oxidation  to  generate  ATP 
as  energy  source. 

4a)  To  determine  if  LDs  are  used  as  an  energy  source  under  stress  conditions.  Plate  OV202NTC,  Sh  clones 
and  SKV03  vector  and  CI7  cells  were  exposed  various  stress  conditions  including  serum  starvation,  glucose 
free  media  for  several  days  (1-5  days).  Cell  survival  was  determined  using  MTT  assay  and  number  of  LDs  by 
staining  with  BODIPY. 

OV202NTC,  Shi  and  Sh2  cells  were  grown  in  the  absence  of  serum  (Fig.4a1-A)  and  in  the  presence  of  2DG 
to  inhibit  glucose  uptake  (Fig.4a1-B).  In  the  absence  of  serum  and  or  glucose  uptake  the  Shi  and  Sh2  cells 
underwent  more  cell  death  than  NTC  cells  (Fig.  4a-1-A  and  B). 

OV202Sh1  cells  were  grown  in  EBSS  medium  for  the  indicated 
time  and  the  amount  of  LDs  were  determined  by  BODIPY  staining. 

The  number  of  LDs  increasingly  became  less  and  less  when  cells 
were  grown  in  EBSS  (no  aminoacids)  (Fig.4a-2)  indicating  that  LDs 
may  be  utilized  for  producing  energy  under  these  conditions.  This 
was  tested  in  the  next  task  indicated  below. 


4b)  Linder  similar  experimental  conditions  and  time  points,  we  had 
proposed  to  assess  mitochondrial  beta  oxidation.  6.0pM  [1-1 4C] 
palmitic  acid  in  HBSS  along  with  a-cyclodextrin  will  be  added  for  30  min  to  0V202NTC  and  Sh  batch  clonal 
cells.  Stop  reaction  with  1M  KOH.  Remove  denatured  protein  and  incubate  supernatant  at  600C  for  1  hr, 
neutralize  with  6M  HCI  and  partition.  The  upper  aqueous  phase  will  be  taken  for  measurement  of  [1-14C] 
palmitic  radioactivity.  14C02,  a  radioactive  end  product  of  the  oxidation  reaction 

Instead  of  the  above  method,  we  used  a  non  radioactive  method  as  described  below  to  measure  beta  oxidation 

with  and  without  AACOCF3  treatment  of  NTC  and  Shi  cells.  Using  the 
XF  analyzer  the  oxygen  consumption  rate  (OCR)  reflective  of  oxygen 
consumption  was  used  to  measure  fatty  acid  oxidation  (FAO)  in 
response  to  AACOCF3  treatment  (17). 

Fatty  acid  oxidation:  Oxygen  consumption  rate  was  measured  using  a 
Seahorse  Bioscience  XF24  flux  analyzer.  5x1 04  cells  were  seeded  per 
well  in  triplicates  in  MEM-a  containing  20%  FBS  in  an  XF24  well  culture 
microplates  and  incubated  overnight  in  a  37°C/10%  C02  incubator.  The 

6 


Mitochondrial  |(-oxidation  in  OV202  NTC  and  Shi  cells 
=  50  1  □  NTC  b  Shi 

e 


o 


-50  J  Figure  4b 


Bodipy  staining  of  OV202Sh1  with  EBSS 
0  2  4 


8  16  24hr 


Figure4a-2 


A  <4 


24  hour  starvation 


<o 

< 

o 


OS#fum 

■  No  Sorun 


OV70?  WTC  OV707  Shi 


assay  medium  for  FAO  is  low-buffered  KHB  buffer  (110  mM  NaCI,  4.7  mM  KCI,  2  mM  MgS04  1.2  mM 
Na2HP04,  2.5  mM  glucose  adjusted  to  pH7.4)  supplemented  with  0.5mM  carnitine.  For  induction  of  FAO, 
BSA  conjugated  palmitate  was  injected  to  a  final  concentration  of  50  pM.  XF  analyses  were  performed  in  the 
XF  Extracellular  Flux  Analyzer  (Seahorse  Bioscience,  Billerica,  MA).  Three  basal  rates  were  measured  prior  to 
automated  injection  of  palmitic  acid  (50  pM)  coupled  to  BSA  vehicle  or  BSA  vehicle  alone.  After  treatment  for 
55  min,  the  cPLA2  inhibitor  AACOCF3  (10  pM),  was  added.  Oxygen  consumption  rates  were  measured  by 
using  time-resolved  method  (Seahorse  Bioscience  XF24)  (21).  Data  were  normalized  to  protein  content 
(assayed  after  completion  of  measurements).  Changes  in  the  FAO  induction  in  Shi  cells  are  compared  with 
that  of  NTC  cells.  AACOCF3  induced  inhibition  of  FAO  in  Shi  cells  are  compared  with  the  FAO  inhibition  in 
NTC  (P<  0.01  and  P<  0.001)  are  shown  in  figure  4b.  The  results  clearly  indicate  that  Shi  cells  have  higher 
level  of  beta  oxidation  which  is  inhibited  when  LD  biogenesis  is  inhibited  with  AACOCF3  treatment. 

4c)  Total  levels  of  intracellular  ATP  will  also  be  examined  to  correlate  with  beta  oxidation  in  mitochondria  using 
ATP  kit  from  Promega.  Higher  levels  of  ATP  in  Shi  and  Sh2  compared  to  NTC  cells,  while  this  phenotype  is 
reversed  in  Sh1-CI7  associated  with  increased  glycolysis  in  these  cells  was  reported  in  a  manuscript 
describing  the  role  of  loss  of  HSulf  -1  promoting  increased  glycolysis  Figs.  2E-F  (18)  (  See  manuscript 
attached). 

Milestone  6:  Task  4:  Determine  the  in  vitro  effects  of  cPLA2a  inhibitors  on  migration  and  invasion  of  OVCA 
cells.  Since  we  saw  no  metastatic  invasion  in  our  invivo  model  (See  Figure  5  below)-  we  did  not  pursue  these 

studies  in  vitro. 


Aim  2:  Examine  the  role  of  lipid  droplets  in  chemoresistance  and  metastasis  using  in  vitro  (Milestones 
7-  task  1)  and  in  vivo  models.  Milestones  8  and  9-  tasks  2  and  3 

Taskl:  Establish  relationship  between  LD  and  chemoresistance  in  ovarian  cancer  using  cell  culture  based  study. 
To  examine  the  role  of  cPLA2a-LD  axis  on  chemotherapy-induced  cytotoxicity,  isogenic  cell  lines  with  and  without 
HSulf-1  (OV202NTC,  Shi,  Sh2  and  SKOV3  vec  and  CI7  cells)  will  be  treated  with  optimal  concentration  of  cPLA2 
inhibitor  AACCOF3  alone  (determined  under  experimental  set  2)  with  increasing  cisplatin  concentration.  Parallel 
experiments  will  determine  the  percentage  of  cells  undergoing  apoptosis  under  these  treatment  conditions. 
Chemotherapy-induced  cytotoxicity  will  be  determined  by  trypan  blue  staining  and  clonogenic  survival  assays  (19). 
Data  will  be  expressed  as  the  percentage  of  total  cells  that  are  apoptotic  after  each  treatment. 


We  did  not  pursue  SKOV3  cells  since  they  do  not  express  p-cPLA2  and  very  low  levels  of  total  cPLA2.We  initially 
established  IC50  values  for  carboplatin,  cisplatin  and  AACCOF3  (F3  in  figure  5)  in  OV202NTC,  Shi  cells  by  MTT 
assays  (48hrs)  (Figs.  5A-F).  Using  the  IC50  for  F3,  we  treated  the  Shi  cells  with  increasing  concentrations  of 
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cisplatin  in  vitro  to  determine  if  F3  treatment  would  sensitize  the  cells  to  cisplatin  induced  cytotoxicity.  To  our 
surprise,  we  did  not  observe  a  significant  difference  in  the  inhibition  of  proliferation  in  the  combination  compared  to 
F3  treatment  alone  (Fig.5G).  Interestingly,  we  did  observe  a  time  dependent  decrease  in  the  levels  of  p-cPLA2 
upon  F3  treatment  in  the  Shi  cells  (Fig.5H).  Although,  we  did  not  see  an  effect  of  the  combination  of  platinum  with 
F3  in  vitro,  our  in  vivo  data  shown  in  Figure  6  clearly  showed  increased  TUNEL  staining  and  decrease  in  Ki67 
staining  in  the  F3  and  the  combination  treated  groups  compared  to  untreated  control  or  carboplatin  treated 
groups  alone  (Fig  6H  -J). 
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Task  2:  Examine  the  in  vivo  efficacy  of  cPLA2  inhibitor  in  potentiating  cisplatin  induced  cytotoxicity  in  athymic 

immunodeficient  mice  using  Sh  clones. 

We  had  proposed  to  examine  the  effect  of  cPLA2a  inhibitor  treatment  in  a  mouse  xenograft  model  of  intraperitoneal 
(i.p.)  carcinomatosis  using  Sh  clone  with  anticipation  that  administration  of  cPLA2  inhibitor  will  attenuate  LDs 
biogenesis  in  vivo  which  will  sensitize  tumor  using  OV202Sh1  cells  as  a  ovarian  cancer  model  shown  in  (20).  All 
experimental  use  of  animals  were  complied  with  the  guidelines  of  Animal  Care  at  the  Mayo  Foundation,  in 
accordance  with  approved  protocols.  5X1 06  Sh  clones  expressing  luciferase  cells  in  serum-free  RPMI  1640  will  be 
injected  intraperitoneally  into  female  athymic  nu/nu  mice  at  4  to  5  wk  of  age  (National  Cancer  Institute-Frederick 
Animal  Production  Area  (Frederick,  MD).  Once  i.p  implant  s  are  visible  via  non-invasive  imaging  (about  day  30), 
animals  will  receive  cisplatin  alone  or  in  combination  with  cPLA2a  inhibitor  randomly  assigned  to  three  groups  of  10 
from  day  1  after  i.p.  tumor  detection  via  non-invasive  imaging. 

To  determine  whether  blocking  cPLA2 
activities  via  AACOCF3  confer  any 
advantage  in  attenuating  tumor  growth 
in  combination  with  chemotherapeutic 
agent  such  as  carboplatin,  we 
performed  testing  of  these  agents  in 
vivo  using  a  mouse  model.  The  effect  of 
AACOCF3  alone  and  in  combination 
with  carboplatin  (CBP)  on  primary  tumor 
growth  was  evaluated  in  OV202Sh1 
cells  bearing  nude  mice  i.p. 

(intraperitoneally).  5X1 06  Sh  clones 
expressing  luciferase  cells  in  serum-free 
RPMI  1640  were  injected 
intraperitoneally  into  female  athymic 
nu/nu  mice  at  4  to  5  wk  of  age  (National 
Cancer  Institute-Frederick  Animal 
Production  Area  (Frederick,  MD).  Once 
i.p  implants  were  visible  via  non- 
invasive  imaging  (4  days  after 


Figure  6:  AACOCF3  alone  and  in  combination  with  carboplatin  suppresses  tumor  growth,  and  inhibits 
lipid  droplet  biogenesis  in  vivo 

A.  Randomized  OV202Sh1  tumor-bearing  mice  (n=10),  were  treated  with  water  or  AACOCF3  (F3)  (lOmg/kg),  or 
carboplatin  (CBP)  (51  mg/kg)  or  with  a  combination  of  AACOCF3  and  CBP  for  4  weeks  and  after  28  days  mice  were 
euthanized.  Representative  images  of  the  mice  prior  to  sacrifice  from  control,  F3  treated,  CBP  treated  and  CBP+F3 
combination  treatment  groups  using  IVIS  luminescence  imaging  system  series  2000.  Color  bar  shows  photon  intensity. 

B.  Representative  images  of  excised  tumors  from  OV202  Shi  xenografts.  C.  Excised  tumor  weights  from  control,  F3, 
CBP  and  CBP+F3  combination  treatment  groups.  D.  Abdominal  circumference  from  different  groups  of  mice  measured 
on  the  day  of  sacrifice.  E.  Total  body  weight  of  untreated  control  and  treated  with  F3,  CBP,  CBP+F3.  F.  Immunoblot 
analysis  of  LC3BI/II  in  lysates  from  untreated,  F3,  CBP  and  F3+CBP  treated  Xenografts  with  GAPDH  as  loading  control. 
G.  Tumor  xenografts  were  stained  with  Bodipy  to  detect  lipid  droplets;  DAPI  was  used  to  stain  the  nuclei.  H. 
Representative  images  of  TUNEL  staining  of  frozen  section  of  xenografts  from  different  treatment  groups  and  untreated 
control.  Green  fluorescence  indicates  TUNEL  and  red  fluorescence  indicates  propidium  iodide  (PI)  staining.  I. 
Representative  images  of  immunohistochemical  staining  of  Ki67,  J.  Bar  graph  showing  the  Ki67  proliferation  index 
(*p=<0.05  and  **p=<0.01).  Proliferation  index  (PI)  of  Ki-67  staining  was  measured  using  ImmunoRatio,  from  public 
domain  image  analysis  software.  K.  phospho-cPLA2  in  FFPE  sections  of  treated  and  untreated  xenografts. 

inoculation),  mice  were  randomized  (10  mice/group)  and  treated  i.p.  with  the  lOmg/kg  body  weight  (b.w.)  of 
cPLA2  inhibitor  AACOCF3  (referred  to  as  F3  in  the  figures)  every  other  day  till  the  end  of  the  study,  51  mg/kg 
b.w.  carboplatin  (CBP)  every  5  days  till  the  end  of  the  study  and  a  combination  of  CBP+F3  every  5  days,  with 
vehicle  treated  group  as  control.  Luciferase  imaging  of  the  representative  mice  from  four  groups  is  shown  in 
Fig.6A.  Higher  luciferase  intensity  in  the  control  group  and  CBP  group  indicates  increased  tumor  volume, 
progression  and  metastasis.  Representative  excised  tumor  from  one  of  the  mouse  in  each  group  is  shown  in 
Fig.6B.  At  necropsy,  the  measured  tumor  weight  and  circumference  of  mice  in  each  group  showed  that  the 
combination  treatment  was  more  effective  at  reducing  cancer  progression  compared  to  all  of  the  other  groups 
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(Figs.6C-D).  There  was  no  significant  body  weight  loss  in  F3  or  combination  treated  group  compared  to 
control  group  suggesting  both  F3  and  CBP  treatment  were  well  tolerated  without  apparent  toxicity  to  the  animal 
(Fig.6E).  More  importantly,  Bodipy  staining  of  frozen  sections  of  the  xenografts  clearly  showed  high  levels  of 
lipid  droplets  (Fig.6G)  in  the  control  and  the  CBP  treated  group  and  significantly  less  LDs  in  the  F3  and 
combination  treated  group  consistent  with  the  in  vitro  data  shown  in  figure  2D,  top  panel.  In  contrast,  TUNEL 
staining  showed  increased  TUNEL  staining  in  the  F3  and  the  combination  treated  groups  compared  to 
untreated  control  or  CBP  treated  groups  alone  (Fig.6H)  Anti-tumor  effects  of  AACOCF3  monotherapy  and 
combination  with  CBP  correlated  with  significant  reductions  in  tumor  cell  proliferation  marker,  Ki67  and 
phospho-cPLA2  as  determined  by  IHC  staining  using  the  specific  antibodies  to  these  proteins  (Figs.  61  and  K). 


Aim  3:  Milestone  10;  Determine  the  expression  of  phospho-cPLA2a  in  primary  ovarian  tumors  and 
correlate  with  clinical  parameters. 

To  determine  if  cPLA2a  is  a  potential  target  in  ovarian  cancer,  we  had  proposed  to  evaluate  the  expression  of 

p-cPLA2a  s505-in  200  ovarian 

cancer  specimens  by 

immunohistochemistry  and 

correlate  to  clinical  outcome  such 
as  progression  free  survival,  overall 
survival,  grade,  stage  and  histology 
and  also  determine  if  p-cPLAas505 
predict  responsiveness  to  adjuvant 
chemotherapy  in  ovarian  cancer. 
Although  we  had  proposed  to 
evaluate  the  expression  of  p- 
cPLA2a  s505-in  200  ovarian 

cancer  specimens,  we  determined 
p-cPLA2  expression  in  512  ovarian 
tumors  by  IHC.  These  tumors  are 
all  from  patients  before  they  commenced  standard  chemotherapy.  We  initially  established  the  conditions  for 
IHC  in  the  Pathology  Research  Core  at  the  Mayo  Clinic.  Immunohistochemistry  was  performed  by  the  Mayo 
Clinic  Pathology  core  on  5  different  TMAs  with  512  tumors  in  triplicates.  The  images  were  digitized  and  the 
expression  levels  were  blindly  read  by  three  different  investigators  including  a  pathologist.  The  scores  were 
averaged  and  the  expression  level  was  correlated  with  clinical  outcome.  To  our  surprise,  we  saw  that  patient 
tumors  with  high  level  of  p-cPLA2  had  a  better  progression  free  survival  and  overall  survival  compared  to 
patient  tumors  with  low  levels  of  p-cPLA2. 


We  also  saw  an  inverse  correlation  in  the  expression  of  HSuf-1  and  p-cPLA2  in  a  subset  of  tumors  that 
expressed  either  of  these  proteins  in  the  same  tumors.  163  out  of  512  tumors  expressed  both  Hsulf-1  and  p- 
cPLA2.  72  of  163  (44%)  of  these  tumors  showed  an  inverse  correlation.  However,  this  inverse  correlation  in 
expression  shown  in  Fig.  7  below  did  not  correlate  with  clinical  response. 


Milestone  11-  Manuscript  Preparation 

Please  see  appendix  for  the  two  manuscripts  where  we  have  acknowledged  the  current  support  from  the  DOD. 
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4.  KEY  RESEARCH  ACCOMPLISHMENTS:  We  have  determined  that 


•  Using  two  different  isogenic  cell  lines  OV202  NTC  vs  Shi  and  OV2008  cells  with  high  levels  of  p- 
cPLA2  vs  OV2008ShcP2  cells  with  genetic  downregulation  of  total  cPLA2,  we  have  determined  that 
increased  ERK  activity  in  OVCA  cells  results  in  increased  cPLA2  activity.  These  results  show  the 
contribution  of  loss  of  HSulf-1  with  high  pERK  leading  to  LD  biogenesis,  also  a  direct  contribution  of  p- 
cPLA2  in  promoting  LD  biogenesis  in  ovarian  cancer. 

•  We  have  also  shown  using  specific  inhibitors  of  p-cPLA2,  AACOCF3  and  MAFP  that  LD  biogenesis  is 
attenuated.  LD  inhibition  by  these  inhibitors  also  sensitizes  these  cells  to  cisplatin  mediated  cytotoxicity 
implicating  LDs  as  a  disease  marker  that  can  be  therapeutically  targeted.  Additionally,  we  have 
determined  that  downregulation  of  p-cPLA2  sensitizes  OV2008  cells  to  cisplatin  mediated  cytotoxicity. 
Taken  together,  these  data  indicate  that  high  levels  of  p-cPLA2  may  promote  cisplatin  resistance. 

•  Our  in  vivo  data  using  a  specific  cPLA2  inhibitor,  AACOCF3  showed  less  tumor  burden  compared  to 
the  untreated  control  indicating  that  a  higher  level  of  cPLA2  itself  can  play  a  crucial  role  in  ovarian 
tumorigenesis.  Interestingly,  AACOCF3  and  carboplatin  combination  therapy  was  more  effective  to 
reduce  tumor  growth  than  carboplatin  alone.  This  is  a  first  report  to  show  that  HSulf-1  loss  associated 
ovarian  tumor  formation  is  regulated  by  higher  activity  of  cPLA2.  More  studies  are  required  to  test  the 
synergistic  effect,  absorption  and  bioavailability  of  AACOCF3  alone  versus  the  combination  with 
carboplatin  in  the  mouse  model.  Our  findings  indicate  that  cPLA2  inhibition  is  a  potent  therapeutic 
approach  to  synergize  with  conventional  chemotherapies  to  reduce  ovarian  tumor  growth. 

•  More  importantly,  we  have  determined  the  expression  level  of  p-cPLA2  in  more  than  500  ovarian 
tumors  on  a  tissue  microarray  by  immunohistochemistry.  These  tumors  are  all  from  patients  before 
they  commenced  standard  chemotherapy.  To  our  surprise,  we  saw  that  patient  tumors  with  high  level 
of  p-cPLA2  had  a  better  progression  free  survival  and  overall  survival  compared  to  patient  tumors  with 
low  levels  of  p-cPLA2. 

Conclusion:  HSulf-1  and  its  tumor  suppressor  activities  are  long  been  a  research  interest  of  our  group.  We 
have  established  the  protective  role  of  HSulf-1  in  the  regulation  of  growth  factor  mediated  tumor  growth  in 
ovarian  cancer.  Recently,  we  have  shown  that  HSulf-1  loss  also  connects  to  metabolic  alterations  which  play  a 
critical  role  for  tumor  growth  (18,  20).  Among  several  metabolites,  lipid  and  glycolytic  metabolites  were  the  one 
which  altered  most  significantly.  When  we  studied  the  specific  changes  among  the  lipid  metabolites,  we 
observed  significant  phenotypical  changes  in  HSulf-1  knockout  cell  lines,  i.e.  lipid  droplets  (LDs).  LDs  are 
unique  organelles  which  primarily  act  as  storage  of  neutral  lipids  and  cholesterol  and  recent  findings  from 
various  investigators  connected  its  role  with  disease  progression  including  cancer  and  metabolic  diseases.  To 
our  surprise,  we  have  observed  that  LDs  are  associated  with  various  proteins  which  are  all  enhanced  in  their 
mRNA  and  protein  level  upon  loss  of  HSulf-1.  Further,  we  targeted  one  specific  LD  associated  protein,  cPLA2 
and  used  its  functional  inhibitors  to  therapeutically  target  the  cells  with  loss  of  HSulf-1  and  higher  LD  content. 
We  have  shown  that,  inhibition  of  cPLA2  activity  not  only  reduces  LDs  but  also  reduces  cell  proliferation, 
release  of  free  fatty  acid,  inhibition  of  beta  oxidation,  and  most  importantly  attenuates  the  tumor  growth  in  vivo. 
These  data  support  our  hypothesis  that  lipid  droplet  can  be  used  as  a  disease  marker  in  ovarian  cancer.  In 
other  cells  with  loss  of  HSulf-1,  other  lipogenic  gene  may  be  playing  a  role  such  as  FASN,  SREBP1C  and 
PLA2G3  playing  a  role  in  promoting  lipid  droplet  biogenesis  or  a  lipogenic  phenotype.  This  is  not  surprising 
since  loss  of  HSulf-1  in  ovarian  cancer  affects  many  genes  involved  in  the  lipogenic  pathway  as  we  have 
shown  in  our  recent  study  (1).  Targeting  the  lipogenic  phenotype  in  ovarian  cancer  with  targeted  inhibitors  of 
the  pathway  or  other  broad  acting  inhibitors  that  can  reduce  the  number  of  lipid  droplets  may  have  a  beneficial 
effect  on  patients  with  this  deadly  disease.  Although  obesity  perse  has  not  been  linked  to  ovarian  cancer,  our 
data  suggest  that  targeting  the  lipogenic  pathway  may  be  of  interest  in  the  future  trials  for  ovarian  cancer.  The 
results  from  this  study  will  be  equally  applicable  to  all  women  with  ovarian  cancer,  both  civil  and  women 
serving  in  the  army  and  their  relatives  afflicted  with  ovarian  cancer. 

Future  Directions 

We  have  extended  our  understanding  on  the  role  of  HSulf-1  in  the  regulation  of  autophagy  and  lipid 
metabolism  that  has  opened  up  newer  research  avenues  to  investigate  the  impact  of  tumor  suppressor 
activities  of  HSulf-1  in  regard  to  autophagic  regulation  od  LDs. 
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6.  PUBLICATIONS,  ABSTRACTS,  AND  PRESENTATIONS: 

(1)  Lay  Press:  Nothing  to  report 

(2)  Peer-Reviewed  Scientific  Journals:  All  attached  in  the  Appendix 

•  Roy,  Mondal  and  Shridhar  et -Loss  of  HSulf-1  promotes  altered  lipid  metabolism  in  ovarian  cancer 
Cancer  and  Metabolism  IN  PRESS,  2014. 

•  Mondal  et.al.,  -  HSulf-1  deficiency  dictates  a  metabolic  reprograming  of  glycolysis  and  TCA  cycle  in 
ovarian  cancer-  Oncotarget,  2015 

•  Manuscript  under  revision  -  Roy  et  al-  Loss  of  HSulf-1 -The  Missing  Link  between  Autophagy  and  Lipid 
Droplets  in  Ovarian  Cancer- Autophagy,  2016 


(3)  Invited  Articles:  Nothing  to  report 

a.  Abstracts:  All  five  abstracts  attached  in  the  appendix  section 

1.  Debarshi  Roy,  Xiaoping  He,  Ashwani  Khurana,  Susmita  Mondal,  Deok-Beom  Jung, 

Sung  Kim,  Thomas  Dierks,  Clifford  Folmes,  Andre  Terzic  and  Viji  Shridhar.  HSulf-1  loss 
in  ovarian  cancer  cells  induces  lipid  droplet  biogenesis  105th  AACR  Annual  Meeting  San 
Diego,  CA,  April  5-9,  2014.Abstract  Number:  1293. 

2.  Debarshi  Roy,  Chen  Wang,  Shailendra  Giri,  Xiaoping  He,  Ashwani  Khurana  ,  Susmita 
Mondal,  Deok-Beom  Jung,  Han  Ihn,  Sung  Kim'  Shaun  Riska,  Ann  Oberg,  Jeremy  Chien 
and  Viji  Shridhar.  Loss  of  HSulf-1  Promotes  Altered  Cellular  Metabolism  of  Ovarian 
Cancer.  The  17  International  Conference  on  Korean  Medicine.  Cancer  Preventive 
Material  Development  Research  Center  (CPMDRC)  &  Institute,  Kyunghee  University, 
2013. 

3.  Viji  Shridhar.  The  Missing  Link  between  Lipid  Droplets  and  Autophagy  in  Ovarian 
Cancer.  1st  annual  Meeting  International  Ovarian  Cancer  Consortium  Tumor 
Microenvironment  and  Drug  Discovery,  Oklahoma  City,  OK. 

4.  Viji  Shridhar.  The  Missing  Link  between  Lipid  Droplets  and  Autophagy  in  Ovarian 
Cancer.  GCC  Global  Cancer  Conference  &  Medicare  summit.  September  15-17,  2014 
Hyderabad  International  Convention  Centre,  India. 

5.  Debarshi  Roy,  Susmita  Mondal,  Ashwani  Khurana,  Xiaoping  He,  Edward  Hammond, 
Keith  Dredge,  and  Viji  Shridhar.  Loss  of  HSulf-1  Promotes  Defective  Autophagy  and 
Increased  Lipid  Droplet  Biogenesis  in  vitro  and  in  vivo  in  Ovarian  Cancer.  American 
Association  of  Cancer  Research,  May  18-22,  2015,  Philadelphia,  PA.  ABSTRACT  #- 
2413. 


b.  Presentations:  Invited  Presentations-Local,  National  and  International 

(1)  2014-09/15-  Title:  The  Missing  Link  between  Lipid  Droplets  and  Autophagy  in  Ovarian 
Cancer.  Global  Cancer  Conference  -  Hyderabad,  India. 

(2)  2014-02/07  Title:  The  Missing  Link  between  Lipid  Droplets  and  Autophagy  in  Ovarian 
Cancer.  1st  annual  Meeting  International  Ovarian  Cancer  Consortium  Tumor 
Microenvironment  and  Drug  Discovery,  Oklahoma  City,  OK. 

(3)  2014-01/28  Title:  The  Link  between  lipid  droplets  and  autophagy  in  Ovarian  Cancer- 
Division  of  Experimental  Pathology,  Mayo  Clinic,  Rochester,  MN. 

(4)  2013-11/14  Title:  Identification  of  HSulf-1  as  a  regulator  of  altered  cellular  metabolism 
in  ovarian  cancer*.  Gyn/Onc  Cancer  Seminar  series,  Mayo  Clinic,  Rochester,  MN. 
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(5)  2013-11/26  Title:  Loss  HSulf-1  promotes  altered  Cellular  Metabolism  of  Ovarian 
Cancer*.  2013,  the  17  International  Conference  on  Korean  Medicine.  Cancer  Preventive 
Material  Development  Research  Center  (CPMDRC)  &  Institute,  Kyunghee  University. 

(6)  2013-11/07  Title:  Story  of  HSulf-1:  From  Genomics  to  Glycomics  to  Metabolomics*. 
University  of  Louisville,  James  Graham  Brown  Cancer  Center,  KY. 

(7)  2013-06:  Title:  Metabolic  Alterations  Mediated  by  Loss  of  HSulf-1  in  Ovarian  Cancer*, 
University  of  Kansas  Medical  Center,  Kansas  City,  Kansas. 

7.  INVENTIONS,  PATENTS  AND  LICENSES:  Nothing  to  Report. 

8.  REPORTABLE  OUTCOMES:  Provide  a  list  of  reportable  outcomes  that  have  resulted  from  this 
research.  Reportable  outcomes  are  defined  as  a  research  result  that  is  or  relates  to  a  product, 
scientific  advance,  or  research  tool  that  makes  a  meaningful  contribution  toward  the  understanding, 
prevention,  diagnosis,  prognosis,  treatment  and  /or  rehabilitation  of  a  disease,  injury  or  condition,  or 
to  improve  the  quality  of  life.  This  list  may  include  development  of  prototypes,  computer  programs 
and/or  software  (such  as  databases  and  animal  models,  etc.)  or  similar  products  that  may  be 
commercialized.  Nothing  to  report 

9.  OTHER  ACHIEVEMENTS: 

a.  We  have  developed  isogenic  non-targeted  control  transduced  OV2008  NTC  cell  line  and  a  batch 
clonal  line  of  ShRNA  downregulated  cPLA2  cell  line  (OV2008ShcP2). 

b. 

c.  Funding  applied  for  based  on  work  supported  by  this  award-  LOI  to  OCRF  postdoctoral  award  by 
Dr.  Roy  in  2013-  But  was  not  invited  back  to  apply  for  the  full  proposal. 


10. REFERENCES: 
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11.  APPENDICES:  Five  Abstracts  and  Three  Manuscripts*  attached 
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Abstract  1:  105th  AACR  Annual  Meeting  San  Diego,  CA,  April  5-9,  2014.Abstract  Number:1293. 

HSulf-1  loss  in  ovarian  cancer  cells  induces  lipid  droplet  biogenesis 

Debarshi  Roy,  Xiaoping  He,  Ashwani  Khurana,  Susmita  Mondal,  Deok-Beom  Jung,  Sung  Kim,  Thomas 

Dierks,  Clifford  Folmes,  Andre  Terzic,  Viji  Shridhar. 

We  have  previously  shown  that  a  frequent  loss  of  endosulfatase  HSulf-1  in  ovarian  cancer  enhances  the 
heparin  binding  growth  factor  signaling  which  subsequently  induces  and  potentiates  tumorigenesis.  Our  recent 
data  indicate  that  knockdown  of  HSulf-1  in  ovarian  cancer  cells  significantly  alter  levels  of  several  metabolites 
which  is  comparable  to  the  deranged  metabolite  profiles  in  primary  ovarian  tumors.  Here  we  report  that  the 
metabolites  associated  with  lipid  metabolism  are  altered  in  such  a  way  that  the  cells  start  manifesting  a 
“lipogenic  phenotype”  characterized  by  an  increased  fatty  acid  synthesis,  storage  and  mitochondrial  beta 
oxidation.  The  lipogenic  phenotype  is  associated  with  an  accumulation  of  cytoplasmic  lipid  droplets  (LDs)  in  a 
stably  HSulf-1  knockdown  cells  [OV202  (Shl/2)  and  OV2223  cells].  Presence  of  LDs  is  evidenced  by  BODIPY 
staining  and  transmission  electron  microscopy  (TEM)  analysis  in  Shl/2  and  a  non-transduced  control  (NTC) 
cells.  Conversely,  enhanced  expression  of  HSulf-1  in  TOV21G,  SKOV3,  and  Cl 3  cells  also  resulted  in 
reducing  the  number  of  LDs.  More  importantly,  MEFs  (mouse  embryonic  fibroblasts)  from  HSulf-1  knockout 
mice  also  showed  an  increased  LD  biogenesis  compared  to  the  wild  type  MEFs.  Microarray  analysis  identified 
cytoplasmic  phospholipase  A2  (cPLA2;  an  enzyme  involved  in  LD  biogenesis)  as  one  of  the  genes  that  was 
significantly  upregulated  in  OV202  Shi  and  Sh2  cells.  We  found  that  Shi  and  Sh2  cells  show  an  increased 
phosphorylation  (ser505)  of  cytoplasmic  phospholipase  A2  (p-cPLA2)  compared  to  NTC  cells.  Phosphorylation 
at  ser505  is  vital  for  the  enzyme  activity  of  cPLA2  and  also  necessary  for  LD  biogenesis.  We  inhibited  the  LD 
biogenesis  in  Shi  cells  by  a  pharmacological  inhibitor  of  cPLA2  (AACOCF3  and  MAFP)  as  well  by  shRNA 
mediated  knockdown  of  endogenous  cPLA2.  Further,  we  stimulated  NTC  and  Shi  cells  with  heparin  binding 
growth  factor  (50ng/ml)  for  different  time  points  and  found  an  elevated  phosphorylation  of  cPLA2  and 
extracellular  signal  related  kinase  (ERK),  which  is  also  inhibited  by  a  pre-treatment  of  cells  with  cPLA2  inhibitor. 
Interestingly,  inhibition  of  phospho-ERK  by  U0126  (MEK  inhibitor)  also  inhibits  the  phosphorylation  of  cPLA2. 
These  findings  strongly  suggest  that  loss  of  HSulf-1  promotes  lipid  metabolism  in  the  ovarian  cancer  cells  via 
activation  of  cPLA2  via  an  ERK  dependent  pathway.  The  present  work  highlights  a  novel  role  HSulf-1  in 
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regulating  lipid  metabolism  which  could  lead  to  a  better  understanding  of  tumorigenesis  and  identification  of 
new  targets  for  tumor  therapy. 


Abstract  2  :  presented  at  the  17  International  Conference  on  Korean  Medicine.  Cancer  Preventive  Material 
Development  Research  Center  (CPMDRC)  &  Institute,  Kyunghee  University,  11/26/2013. 

Loss  of  HSulf-1  Promotes  Altered  Cellular  Metabolism  of  Ovarian  Cancer 
Debarshi  Roy,  Chen  Wang,  Shailendra  Giri,  Xiaoping  He,  Ashwani  Khurana  ,  Susmita  Mondal,  Deok-Beom 
Jung,  Han  Ihn,  Sung  Kim'  Shaun  Riska,  Ann  Oberg,  Jeremy  Chien  and  Viji  Shridhar* 

Loss  of  the  endosulfatase  HSulf-1  is  common  in  ovarian  cancer,  upregulates  heparin  binding  growth  factor 
signaling  and  potentiates  tumorigenesis  and  angiogenesis.  Since  growth  factor  signaling  is  closely  tied  to 
metabolic  signaling,  we  determined  the  extent  to  which  HSulf-1  loss  affects  cancer  cell  metabolism.  We 
performed  untargeted  global  metabolite  profiling  in  the  HSulf-1  ShRNA  knockdown  cells  compared  to  isogenic  non- 
targeted  control  transduced  cells  and  identified  ~  338  metabolites  using  gas  chromatography  and  mass 
spectrometry  (GC/MS)  and  liquid  chromatography  tandem  mass  spectrometry  (LC/MS/MS  platforms). 
Knockdown  of  HSulf-1  in  OV202  cells  led  to  significant  changes  in  156  metabolites  pertaining  to  several 
metabolic  pathways  including  amino  acid,  glycolysis,  lipids,  glutathione  and  extracellular  matrix  remodeling. 
Combining  ingenuity  pathway  analysis  of  the  gene  expression  changes  of  these  isogenic  cell  lines  and  Kyoto 
Encyclopedia  of  Genes  and  Genomics  (KEGG)  database  analysis  of  commonly  altered  metabolites  indicated 
that,  specifically  changes  in  the  metabolites  associated  with  lipid  metabolism  leading  to  a  “lipogenic 
phenotype”  characterized  by  an  increase  in  fatty  acid  synthesis,  storage  and  mitochondrial  beta  oxidation. 
Evaluation  of  the  elevated  lipid  metabolites  upon  HSulf-1  knockdown  in  ovarian  cancer  cells  indicated  that 
these  were  similar  to  the  metabolites  reported  to  be  elevated  in  primary  ovarian  tumors,  thereby  alluding  to 
the  possibility  that  loss  of  HSulf-1  could  potentially  contribute  to  the  metabolic  alterations  associated  with  the 
progression  of  ovarian  pathogenesis,  specifically  impacting  the  lipogenic  phenotype  of  ovarian  cancer  cells. 


Abstract  3:  1st  annual  Meeting  International  Ovarian  Cancer  Consortium  Tumor  Microenvironment  and  Drug 
Discovery,  Oklahoma  City,  OK,  02/07/2014 

The  Missing  Link  between  Lipid  Droplets  and  Autophagy  in  Ovarian  Cancer 

Viji  Shridhar,  Ph.D. 

Department  of  Experimental  Pathology  and  Laboratory  Medicine,  Mayo  Clinic,  Rochester,  MN. 

Excess  free  fatty  acids  are  stored  as  triglycerides  in  lipid  droplets  (LD)  to  be  used  for  energy  under  nutrient 
deprived  conditions.  A  second  cellular  response  to  starvation  is  autophagy,  in  which  the  cell  digests  its  own 
components  to  provide  nutrients.  Stable  downregulation  of  HSulf-1  in  OV202  (Shl/2)  and  OV2223  lead  to 
increased  LDs  as  evidenced  by  BODIPY  staining  and  transmission  electron  microscopy  (TEM)  analysis 
compared  to  non-transduced  control  (NTC)  cells.  Conversely,  enhanced  expression  of  HSulf-1  in  TOV21G, 
and  SKOV3  cells  resulted  in  reducing  the  number  of  LDs  compared  to  vector  transfected  controls.  More 
importantly,  MEFs  from  HSulf-1  KO  mice  also  showed  increased  LD  biogenesis  compared  the  WT  MEFs. 
Microarray  analysis  identified  cPLA2,  a  protein  involved  in  LD  biogenesis  as  one  of  the  genes  that  was 
significantly  upregulated  OV202  Shi  and  Sh2  cells.  Pharmacological  inhibition  and  ShRNA  downregulation  of 
cPLA2  activity  resulted  in  decreasing  the  #  of  LDs  in  OV202Sh1  cells.  TEM  analysis  showed  that  all  HSulf-1 
proficient  cells  displayed  increased  number  of  autophagic  vesicles  compared  to  isogenic  HSulf-1  deficient 
cells.  Pharmacological  inhibition  of  cPLA2  promoted  autophagy  in  OV202  Shi  cells  and  elevated  DAPK  and 
LC3BII  levels  suggesting  a  novel  cross-talk  between  LD  biogenesis  and  autophagy  in  ovarian  cancer  cells. 
Enhanced  expression  of  WT  DAPK  reduced  LD  formation  and  induced  autophagy  related  markers  such  as 
LC3BII,  Beclin,  ATG12  and  5  in  OV202Sh1  cells.  HS  mimetic  PG545,  a  tumor  microenvironment  targeting 
drug  that  mimics  the  action  of  HSulf-1,  reduced  LDs  and  promoted  autophagy  in  Shi  and  Sh2  cells. 
Collectively,  these  results  identify  the  critical  role  of  HSulf-1,  a  major  regulator  of  growth  factor  mediated 
signaling  in  the  tumor  microenvironment  as  the  missing  link  in  regulating  both  autophagy  and  LDs  in  ovarian 
cancer. 
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Abstract 

Dysregulation  of  autophagy  and  altered  metabolic  pathways  are  frequently  observed  in  cancer.  Due  to  these 
alterations,  pharmacological  targeting  of  these  two  pathways  simultaneously  could  provide  a  viable  therapeutic 
option.  Although  the  association  between  these  two  pathways  is  well  characterized  in  metabolic  disorders,  it  is 
not  well  defined  in  ovarian  cancer  (OVCA).  In  this  regard,  we  found  that  loss  of  endosulfatase  HSulf-1,  a 
known  putative  tumor  suppressor,  suppresses  LC3-GFP  foci  formation  and  promotes  increased  lipid  droplet 
(LD)  biogenesis  suggesting  that  absence  of  HSulf-1  in  OVCA  affects  both  autophagy  and  lipid 
metabolism.  While  isogenic  cells  with  genetic  ablation  of  HSulf-1  (OV202Sh1/2  and  TOV2223Sh1  cells) 
displayed  LDs,  the  nontargeted  control  transduced  (NTC)  OV202  and  TOV2223  cells  had  significantly  less 
LDs.  In  contrast,  Transmission  Electron  Micrographs  (TEMs)  showed  that  OV202  and  TOV2223  NTC  cells 
had  significantly  more  autophagic  vacuoles  (AVs)  compared  to  their  isogenic  ShRNA  targeted  cells. 
Conversely,  ectopic  expression  of  HSulf-1  in  SKOV3  cells  decreased  the  number  of  LDs  and  increased  the 
number  of  AVs  compared  to  vector  transfected  controls.  Here  we  report  that  OV202Sh1  cells  and  HSulf-1 
deficient  OV2008  cells  have  increased  p-cPLA2a(ser505)  levels  that  are  associated  with  biogenesis  of  large 
number  of  LDs  with  reduced  AVs.  Interestingly,  pharmacological  inhibition  of  cPLA2a  with  AACOCF3  in 
OV202Sh1  cells  resulted  in  reduced  LD  biogenesis,  inhibited  colony  formation  and  reduced  tumorigenesis  in 
vivo.  More  importantly,  treatment  of  HSulf-1  deficient  cells  with  HS  mimetic  PG545  which  can  compensate  for 
loss  of  HSulf-1,  reduced  LD  biogenesis,  promoted  autophagy  and  inhibited  tumor  growth  in  vivo.  Collectively, 
these  results  show  a  critical  role  of  HSulf-1  in  regulating  both  autophagy  and  LD  biogenesis  in  ovarian  cancer. 
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It  is  increasingly  being  recognized  that  altered  lipid  metabolism  is  an  early  event  in  carcinogenesis  and  a 
central  hallmark  of  many  cancers.  Under  nutrient  deprived  conditions,  excess  free  fatty  acids  are  stored  as 
triglycerides  in  lipid  droplets  (LD)  to  be  used  for  energy.  A  second  cellular  response  to  starvation  is  autophagy, 
in  which  the  cell  digests  its  own  components  to  provide  nutrients.  Microarray  and  Metabolomics  profiling 
identified  the  lipid  pathway  as  one  the  major  pathways  modulated  by  loss  of  HSulf-1  in  ovarian  cancer  (OVCA). 
HSulf-1  deficient  cells  (OV202  Sh1/sh2,  OV2223  and  Sulf-1  KO  MEFs)  possess  high  levels  of  lipid  droplets 
(LD)  that  are  absent  in  the  HSulf-1  proficient  isogenic  cells  (TOV21G  and  SKOV3).  More  importantly,  TEM 
analysis  showed  that  all  HSulf-1  proficient  cells  displayed  increased  number  of  autophagic  vesicles  compared 
to  isogenic  HSulf-1  deficient  cells.  Pharmacological  inhibition  and  ShRNA  downregulation  of  cPLA2  activity,  a 
protein  involved  in  LD  biogenesis  resulted  in  decreasing  the  #  of  LDs  and  promoted  autophagy  in  OV202 
Shlcells  and  elevated  the  autophagy  inducing  protein  DAPK  and  autophagy  related  markers  including  LC3BII 
levels  suggesting  a  novel  cross-talk  between  LD  biogenesis  and  autophagy  in  ovarian  cancer  cells.  HS 
mimetic  PG545,  a  tumor  microenvironment  targeting  drug  that  mimics  the  action  of  HSulf-1,  reduced  LDs  and 
promoted  autophagy  in  Shi  and  Sh2  cells.  Collectively,  these  results  identify  the  critical  role  of  HSulf-1,  a 
major  regulator  of  growth  factor  mediated  signaling  in  the  tumor  microenvironment  as  the  missing  link  in 
regulating  both  autophagy  and  LDs  in  OVCA. 
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Abstract 

Background:  Loss  of  the  endosulfatase  HSulf-1  is  common  in  ovarian  cancer,  upregulates  heparin  binding  growth 
factor  signaling  and  potentiates  tumorigenesis  and  angiogenesis.  However,  metabolic  differences  between  isogenic 
cells  with  and  without  HSulf-1  have  not  been  characterized  upon  HSulf-1  suppression  in  vitro.  Since  growth  factor 
signaling  is  closely  tied  to  metabolic  alterations,  we  determined  the  extent  to  which  HSulf-1  loss  affects  cancer  cell 
metabolism. 

Results:  Ingenuity  pathway  analysis  of  gene  expression  in  HSulf-1  shRNA-silenced  cells  (Shi  and  Sh2  cells)  compared 
to  non-targeted  control  shRNA  cells  (NTC  cells)  and  subsequent  Kyoto  Encyclopedia  of  Genes  and  Genomics  (KEGG) 
database  analysis  showed  altered  metabolic  pathways  with  changes  in  the  lipid  metabolism  as  one  of  the  major 
pathways  altered  inShl  and  2  cells.  Untargeted  global  metabolomic  profiling  in  these  isogenic  cell  lines  identified 
approximately  338  metabolites  using  GC/MS  and  LC/MS/MS  platforms.  Knockdown  of  HSulf-1  in  OV202  cells  induced 
significant  changes  in  156  metabolites  associated  with  several  metabolic  pathways  including  amino  acid,  lipids,  and 
nucleotides.  Loss  of  HSulf-1  promoted  overall  fatty  acid  synthesis  leading  to  enhance  the  metabolite  levels  of  long 
chain,  branched,  and  essential  fatty  acids  along  with  sphingolipids.  Furthermore,  HSulf-1  loss  induced  the 
expression  of  lipogenic  genes  including  FASN,  SREBF1,  PPARy,  and  PLA2G3  stimulated  lipid  droplet  accumulation. 
Conversely,  re-expression  of  HSulf-1  in  Shi  cells  reduced  the  lipid  droplet  formation.  Additionally,  HSulf-1  also  enhanced 
CPT1 A  and  fatty  acid  oxidation  and  augmented  the  protein  expression  of  key  lipolytic  enzymes  such  as  MAGL,  DAGLA, 
HSL,  and  ASCL1 .  Overall,  these  findings  suggest  that  loss  of  HSulf-1  by  concomitantly  enhancing  fatty  acid  synthesis 
and  oxidation  confers  a  lipogenic  phenotype  leading  to  the  metabolic  alterations  associated  with  the  progression  of 
ovarian  cancer. 

Conclusions:  Taken  together,  these  findings  demonstrate  that  loss  of  HSulf-1  potentially  contributes  to  the  metabolic 
alterations  associated  with  the  progression  of  ovarian  pathogenesis,  specifically  impacting  the  lipogenic  phenotype  of 
ovarian  cancer  cells  that  can  be  therapeutically  targeted. 

Keywords:  HSulf-1,  Lipogenesis,  Lipolysis,  Lipid  droplets,  Microarray  and  metabolite  profiling 
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Background 

Metabolic  re-programming  has  recently  emerged  as  a  new 
hallmark  of  cancer.  Alteration  of  cellular  metabolism  in 
cancer  cells  is  proposed  to  increase  the  availability  of  es¬ 
sential  building  blocks  that  support  uncontrolled  cellular 
proliferation  [1].  Most  cancer  cells,  although  diversified  by 
etiology  and  type,  reprogram  their  metabolism  to  accumu¬ 
late  metabolic  intermediates  as  sources  of  building  blocks 
[2].  The  Warburg  effect  is  one  of  the  most  important  meta¬ 
bolic  alteration  in  cancer,  in  which  neoplastic  cells  exhibit 
higher  glucose  uptake  and  utilization  by  altering  glucose 
metabolism  even  in  the  presence  of  oxygen  to  produce  lac¬ 
tate  from  glucose  and  thereby  decouples  glycolysis  from 
mitochondrial  oxidation  [3].  Although  alterations  in 
fatty  acids  (FAs)  and  lipid  metabolism  have  received  less 
attention,  recently,  their  importance  in  cancer  metabolism 
is  being  increasingly  recognized.  The  total  lipid  pools  re¬ 
quired  for  membrane  synthesis  of  dividing  cells  are  derived 
mainly  from  FAs  and  in  part  from  acetyl  CoA  [4].  Cancer 
cells  meet  their  FAs  demand  mainly  by  increasing  de  novo 
FA  synthesis  rather  than  from  exogenous  sources  and 
partly  from  acetyl  CoA  [5,6] .  The  activated  FAs  can  then  be 
utilized  in  the  synthesis  of  membrane  phospholipids  (PLs) 
including  phosphatidylcholine  (PC),  phosphatidylethanol- 
amine  (PE)  in  addition  to  sterols,  sphingolipids,  and  lysoli- 
pids  to  meet  the  energy  demands  and  proliferation.  The 
rest  of  the  activated  FAs  then  function  as  signaling  mole¬ 
cules  or  esterified  as  triglycerols  or  sterol  esters  and  stored 
in  lipid  droplets  (LDs)  [7,8]. 

Accumulating  evidence  suggest  that  activation  of  onco¬ 
genes  such  as  MYC,  NF-kB,  K-RAS  [2,9]  and  loss  of  tumor 
suppressor  genes  (P.53,  LKB1/AMPK)  [10-12]  forms  a 
basis  for  altered  metabolism  of  cancer  cells.  The  onco¬ 
genic  activation  of  c-MYC  turns  on  the  Ras-Raf-MAPK 
signaling  pathway  along  with  HIFla  and  PI3K-Akt-mTOR 
axis  which  transcriptionally  stimulate  the  expression  of 
most  glycolytic  and  glutaminolytic  genes  and  subsequently 
activate  lipid  metabolism  [13].  Similarly,  loss  of  tumor 
suppressors  forms  the  basis  of  the  Warburg  effect  leading 
to  carcinogenesis.  For  example,  p53  activates  TIGAR  to 
reduce  the  cellular  accumulation  of  fructose  2,6  bispho- 
sphate,  an  allosteric  activator  of  phosphofructo  kinase,  a 
critical  control  point  in  glycolysis  leads  to  decreased  gly¬ 
colysis  [10].  Moreover,  growth  factor- mediated  phosphor¬ 
ylation  of  pyruvate  kinase  isoform  2  (PKM2)  [14]  and 
mutations  of  genes  in  IDH1/2  [15]  in  the  metabolic  path¬ 
ways  such  as  glycolysis  and  the  TCA  cycle  respectively 
have  been  identified  that  support  growth  of  the  rapidly 
proliferating  cells  and  survive  metabolic  stress  [16,17]. 

Human  Sulfatase  1  (HSulf-1),  an  endosulfatase  estab¬ 
lished  as  a  putative  tumor  suppressor  in  ovarian  cancer,  has 
been  shown  to  modulate  the  signaling  of  growth  factors 
and  cytokines  in  tumor  microenvironment  [18].  Our  previ¬ 
ous  work  demonstrated  that  loss  of  HSulf-1  modulates 


heparin-binding  growth  factors  such  as  bFGF,  VEGF,  HGF, 
PDGF,  and  heparin  binding  EGF  (HB-EGF)  signaling, 
which  plays  an  important  role  in  tumor  progression,  me¬ 
tastasis,  and  angiogenesis  [19-22].  Moreover,  serous  tu¬ 
mors  with  moderate  to  high  levels  of  HSulf-1  had  better 
prognosis  in  terms  of  overall  survival,  implicating  its  crit¬ 
ical  role  in  the  progression  of  ovarian  cancer  [23].  Our 
more  recent  data  demonstrated  that  HSulf-1  knockdown 
clones  in  the  OV202  ovarian  cancer  cell  line  (OV202Shl 
and  Sh2  cells)  have  significantly  increased  the  ability  to 
form  anchorage-independent  colonies  in  vitro  and  en¬ 
hanced  tumorigenicity  in  vivo  [24].  Also,  in  breast  cancer, 
HSulf-1  is  negatively  regulated  by  HIFla,  but  positively  by 
von  Hippel-Lindau  tumor  suppressor  gene  [23,25,26]. 
These  findings  led  us  to  hypothesize  that  HSulf-1  might 
play  a  unique  role  to  alter  tumor  microenvironment  rais¬ 
ing  the  possibility  that  its  loss  might  alter  the  cellular  me¬ 
tabolism  and  levels  of  the  resulting  metabolites  as  a 
downstream  effect  of  altered  growth  factor  signaling. 

In  the  current  study,  global  changes  in  metabolism  were 
investigated  in  HSulf-1  silenced  OV202  cells  by  microarray, 
metabolic  data  analysis,  and  Western  blotting.  Here,  we  for 
the  first  time  report  that  loss  of  HSulf-1  promotes  overall 
fatty  acid  synthesis  and  oxidation  leading  to  a  lipogenic 
phenotype  to  promote  cancer  growth  in  ovarian  cancer. 

Methods 
Cell  culture 

OV202  cell  line  was  low-passage  primary  line  established  at 
the  Mayo  Clinic  [27].  OV202NTC,  Shi,  Sh2,  and  Cl  11  cells 
were  cultured  in  5%  CC>2-95%  air  humidified  atmosphere  at 
37°C  with  minimal  essential  medium  supplemented  with 
20%  fetal  bovine  serum  and  1  pg/ml  puromycin,  with  non- 
essential  amino  acids.  All  cell  lines  were  tested  using  a 
PCR-based  assay  and  found  to  be  free  of  Mycoplasma 
contamination. 

ShRNA 

HSulf-1  short-hairpin  (sh)  RNA1  (Shi-  AGCTACCCT 
GGGTTCCTTTGT)  which  targets  the  3 '-untranslated 
region  (UTR)  was  cloned  into  lentiviral  vector  pLKO.l- 
puro  as  described  previously  [27].  HSulf-1  shRNA2  (Sh2- 
CGTCGAATTTGAAGGTGAAAT)  and  nontargeted 
control  shRNAs  (NTC  shRNA-  ACTTACGAGTGAC 
AGTAGATT)  cloned  into  the  lentivirus  vector  pLKO.l- 
puro  were  chosen  from  the  human  library  (MISSION 
TRC-Hs  1.0)  and  purchased  as  glycerol  stock  from  Sigma. 
Transfection  with  Fugene  (Roche)  was  performed  accord¬ 
ing  to  the  manufacturer's  instructions.  Transducted  cells 
were  selected  with  1  pg/ml  puromycin. 

Rescue  of  HSulf-1  in  Shi  cells 

pcDNA-HSulf-1  plasmid  was  cloned  as  described  previ¬ 
ously  [19].  Since  Shi  shRNA  targeted  the  3 'UTR  of  HSulf- 
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1,  we  rescued  the  expression  of  HSulf-1  in  this  cell  line 
with  CMV-driven  WT  expression  construct  and  selected 
stable  clone  Cll  as  previously  described  [19].  Vector  only 
transfected  cells  served  as  controls. 

Microarray  expression  data  analysis 

OV202  NTC,  Shi,  and  Sh2  cells  in  triplicates  were  profiled 
using  Illumina  Human  HT-12  3.0  Expression  Beadchip 
array  as  previously  described  [28].  Microarray  expression 
data  were  analyzed  on  the  log2  scale.  Data  quality  was 
assessed  via  box  and  whisker  plots  along  with  residual  and 
pair-wise  MVA  plots  before  and  after  normalization 
[29,30].  All  arrays  were  normalized  together  using  fastlo,  a 
non-linear  normalization  similar  to  cyclic  loess  which  runs 
in  a  fraction  of  the  time  [31].  Both  supervised  and  un¬ 
supervised  analyses  were  performed.  Supervised  analysis 
to  determine  differentially  expressed  genes  was  performed 
using  Significance  Analysis  of  Microarrays  (SAM)  [32]. 
For  SAM  analysis  for  unpaired  samples,  Biometric  Re¬ 
search  Branch  (BRB)-ArrayTools  (Version  3.7.0,  devel¬ 
oped  by  Dr.  Richard  Simon  and  Amy  Peng  Lam.)  was 
used  with  Delta  set  to  0.822,  resulting  in  false  discovery 
rate  <5%.  Unsupervised  clustering  was  performed  using 
the  one  minus  correlation  metric  with  average  linkage. 
Heat  maps  were  generated  for  visualization.  Pathway  ana¬ 
lysis  was  performed  using  Ingenuity  Pathway  Analysis  (In¬ 
genuity8  Systems,  www.ingenuity.com). 

Liquid  chromatography/mass  spectrometry 
(LC/MS,  LC/MS2) 

The  LC/MS  portion  of  the  platform  was  based  on  a  Waters 
ACQUIT Y  UPLC  and  a  Thermo-Finnigan  LTQ  mass  spec¬ 
trometer,  which  consisted  of  an  electrospray  ionization 
(ESI)  source  and  linear  ion-trap  (LIT)  mass  analyzer.  The 
sample  extract  was  split  into  2  aliquots,  dried,  then  recon¬ 
stituted  in  acidic  or  basic  LC-compatible  solvents,  each  of 
which  contained  11  or  more  injection  standards  at  fixed 
concentrations.  One  aliquot  was  analyzed  using  acidic 
positive  ion  optimized  conditions  and  the  other  using  basic 
negative  ion  optimized  conditions  in  two  independent  in¬ 
jections  using  separate  dedicated  columns.  Extracts  recon¬ 
stituted  in  acidic  conditions  were  gradient  eluted  using 
water  and  methanol  both  containing  0.1%  Formic  acid, 
while  the  basic  extracts,  which  also  used  water/methanol, 
contained  6.5  mM  ammonium  bicarbonate.  The  MS  ana¬ 
lysis  alternated  between  MS  and  data-dependent  MS2 
scans  using  dynamic  exclusion. 

Gas  chromatography/mass  spectrometry  (GC/MS) 

The  samples  destined  for  GC/MS  analysis  were  re-dried 
under  vacuum  desiccation  for  a  minimum  of  24  h  prior  to 
being  derivatized  under  dried  nitrogen  using  bistrimethyl- 
silyl-triflouroacetamide  (BSTFA).  The  GC  column  was  5% 
phenyldimethyl  silicone  and  the  temperature  ramp  is  from 


40°C  to  300°C  in  a  16-min  period.  Samples  were  analyzed 
by  a  Thermo-Finnigan  Trace  DSQ  fast-scanning  single- 
quadrupole  mass  spectrometer  using  electron  impact 
ionization.  The  instrument  was  tuned  and  calibrated  for 
mass  resolution  and  mass  accuracy  on  a  daily  basis.  The 
information  output  from  the  raw  data  files  was  automatic¬ 
ally  extracted  as  discussed  below. 

Data  extraction  and  compound  identification 

Peaks  were  identified  using  Metabolon's  proprietary  peak 
integration  software.  Compounds  were  identified  by  com¬ 
parison  to  library  entries  of  purified  standards  or  recurrent 
unknown  entities.  Identification  of  known  chemical  en¬ 
tities  was  based  on  comparison  to  metabolomic  library  en¬ 
tries  of  purified  standards.  As  of  this  writing,  more  than 
2,600  commercially  available  purified  standard  com¬ 
pounds  had  been  identified  and  registered  into  LIMS  for 
distribution  to  both  the  LC  and  GC  platforms  for  deter¬ 
mination  of  their  analytical  characteristics.  The  combin¬ 
ation  of  chromatographic  properties  and  mass  spectra 
gave  an  indication  of  a  match  to  the  specific  compound  or 
an  isobaric  entity.  Metabolon  data  analysts  use  proprietary 
visualization  and  interpretation  software  to  confirm  the 
consistency  of  peak  identification  among  the  various  sam¬ 
ples.  Library  matches  for  each  compound  were  checked 
for  each  sample  and  corrected  if  necessary. 

Normalization 

Raw  data  from  each  sample  was  normalized  to  protein 
concentration  as  measured  by  Bradford  assay  prior  to 
statistical  analysis. 

Sample  accessioning 

Each  sample  received  was  accessioned  into  the  Metabolon 
LIMS  system  and  was  assigned  by  the  LIMS,  a  unique 
identifier,  which  was  associated  with  the  original  source 
identifier  only.  This  identifier  was  used  to  track  all  sample 
handling,  tasks,  and  results.  The  samples  (and  all  derived 
aliquots)  were  bar-coded  and  tracked  by  the  LIMS  system. 
All  portions  of  any  sample  were  automatically  assigned 
their  own  unique  identifiers  by  the  LIMS  when  a  new  task 
was  created;  the  relationship  of  these  samples  was  also 
tracked.  All  samples  were  maintained  at  -80°C  until 
processed. 

Sample  preparation 

The  sample  preparation  process  was  carried  out  using 
the  automated  MicroLab  STAR8  system  from  Hamilton 
Company  (Reno,  NV,  USA).  Recovery  standards  were 
added  prior  to  the  first  step  in  the  extraction  process  for 
QC  purposes.  Sample  preparation  was  conducted  using 
a  proprietary  series  of  organic  and  aqueous  extractions  to 
remove  the  protein  fraction  while  allowing  maximum  re¬ 
covery  of  small  molecules.  The  resulting  extract  was 
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divided  into  four  fractions;  two  for  analysis  by  LC  and  one 
for  analysis  by  GC  and  a  forth  as  a  spare.  Samples  were 
placed  briefly  on  a  Turbo Vap'  (Zymark,  Hopkinton,  MA, 
USA)  to  remove  the  organic  solvent.  Each  sample  was  then 
frozen  and  dried  under  vacuum.  Samples  were  then  pre¬ 
pared  for  the  appropriate  instrument,  either  LC/MS  or  GC/ 
MS. 

QA/QC 

For  QA/QC  purposes,  a  number  of  additional  samples  are 
included  with  each  day's  analysis.  Furthermore,  a  selection 
of  QC  compounds  is  added  to  every  sample,  including 
those  under  test.  These  compounds  are  carefully  chosen 
so  as  not  to  interfere  with  the  measurement  of  the  en¬ 
dogenous  compounds.  These  QC  samples  are  primarily 
used  to  evaluate  the  process  control  for  each  study  as  well 
as  aiding  in  the  data  curation. 

Metaboanalyst 

Differently  expressed  metabolites  between  Shl/Sh2  and 
baseline  conditions  were  firstly  mapped  to  KEGG  me¬ 
tabolites  IDs  according  to  Human  Metabolome  Database 
(HMDB;  URL:  http://www.hmdb.ca/)  [33].  Then,  pathway 
analysis  was  performed  to  highlight  relevant  metabolic 
pathways  defined  in  KEGG  database  (http://www.genome. 
jp/kegg/),  using  an  on-line  tool  named  MetaboAnalyst 
(http://www.metaboanalyst.ca/)  [34,35].  Specifically,  two 
types  of  pathway  analysis  were  done:  one  is  over¬ 
representation  analysis  using  hypergeometric  test  [35], 
asking  if  differentially  expressed  metabolites  are  par¬ 
ticularly  enriched  in  a  same  pathway;  the  other  is  path¬ 
way  topology  analysis  summarizing  relative-betweeness 
centrality  [34],  investigating  potential  pathway  impact 
of  observed  metabolite  changes  based  on  known  path¬ 
way  topology  relationships. 

Western  blot  analysis 

Western  blot  analysis  was  performed  as  described  previ¬ 
ously  [36].  Whole  cell  lysates  were  analyzed  with  the  fol¬ 
lowing  antibodies:  FASN,  ASCL1  (Cell  signaling),  SREBPlc, 
PLA2G3,  HSulf-1  (Abeam,  AB96533),  CPT1A,  HSL, 
DAGLA,  [3-tubulin  (GeneTex)  and  (3-actin  (Sigma-Aldrich). 

Real-time  PCR 

Quantitative  real-time  PCR  (qRT-PCR)  was  carried  out 
using  SYBR-Green  PCR  Master  Mix  (Applied  Biosys¬ 
tems,  Foster  City,  CA,  USA),  with  specific  primers  for 
the  genes  shown  in  this  study.  GAPDH  or  18S  ribosomal 
subunit  (Applied  Biosystems)  were  used  as  internal  control 
in  a  Light  Cycler  kit  (BioRad  Chromo  4).  Normalization 
across  samples  was  carried  out  using  the  average  of  the 
constitutive  human  gene  18S  and/or  GAPDH  primers  and 
calculated  as  previously  described  [18].  Binding  efficiencies 


of  primer  sets  for  both  target  and  reference  genes  were 
similar. 

Bodipy  staining 

Cells  (50,000)  were  seeded  on  a  coverslip  in  a  24-well  plate 
and  were  grown  for  24  hours  in  the  presence  of  complete 
growth  medium.  Cells  were  washed  and  fixed  in  4%  para¬ 
formaldehyde  for  10  min  at  room  temperature  before 
staining  with  1  pg/ml  BODIPY  (493/503;  Sigma,  St.  Louis, 
MI,  USA)  in  PBS  for  10  min  at  room  temperature.  Cover- 
slips  were  washed  with  PBS  and  mounted  in  a  slide  with 
Prolong  Gold  Antifade  Reagent  (Invitrogen).  BODIPY 
stained  cells  were  examined  under  inverted  confocal  fluor¬ 
escence  microscope  (Zeiss). 

Transient  transfection 

To  determine  the  effect  of  PLA2G3  on  lipid  droplet  bio¬ 
genesis  OV202  NTC  cells  were  transiently  transfected 
with  plasmids  containing  empty  vector  or  cDNA  encoding 
PLA2G3.  After  24  h  of  transfection,  we  performed  BOD¬ 
IPY  staining  to  visualize  lipid  droplets.  PLA2G3  plasmid 
was  obtained  on  a  MTA  from  Addgene. 

Fatty  acid  synthesis 

Cells  were  washed  twice  in  cold  PBS  and  resuspend  in 
lysis  buffer  (50  mM  Tris-HCl,  pH  7.4,  1  mM  EDTA,  150 
mM  NaCl  and  PMSF.  Cells  were  sonicated  and  homoge¬ 
nized  by  dounce  homogenizer  followed  by  centrifugation 
at  13,000  rpm  for  15  min  at  4°C.  The  supernatant  was  col¬ 
lected  and  proteins  were  measured  by  Bradford  assay,  and 
100  pg  of  protein  was  used  to  conduct  FASN  activity 
assay.  FASN  activity  was  measured  by  protocol  described 
by  Vazquez-Martin  et  al.  [37].  Briefly,  100  pg  of  protein 
was  incubated  with  240  pM  NADPH,  30  pM  acetyl  CoA 
and  50  pM  malonyl  CoA  in  assay  buffer  (200  mM  potas¬ 
sium  phosphate,  pH  6.6,  1  mM  DTT,  1  mM  EDTA)  and 
oxidation  of  NADPH  was  measured  by  monitoring  the  ab¬ 
sorbance  at  340  nM  over  the  period  of  the  time.  Results 
presented  here  compares  FASN  activity  monitored  for  10 
min.  Values  are  presented  as  nanomolar  NADPH  oxidized 
per  minute  per  milligram  of  protein. 

Fatty  acid  oxidation 

Oxygen  consumption  rate  was  measured  using  a  Seahorse 
Bioscience  XF24  flux  analyzer.  5  x  104  cells  were  seeded 
per  well  in  triplicates  in  MEM-k  containing  20%  FBS  in 
an  XF24  well  culture  microplates  and  incubated  overnight 
in  a  37°C/10%  C02  incubator.  The  assay  medium  for  FAO 
is  low-buffered  I<HB  buffer  (110  mM  NaCl,  4.7  mM  I<C1, 
2  mM  MgS04,  1.2  mM  Na2HP04,  2.5  mM  glucose  ad¬ 
justed  to  pH  7.4)  supplemented  with  0.5  mM  carnitine. 
For  induction  of  FAO,  BSA  conjugated  palmitate  was 
injected  to  a  final  concentration  of  50  pM.  XF  analyses 
were  performed  in  the  XF  Extracellular  Flux  Analyzer 
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(Seahorse  Bioscience,  Billerica,  MA,  USA).  Three  basal 
rates  were  measured  prior  to  automated  injection  of  pal¬ 
mitic  acid  (50  pM)  coupled  to  BSA  vehicle  or  BSA  vehicle 
alone.  After  treatment  for  55  min,  the  carnitine  palmitoyl 
transferase-1  inhibitor,  Etomoxir  (ETO,  50  pM),  was 
added.  Oxygen  consumption  rates  were  measured  by 
using  time-resolved  method  (Seahorse  Bioscience  XF24) 
(21).  Data  were  normalized  to  protein  content  (assayed 
after  completion  of  measurements). 

Proliferation  assay 

Equal  number  of  cells  (1  x  10s)  was  plated  in  triplicate  in 
12-well  plates.  OV202NTC,  Shi,  and  Sh2  cells  were 
counted  after  24,  48,  and  72  h  using  a  cellometer  (Nexelom, 
Lawrence,  MA,  USA).  For  Etomoxir  treatment,  Equal  num¬ 
ber  of  cells  (1  x  105)  were  seeded  in  12-well  plates  in  tripli¬ 
cate  and  treated  with  increasing  concentration  of  Etomoxir 
(0  to  100  pM)  for  24  h  and  total  cell  numbers  were  counted 
using  cellometer. 

Results 

Loss  of  HSulf-1  comprehensively  altered  major  metabolic 
pathways 

We  recently  reported  that  HSulf-1  knockdown  clones 
in  the  OV202  ovarian  cancer  cell  line  (OV202Shl  and 
Sh2  cells)  have  significantly  increased  ability  to  form 
anchorage-independent  colonies  in  vitro  and  enhanced 
tumorigenicity  in  vivo  [24] .  Consistent  with  these  observa¬ 
tions,  our  growth  assays  showed  enhanced  growth  rate  in 
Shi  and  Sh2  cells  compared  to  NTC  cells  (Additional 
file  1:  Figure  SI).  Here,  to  elucidate  the  function  of  puta¬ 
tive  tumor  suppressor  HSulf-1  in  the  metabolism  of  ovar¬ 
ian  cancer,  we  performed  gene  expression  profiling  of 
stably  knockdown  HSulf-1  clonal  lines  OV202  Shi  and 
Sh2  cells  (referred  to  from  hereon  as  Shi  and  Sh2)  com¬ 
pared  to  HSulf-1  expressing  non-targeted  control  cells 
(OV202NTC,  referred  to  as  NTC)  [26]  in  triplicates  using 
Illumina  HumanHT  12  v3  platform  [28].  Unbiased  hier¬ 
archical  clustering  and  heat  maps  showed  that  genes  in 
several  different  pathways  were  differentially  expressed  in 
Shi  and  Sh2  compared  to  NTC  cells  (Figure  1A).  We 
found  that  over  1,645  and  780  genes  were  differentially 
expressed  in  Shi  and  Sh2  cells,  respectively,  compared  to 
NTC  at  2.6  FC  (p  and  FDR  <0.0001).  We  identified  500 
and  280  altered  genes  in  Shi  and  Sh2  cells  respectively 
by  significance  analysis  of  microarrays  (SAM)  [32]  from 
the  comprehensive  list  of  2,752  genes  which  encoded  all 
known  human  metabolic  enzymes  and  transporters  re¬ 
ported  by  Possemato  et  al.  [38].  Ingenuity  pathway 
analysis  (http://www.ingenuity.com)  for  these  genes 
showed  that  most  genes  were  differentially  regulated  in 
the  fatty  acid/lipid  pathways  in  Shi  and  Sh2  cells  com¬ 
pared  to  NTC  cells  (Additional  file  2:  Table  SI,  Figure  IB). 
We  next  explored  the  271  genes  in  the  lipid  related 


pathways  from  our  microarray  data  by  unsupervised  clus¬ 
tering  (Figure  1C)  and  found  that  26%  (73  of  271)  genes 
were  differentially  expressed  in  Shi  and  Sh2  cells  com¬ 
pared  to  NTC  cells  (Figure  ID),  indicating  that  tumor  sup¬ 
pressor  HSulf-1  possibly  regulates  the  lipid  metabolism  in 
ovarian  cancer  cells. 

Loss  of  HSulf-1  altered  global  metabolic  profile  in  ovarian 
cancer  cells 

To  determine  whether  loss  of  HSulf-1  has  effect  in  cellu¬ 
lar  metabolism,  we  performed  unbiased  global  metabolic 
profiling  using  the  Metabolon  platform  (Metabolon  Inc, 
Durham,  NC,  USA)  in  Shi  and  Sh2  cells  compared  to 
NTC.  The  samples  were  extracted  using  Metabolon’s 
standard  solvent  extraction  method  from  cells  in  loga¬ 
rithmic  phase  with  5  biological  replicates  for  each  sam¬ 
ple  and  distributed  into  equal  parts  for  analysis  on  the 
GC/MS  and  LC/MS/MS  platforms.  Hierarchical  cluster¬ 
ing  revealed  that  metabolite  levels  of  Shi  and  Sh2  cells 
cluster  together  separately  from  control  NTC  (Figure  2A). 
Our  initial  principal  component  analysis  (PCA)  [39]  re¬ 
vealed  that  Shi  and  Sh2  cells  had  specific  group  of  metab¬ 
olites  which  were  different  from  NTC  cells  (Figure  2B). 
Additionally,  t  test  in  each  principal  component  dimen¬ 
sion  showed  that  PCA  differences  between  Shl-Sh2  and 
NTC  cells  mainly  reside  in  the  first  and  second  principal 
components,  with  statistically  significant  p  values  of  1.8e 
-  10  and  4.5e  -  6,  respectively.  Consistent  with  this,  two- 
dimensional  PCA  plots  showed  differences  between  Shl- 
Sh2  and  NTC  are  mainly  in  the  first  and  second  PC  di¬ 
mensions  (Figure  2C,D,E).  We  also  performed  linear  re¬ 
gression  analysis  to  evaluate  association  of  each  principal 
component  vs.  major  metabolite  classes,  to  determine  if 
lipid  class  is  one  of  the  dominant  factors  determining 
the  first  and  second  principal  components.  The  associ¬ 
ation  was  evaluated  in  a  multivariate  regression  model 
(Additional  file  3:  Figure  S2),  where  coefficients  of  first/ 
second  PCs  were  treated  as  dependent  variable,  and 
metabolite  class  labels  were  regarded  as  independent 
variables.  This  analysis  showed  that  the  lipid  class  has 
statistically  significant  association  with  both  first  and 
second  PC  dimensions,  with  p  values  =  5.0e  -  3  and 
1.4e  -  3,  respectively.  Interestingly,  the  peptide  class, 
despite  its  relatively  small  size,  also  has  significant  asso¬ 
ciations  with  both  first  and  second  PC  dimensions, 
while  the  amino  acid  class  has  significant  association 
with  the  first  principal  component. 

Moreover,  this  analysis  showed  a  total  of  338  known 
metabolites  altered  by  the  loss  of  HSulf-1  in  Shi  and 
Sh2  cells.  Among  them,  Shi  and  Sh2  cells  had  a  total  of 
193  and  188  biochemical,  respectively,  which  were  sig¬ 
nificantly  altered  {p  <  0.05,  Welch’s  t  test)  compared  to 
NTC  cells  (Figure  2F,  Additional  file  4:  Table  S2,  and 
Additional  file  5:  Table  S3).  Additionally,  with  the  loss  of 
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Figure  1  (See  legend  on  next  page.) 
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(See  figure  on  previous  page.) 

Figure  1  Microarray  analysis  of  differentially  expressed  genes  in  OV202NTC,  Shi,  and  Sh2  cells.  (A)  Unsupervised  hierarchical  clustering 
normalized  expression  values  of  20,090  selected  probe  sets  for  NTC,  Shi,  and  Sh2.  Each  class  is  represented  by  three  biological  replicates.  Red: 
expression  values  above  the  average  across  all  samples;  blue:  expression  values  below  the  average  across  all  samples.  (B)  Ingenuity  pathway  analysis  of 
metabolic  genes.  The  most  statistically  significant  metabolic  pathways  identified  in  the  confirmed  cell-specific  marker  list  are  listed  according  to  their  p 
value  (-Log)  ( blue  bars)  and  the  ratio  of  list  genes  found  in  each  pathway  over  the  total  number  of  genes  in  that  pathway  (Ratio,  orange  squares).  The 
threshold  line  corresponds  to  a  pathway  enrichment  p  value  of  0.05.  (C)  Unsupervised  hierarchical  clustering  normalized  expression  values  of  271 
selected  lipid  pathway  related  probe  sets  for  NTC,  Shi,  and  Sh2.  Each  class  is  represented  by  three  biological  replicates.  Red:  expression  values  above 
the  average  across  all  samples;  blue:  expression  values  below  the  average  across  all  samples.  (D)  The  most  differentially  expressed  lipid  pathway  related 
genes  (73  of  the  271  genes  in  A,  FDR  =  <1%).  The  statistical  difference  was  tested  according  to  t  test  using  'genefiltefi  package  in  R,  and  multi 
hypothesis-testing  corrected  FDR  was  estimated  using  'fdrtool'  package  in  R.  The  black  arrow  indicates  genes  analyzed  by  real-time  and/or  western  blot 
analysis  in  this  study. 


n 


Type 

■  AmmoacxJ 

■  Carbohydrate 

■  Cofactors  and  vitamns 

■  Energy 

■  Up»d 
Nudeotide 

■  Peptide 
Xenobiotic$ 


F 


Statistical  Comparisons 
Welch's  Two-Sample  t-Test 

Significantly 

Altered 

Biochemicals 

Shi 

NTC 

Sh2 

NTC 

Total 

biochemical* 

0SO.O5 

193 

188 

Biochemical* 

(li) 

64| 129 

74|114 

Total 

biochemical* 

0.05<p<0.10 

28 

19 

Biochemical* 

(li) 

10|  18 

9|  10 

c«-» ) 


G 


Others  Lipids 

20%  28% 


Figure  2  Metabolite  profile  of  ovarian  cell  line  expressing  HSulf-1  (NTC)  and  HSulf-1  downregulates  Shi  and  Sh2  cells.  (A)  Heatmap 
showing  338  biochemicals  in  lysates  from  5  replicates  each  of  OV202  cells  expressing  HSulf-1  (NTC),  HSulf-1  Shi  and  Sh2  cells  arranged  by  unsupervised 
clustering.  (B)  Three-dimensional  sample  PCA  plot  of  log-transformed  normalized  concentration  of  338  biochemicals,  where  samples  were  colored  by  NTC, 
Shi,  and  Sh2  cell-types.  (C,  D,  E)  2-D  PCA  plots:  two-dimensional  PCA  plots  between  paired  principal  components  (first  PC  vs.  second  PC;  second  PC  vs. 
third  PC;  first  PC  vs.  third  PC).  x-/y-axis  corresponds  to  loadings  of  specified  PC  dimension  and  each  point  is  a  metabolite  sample,  with  black,  green,  and  red 
colors  indicating  NTC,  Shi,  and  Sh2  conditions,  respectively.  (F)  Table  showing  significantly  altered  biochemical  in  Shi  and  Sh2  cells  compared  to  NTC  cells. 
Welch's  two-sample  f  test  was  used  to  identify  biochemicals  that  differed  significantly  between  experimental  groups.  A  summary  of  the  numbers  of 
biochemicals  that  achieved  statistical  significance  (p  <  0.05),  as  well  as  those  approaching  significance  (0.05  <p<  0.10),  are  shown  along  with  the  Venn 
diagram  of  shared  metabolites  between  Shi  and  Sh2  cells.  (G)  Pie  chart  of  percentage  of  common  altered  metabolites  in  Shi  and  Sh2  cells  in 
major  pathways. 
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HSulf-1,  129  and  114  metabolites  were  downregulated, 
whereas  64  and  74  metabolites  were  upregulated  in  Shi 
and  Sh2  cells,  respectively.  Also,  40%  of  these  biochemi¬ 
cals  were  altered  in  the  same  direction  in  Shi  and  Sh2 
cells  compared  to  NTC  cells  (p  =  <0.05,  with  false  discov¬ 
ery  rate  (FDR)  <0.05,  Additional  file  6:  Table  S4).  Interest¬ 
ingly,  we  found  that  among  all  the  metabolites  altered, 
28%  were  lipids,  22%  were  amino  acids,  whereas  sugar 
and  dipeptide  comprised  7%  (Figure  2G),  implicating  that 
loss  of  HSulf-1  mediates  a  global  metabolic  alteration  in 
ovarian  cancer  with  changes  in  the  lipid  class  being  a 
major  contributor. 

However,  some  metabolites  were  differentially  altered  in 
the  Shi  and  Sh2  cells.  The  differentially  changed  metabo¬ 
lites  between  these  ShRNAs  could  also  be  a  function  of  the 
level  of  knockdown  and  where  the  ShRNAs  could  have  in¬ 
tegrated.  The  extent  of  knockdown  in  Shi  targeting  the 
3'UTR  was  close  to  100%.  However,  with  Sh2  RNA  target¬ 
ing  the  open  reading  frame,  there  was  still  some  level  of 
HSulf-1  present  (Figure  1A,  [24]). 

Alteration  of  lipid  metabolites  upon  HSulf-1  loss 

The  lipogenic  phenotype  characterized  by  the  activation 
of  lipid  metabolism  is  recognized  as  a  universal  feature 
of  most  cancers  [40,41].  Apart  from  the  fatty  acid  (FA) 
uptake,  cancer  cells  requires  de  novo  FA  biosynthesis 
to  synthesize  new  membranes,  to  store  energy  in  lipid 
droplets  and  to  form  the  lipidic  platform  for  signaling 
in  membrane  level  in  lipid  rafts  for  increased  signaling 
of  cell  growth  receptors  [42,43].  Moreover,  circulating 
lipids  also  play  a  significant  role  in  cancer  cell  growth, 
migration,  and  invasion  [44,45].  Of  the  156  common 
metabolites  altered  in  Shi  and  Sh2  cells,  44  (28%)  me¬ 
tabolites  belonged  to  the  lipid  class  including  long- 
chain  FAs,  lysolipids,  sphingolipids,  glycerolipids,  ei- 
cosanoids,  and  carnitine  (Additional  file  6:  Table  S4), 
in  which  20  metabolites  were  upregulated  (45%).  Using 
the  Kyoto  Encyclopedia  of  Genes  and  Genomics  (KEGG) 
database,  we  mapped  these  metabolites  to  major  pathways 
impacted  with  alterations  in  the  key  junctions  [34,35]. 
These  included  linoleic  acid,  glycerophospholipid,  arachi- 
donic  acid  (Additional  file  7:  Figure  S3),  and  sphingolipid 
pathways  [46,47].  Consistent  with  these  data,  we  found 
that  among  the  20  long-chain  fatty  acid  metabolites  de¬ 
tected,  13  metabolites  were  significantly  increased  by  2-6- 
fold  including  palmitate,  stearate,  and  oleate,  while  the 
remaining  7  long-chain  FAs  did  not  show  any  significant 
differences  in  both  Shi  and  Sh2  cells  (Figure  3A).  Of  note, 
docosadienoate  (22:2n6),  10-nonadeconoate  and  eicoseno- 
ate  (20:ln9)  levels  were  increased  approximately  sevenfold 
compared  to  NTC.  Additionally,  all  the  detected  branched 
fatty  acid  metabolites  were  augmented  (Figure  3B)  while 
five  of  six  essential  fatty  acids  detected  were  increased 
2-5.5-fold  (Figure  3C). 


The  key  structural  lipids  in  cell  membranes  are  the 
glycerol-phospholipids  including  phosphatidyl-choline, 
phosphatidyl-ethanolamine,  phosphatidyl-serine,  phos¬ 
phatidyl-inositol,  and  phosphatidic  acid  in  addition  to 
other  lipids,  such  as  sterols,  sphingolipids,  and  lyso- 
phospholipids.  Sphingolipids  such  as  ceramide,  sphingo- 
sine,  and  sphingosine-l-phosphate  are  bioactive  lipids 
which  can  dictate  the  signaling  including  growth  factor 
responses,  inflammation,  apoptosis,  and  proliferation 
[48].  Our  metabolite  analysis  revealed  that  HSulf-1 
deficiency  was  closely  associated  with  increased  levels 
of  sphinganine  (FC  =  21.79,  p  =  <0.001),  sphingosine 
(FC  =  5.08,  p  =  <0.001),  palmitoyl  sphingomyelin  (FC  =  1.62, 
p  =  <0.0047,  and  FC=1.52,  p  =  <0.001),  and  stearoyl 
sphingomyelin  (FC  =  6.24,  p  =  <0.001  and  FC  =  5.30, 
p  =  <0.001)  (Figure  3D).  Additionally,  we  found  that 
knockdown  of  HSulf-1  expression  also  largely  affected  the 
metabolite  levels  of  choline/inositol  pathway  (Figure  3E) 
and  lysophospholipids  (Figure  3F)  including  reduced 
glycerophosphorylcholine  (GPC)  (FC  =  -1.5 ,  p  =  0.0940  in 
Shi  cells  only)  and  glycerol  3-phosphate  (G3P)  (FC  =  -1.6, 
p- 0.0154  and  FC  =  -1.9,  p  =  <0.001).  Altogether,  these 
results  implicated  a  key  role  of  HSulf-1  in  increased  lipid 
metabolism  and  signaling. 

Deficiency  of  HSulf-1  in  ovarian  cancer  induced  higher 
expression  of  'lipogenic  genes' 

All  the  lipid  molecules  in  cells  are  derived  in  part  from 
acetyl  CoA,  and  many  contain  FAs.  These  FA  building 
blocks  come  from  either  exogenous  sources  or  from  de 
novo  FA  synthesis.  Thus,  malignant  cells  synthesize 
their  own  FA  de  novo  and  thereby  exhibit  a  preference 
over  exogenous  FA  uptake,  while  most  normal  human 
cells  prefer  exogenous  sources  [49].  Our  microarray 
analysis  showed  that  mRNA  levels  of  lipogenic  en¬ 
zymes  fatty  acid  synthase  ( FASN ),  sphingosine  kinase  1 
(SPHK1),  PLA2G4A,  PLA2G3,  sterol  regulatory  element¬ 
binding  transcription  factor  1  ( SREBF1 ),  and  peroxisome 
proliferator-activated  receptor  ( PPARy )  are  also  upregu¬ 
lated  upon  loss  of  HSulf-1  (Figure  4A).  Consistent  with 
the  increased  levels  of  these  genes  and  lipid  metabolites, 
qRT  PCR  and  immunoblotting  showed  enhanced  mRNA 
(Figure  4B)  and  protein  (Figure  4C)  expressions  of  FASN, 
SPHK1,  PLA2G4A,  PLA2G3,  SREBF1,  and  PPARy  in  the 
Shi  and  Sh2  cells  compared  to  the  NTC  cells.  The  in¬ 
crease  in  the  mRNA  and  protein  expression  of  FASN 
along  with  the  increased  production  of  several  long  chain 
FAs  in  Shi  and  Sh2  cells  indicated  that  fatty  acid  synthesis 
was  enhanced  in  HSulf-1  silenced  cells.  To  confirm,  we 
also  measured  the  enzymatic  activity  of  FASN  and  found 
almost  twofold  higher  activity  of  FASN  in  Shi  and  Sh2 
cells  compared  to  NTC  cells  (Figure  4D).  These  data  sug¬ 
gest  that  loss  of  HSulf-1  can  increase  the  activity  of  FASN 
to  enhance  FA  synthesis. 


Roy  ef  al.  Cancer  &  Metabolism  201 4,  2:1 3 
http://www.cancerandmetabolism.eom/content/2/1/13 


Page  9  of  1 7 


Essential  Fatty  Acids 


6.00 

5.00 

4.00 

3.00 

S?  2.00 
■ 

JS 

u  loo 
®  0.00 
-1.00 
-2.00 
-3.00 
-4.00 


|  Shl/NTC  □  SM2/NTC 


*!*  j  j 

I  *** 

*  *1!  f| 

i  I  I  I 


1  z 


II 


5  S 


Choline/Inositol  pathway 

]  Sh2/NTC 


\ 


iBTyirF 


I 


1 1 
E  * 
*  r! 
S  " 

|| 

II 


1? 

If 


J 


I 

I 


I  I 


Figure  3  (See  legend  on  next  page.) 


D 


Sphingolipids 


7.00 

6.00 

gj  5.00 

jS  4.00 

^  3.00 

!=  2.00 

is, 

1.00 

0.00 

-1.00 

-2.00 


|  Shl/NTC  Q  Sh2/NTC 
*** 


£■ 


*** 

■*** 


Roy  ef  al.  Cancer  &  Metabolism  201 4,  2:1 3 
http://www.cancerandmetabolism.eom/content/2/1/13 


Page  10  of  17 


r  \ 
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Figure  3  Loss  of  HSulf-1 -mediated  increase  in  fatty  acids,  sphingolipids,  and  lysolipids.  The  samples  were  extracted  using  Metabolon's 
standard  solvent  extraction  method  from  cells  in  logarithmic  phase  having  five  biological  replicates  for  each  sample  and  distributed  into  equal  parts 
for  analysis  on  the  GC7MS  and  LC/MS/MS  platforms.  (A)  Fold  increase  of  long-chain  fatty  acids  were  calculated  by  the  average  metabolite  level  of 
Shl/NTC  and  Sh2/NTC.  *p  =  0.02  to  0.09;  **p  =  0.002  to  0.01;  ***p  <  0.001  compared  to  NTC.  (B)  Branched  chain  fatty  acids  and  (C)  essential  fatty  acids 
were  calculated  by  the  average  metabolite  level  of  Shl/NTC  and  Sh2/NTC.  p  values  for  both  branched  chain  fatty  acids  and  essential  fatty  acids  were 
*p  =  0.01  to  0.02;  **p  =  0.002  to  0.01 ;  ***p  <  0.001  compared  to  NTC.  (D)  Fold  increase  of  sphingolipids  were  calculated  by  the  average  metabolite  level 
of  Shi  /NTC  and  Sh2/NTC.  *p  =  0.01  to  0.09;  **p  =  0.002  to  0.008;  ***p  <  0.001  compared  to  NTC.  (E)  Choline/inositol  pathway  metabolites'  fold  change 
and  (F)  membrane  lysolipids  were  calculated  as  mentioned  earlier.  P  values  for  choline/inositol  pathway  was  *p  =  0.009  to  0.09;  **p  =  0.002  to  0.008; 

***p  <  0.001  whereas  *p  =  0.01  to  0.02;  **p  =  0.002  to  0.009;  ***p  <  0.001  for  membrane  lysolipids  compared  to  NTC. 


Interestingly,  FASN,  SPHK1,  PLA2G4A,  PLA2G3,  SRE 
BF1,  and  PPARy  are  the  key  enzymes  which  are  localized 
and  involved  in  the  biogenesis  of  lipid  droplets  (LD) 
[50-54] .  The  FAs  are  activated  by  covalent  modification  by 
CoA  via  fatty-acyl-CoA  synthetases  and  esterified  to  gly¬ 
cerol  generating  triglycerides  or  sterol  esters  that  are  stored 
in  lipid  droplets  [42],  Accumulation  of  cytoplasmic  LDs 
forms  a  basis  of  increased  growth  and  chemoresistance  in 
neoplastic  cells  [40,55].  To  explore  whether  increased  ex¬ 
pression  of  above  enzymes  resulted  in  LD  formation,  we 
imaged  the  LDs  by  bodipy  staining  which  clearly  showed 
increased  LDs  in  Shi  and  Sh2  cells.  As  expected,  there  were 
very  few  LDs  in  the  NTC  cells  (Figure  4E).  To  confirm 


whether  HSulf-1  was  directiy  involved  in  LD  biogenesis,  we 
also  rescued  the  expression  of  HSulf-1  in  Shi  cells  by  stable 
transfection  of  CMV-driven  HSulf-1  expression  construct 
(clone  11)  with  vector-transfected  Shi  cells  served  as  con¬ 
trols  (Figure  4F,G).  Bodipy  staining  revealed  that  rescue  of 
HSulf-1  significantly  reduced  the  number  of  lipid  droplets 
in  Shi  (Shi  clone  11)  compared  to  Shi  vector  control  with 
many  LDs  (Figure  4G)  demonstrating  that  loss  of  HSulf-1 
promotes  LD  biogenesis.  Additionally,  transient  transfec¬ 
tion  of  PLA2G3  into  OV202NTC  cells  showed  increased 
LDs  in  these  cells  compared  to  vector-transfected  controls 
(Figure  4H),  suggesting  a  role  of  PLA2G3  in  loss  of  HSulf- 
1  mediated  LD  biogenesis. 
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Figure  4  Enhanced  expression  of  lipogenic  genes  in  ovarian  cancer.  (A)  Heat  map  of  a  subset  of  significantly  altered  lipid  pathway-related 
genes  by  supervised  clustering  {Red,  overexpressed  and  green,  downregulated  genes  in  Shi  and  Sh2  compared  to  NTC  cells).  (B)  Normalized  levels  of 
SREBPlc  { SREBF1 ),  SPHK1,  PLA2G3,  PLA2G4A  ( C-PLA2 )  and  PPARy  mRNA  by  quantitative  RT-PCR  in  NTC  and  Shi  and  Sh2  cells.  *p  <  0.05;  **p  <  0.01; 

***p  <  0.001  compared  to  NTC  cells.  (C)  Western  blot  analysis  of  FASN,  SREBF1,  PPARy,  PLA2G3  with  beta  tubulin  as  loading  control.  (D)  FASN  activity  is 
expressed  as  fold  change  in  Shi  and  Sh2  cells  compared  to  NTC  cells.  Results  are  means  ( columns )  of  two  independent  experiments  made  in  triplicate. 
One-factor  ANOVA  was  used  to  analyze  the  differences  in  FASN  activity  between  each  experimental  condition.  All  statistical  tests  were  two-sided 
(**p<  0.01).  (E)  Bodipy  {green)  and  DAPI  {blue)  staining  of  the  lipid  droplets  in  NTC,  Shi,  and  Sh2  cells  imaged  in  Carl  Zeiss  LSM  51  OS  confocal  microscope. 
(F)  Immunoblot  analysis  of  HSulf-1  in  Shi  vec,  and  Shi  Cl  1 1,  where  (S-actin  was  used  as  loading  control.  (G)  Bodipy  {green)  and  DAPI  {blue)  staining  of 
lipid  droplets  and  nuclei  respectively  in  Shi  vec  and  Shi  Cl  1 1  cells.  (H)  Bodipy  staining  of  LDs  following  enhanced  expression  of  empty  vector  compared 
to  PLA2G3  construct  by  transient  transfection  in  OV202NTC  cells  shows  LDs  only  in  PLA2G3  transfected  cells. 
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Loss  of  HSulf-1  facilitates  enhanced  (J-oxidation  and 
lipolysis 

The  analysis  of  metabolite  data  additionally  revealed  the 
increased  levels  of  acetylcarnitine  (FC  =  3.74,  p  =  <0.001 
and  FC  =  3.22,  p  =  <0.001),  butyrylcarnitine,  hexanoylcar- 
nitine  and  octanoylcarnitine  in  HSulf-1  silenced  cells. 
Moreover,  increased  levels  of  oleolycarnitine  were  16.84- 
and  6.64-fold  in  Shi  and  Sh2,  respectively  (Figure  5A). 
Carnitine  is  important  for  shuttling  FAs  across  mitochon¬ 
drial  membranes  for  oxidation  and  so  the  escalation  of 
carnitine  level  enhances  [i-oxidation  support  the  increased 
growth  of  ovarian  cancer.  To  determine  if  (1-oxidation  was 
enhanced  in  Shi  and  Sh2  cells,  we  measured  the  mRNA 
and  protein  level  of  carnitine  palmitoyltransferase  1 
( CPT1A )  and  the  results  showed  a  higher  level  of  CPT1A 
in  both  in  Shi  and  Sh2  cells  (Figure  5B).  The  increase  in 
the  expression  of  CPT1A  mRNA  along  with  increased 
production  of  several  of  the  long  chain  FAs  in  Shi  and 
Sh2  cells  indicated  that  these  cells  may  utilize  long  chain 
FAs  by  (1-oxidation  to  generate  more  ATP  to  accommo¬ 
date  increase  proliferation.  Next,  to  confirm  this  notion, 
we  determined  the  fatty  acid  oxidation  (FAO),  and  results 
showed  that  FAO  was  higher  both  in  Shi  and  Sh2  com¬ 
pared  to  NTC  cells  upon  addition  of  palmitate  (Figure  5C). 
Subsequently,  we  demonstrated  that  Etomoxir,  a  spe¬ 
cific  inhibitor  of  CPT1A,  decreased  FAO  more  sig¬ 
nificantly  in  Shi  and  Sh2  cells  compared  to  NTC 
cells.  We  also  demonstrated  that  Shi  and  Sh2  cells 
were  more  sensitive  toward  Etomoxir  treatment  than 
NTC  confirming  a  major  role  of  FAO  in  their  survival 
(Additional  file  8:  Figure  S4). 

Since  FAs  are  stored  in  lipid  droplets  and  released  by 
the  action  of  lipolytic  enzymes,  we  determined  the  expres¬ 
sion  of  monoacylglycerol  lipase  (MAGL),  diacylglycerol 
lipase  alpha  (DAGLA),  long-chain  acyl-CoA  synthetase 
(ACSL1),  and  hormone-sensitive  lipase  (HSL).  Immuno- 
blot  analysis  showed  that  the  expression  of  these  enzymes 
was  upregulated  in  Shi  and  Sh2  cells  compared  to  NTC 
cells  (Figure  5D),  demonstrating  that  loss  of  HSulf-1  acti¬ 
vates  beta  oxidation  and  lipolysis. 

Altered  metabolism  in  amino  acid  and  peptide  super 
pathways  in  cells  with  loss  of  HSulf-1 

In  addition  to  the  altered  lipid  metabolism,  the  metabolo- 
mic  analysis  revealed  that  amino  acids  and  their  deriva¬ 
tives  were  the  second  most  class  of  metabolites  altered  in 
HSulf-1  knockdown  cells  (Additional  file  4:  Table  S2). 
Most  of  the  amino  acids  including  serine,  threonine,  as¬ 
partate,  asparagine,  alanine,  phenylalanine,  tyrosine,  tryp¬ 
tophan,  arginine,  ornithine,  proline,  and  methionine  were 
decreased  in  the  range  of  -1.36-  to  -3.33-fold  (p  <  0.05) 
in  both  OV202  Shi  and  Sh2  cells.  In  contrast,  glycine 
was  increased  marginally  in  Shi  and  Sh2  cells  (range 
1.08-  to  1.1-fold;  p  <  0.001  for  each).  Interestingly,  we 


found  significant  upregulation  of  cysteine  (eight  to 
tenfold,  p  <  0.001)  and  its  derivatives,  hypotaurine 
(FC  =  2.38,  p  =  <0.001  and  FC  =  1.99,  p  =  < 0.0244 )  and  tau¬ 
rine  (FC  =  9.91,  p  =  <0.0032  and  FC  =  5.38,  p  =  < 0.0090 )  in 
both  Shi  and  Sh2  cells  (Additional  file  4:  Table  S2).  We 
also  observed  alterations  of  AAacetylated  amino  acids 
including  AT-acetyl-alanine,  serine,  and  threonine  resulted 
from  the  action  of  AT-acetyltransferases  on  acetyl-CoA 
and  L-amino  acids  where  all  the  six  AAacetylated  amino 
acid  tested  were  downregulated.  The  amino  acid- 
derived  antioxidant,  both  reduced  and  oxidized  gluta¬ 
thione  levels  were  lower  in  HSulf-l-deficient  cells. 
Additionally,  among  different  dipeptide  molecules,  gamma- 
glutamylglutamate  (FC  =  3.73,  2.77  in  Shi  and  Sh2, 
p  =  <.0001)  and  gamma-glutamylmethionine  (FC  =  7.31, 
6.23  in  Shi  and  Sh2  respectively,  p  =  <.0001)  levels  were 
increased  upon  loss  of  HSulf-1.  This  enhancement  may 
be  particularly  relevant  as  increased  expression  of  gamma- 
glutamylglutamate  is  reported  to  be  associated  with  tumor 
progression  and  drug  resistance  observed  in  human 
malignancies  [56].  Additionally,  a  significant  increase  of 
pro-hydroxy-proline,  (FC  =  2.05,  p  =  0.0057  and  FC  =  1.92, 
p  =  0.0088)  a  dipeptide,  also  demonstrated  in  Shi  and 
Sh2  cells.  Pro-hydroxy-proline  is  a  marker  of  collagen 
and  extracellular  matrix  degradation  (Additional  file  4: 
Table  S2). 

Discussion 

Aberrant  cellular  metabolism  in  cancer  is  now  well 
known  and  is  directly  related  to  tumorigenesis  in  most 
of  the  cancers  [57,58].  Multiple  signaling  pathways  and 
several  molecules  are  involved  in  the  synthesis  and  deg¬ 
radation  of  the  lipids  and  also  the  activities  of  lipid  me¬ 
tabolizing  enzymes  are  regulated  by  a  complex  interplay 
between  metabolic,  tumor  suppressor,  and  oncogenic 
signaling  [59].  Though  loss  of  HSulf-1,  a  putative  tumor 
suppressor  gene,  was  well  known  to  promote  tumorigen¬ 
esis,  angiogenesis  [22,23],  and  invasion  [20]  in  breast 
[18,21]  and  ovarian  [23,24]  cancers,  the  role  of  HSulf-1 
was  never  elucidated  in  altered  metabolism  of  ovarian 
cancer  cells  so  far.  Thus,  in  the  present  study,  we  have 
shown  that  loss  of  the  putative  tumor  suppressor,  HSulf- 
1  promotes  altered  metabolic  pathways  including  lipid, 
amino  acid,  and  nucleotide.  Of  the  several  pathways  al¬ 
tered  by  loss  of  HSulf-1,  the  lipid  and  amino  acid 
pathway-related  metabolites  accounted  for  50%  of  the 
total  altered  metabolites  identified.  In  ovarian  cancer,  al¬ 
tered  lipid  metabolism  was  detected  in  patients  during 
early  and  late  stages  of  disease  compared  to  healthy  con¬ 
trols  [60,61].  In  contrast  to  lipid  pathway-related  metabo¬ 
lites,  there  are  very  few  reports  on  the  alteration  in  the 
amino  acid  levels  in  OVCa.  We  found  a  significant  down- 
regulation  of  most  of  the  amino  acids  with  loss  of  HSulf-1 
consistent  with  Zhang  et  al.'s  report  where  they  showed 
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Figure  5  Loss  of  HSulf-1  induced  enhanced  p-oxidation  and  lipolytic  enzymes.  (A)  The  samples  were  extracted  using  Metabolon's  standard 
solvent  extraction  method  from  cells  in  logarithmic  phase  having  five  biological  replicates  for  each  sample  and  distributed  into  equal  parts  for  analysis 
on  the  GC/MS  and  LC/MS/MS  platforms.  Fold  increase  of  carnitine  and  its  derivatives  were  calculated  by  the  average  metabolite  level  of  Shl/NTC  and 
Sh2/NTC.  *p  =  0.01  to  0.09;  **p  =  0.001  to  0.009;  ***p  <  0.001  compared  to  NTC.  (B).  Levels  of  CPT1A  in  NTC,  Shi ,  and  Sh2  cells  determined  by  real-time 
PCR  and  Western  blot  analysis.  *p  <  0.05;  **p  <  0.01 ;  ***p  <  0.001  compared  to  NTC.  (C)  FAO  in  terms  of  OCR  in  pmol/min/mg  of  protein  was  monitored 
using  a  Seahorse  Bioscience  Extracellular  Flux  Analyzer  in  real  time  (mean  ±  S.D.,  n  =  3).  Cells  treated  with  etomoxir  (50  pM),  an  inhibitor  of  carnitine 
palmitoyltransferase  1,  served  as  a  positive  control.  Changes  in  the  FAO  induction  in  Shi  and  Sh2  cells  are  compared  with  that  of  NTC  cells. 
Etomoxir-induced  inhibition  of  FAO  in  Shi  and  Sh2  cells  are  compared  with  the  FAO  inhibition  in  NTC  (p  <  0.01  and  p  <  0.001).  (D)  Immunoblot  analysis 
of  MAGL,  DAGLA,  HSL,  ASCL1  in  NTC,  Shi,  and  Sh2  where  tubulin  used  as  internal  loading  control. 
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presence  of  lower  levels  of  amino  acids  in  the  serum  of  pa¬ 
tients  with  esophageal  adenocarcinomas  [62].  This  down- 
regulation  could  be  due  to  differences  in  uptake  from  the 
media,  increased  utilization,  and/or  catabolism  since  in¬ 
creased  demand  for  utilization  of  amino  acids  has  been  re¬ 
ported  in  other  cancers  [63,64].  Other  studies  have  shown 
an  increase  in  amino  acid  levels  in  the  serum  of  patients 
with  colon  and  breast  cancers  [65,66].  While  the  majority 
of  amino  acids  were  downregulated,  we  found  significant 
upregulation  of  cysteine  and  its  derivatives  hypotaurine 
and  taurine  in  Shi  and  Sh2  cells.  Significant  upregulation 
in  cysteine,  taurine,  and  hypotaurine  could  be  the  result  of 
increased  levels  of  cysteine  dioxygenase  type  I  (CDOl) 
[67-69]  in  Shi  and  Sh2  cells  (data  not  shown).  While  tau¬ 
rine  is  reported  to  be  downregulated  in  cancer  [70],  it  is 
also  reported  to  be  increased  in  tumors  of  the  prostate, 
squamous  cell  carcinoma  and  liver  metastasis  [71,72]. 
Additionally,  levels  of  myo-inositol  and  taurine  concentra¬ 
tions  both  in  vivo  and  in  vitro  are  correlated  with  cell 
density  of  the  tumors  [73]. 

Additionally,  we  also  saw  an  increase  in  the  dipeptide 
pro-hydroxy-pro  in  Shi  and  Sh2  cells.  Increase  in  pro¬ 
hydroxy-proline  is  consistent  with  the  report  that  major 
cartilage  ECM  proteins  type  II  collagen  and  aggrecan 
were  significantly  lower  in  HSulJ ^  ~  chondrocytes  sug¬ 
gesting  that  loss  of  HSulf-1  may  regulate  the  overall  bal¬ 
ance  of  cartilage  matrix  synthesis  and  degradation  [74]. 
Consistent  with  this,  significant  alterations  in  several 
UDP-glycosylation  moieties,  including  UDP-acetylglu- 
cosamine/UDP-acetylgalactosamine  UDP-glucose,  UDP- 
glucuronate,  and  UDP-galactose  are  also  supportive  of 
changes  in  extracellular  matrix  remodeling  with  reduced 
HSulf-1  expression. 

The  major  pathway  identified  by  both  ingenuity  pathway 
analysis  and  metabolic  profiling  was  the  alterations  in  the 
lipid  pathway.  Lipids  are  as  important  building  blocks  as 
carbohydrates  to  form  the  basic  skeleton  of  rapidly  dividing 
cells,  and  therefore,  large  amounts  of  FAs  are  required  to 
accommodate  high  rates  of  proliferation  in  cancer  cells  [4], 
Additionally,  the  source  of  FAs  may  determine  the 
phospholipid  composition  of  membranes.  Moreover,  it  was 
reported  that  in  high-grade  ovarian  cancer,  long-chain  fatty 
acids  were  elevated  [75].  Our  results  indicate  that  with 
the  loss  of  HSulf-1,  there  was  a  significant  increase  in 
the  long-chain  FAs  along  with  branched  and  essential 
FAs  in  ovarian  cancer.  These  high  levels  of  FAs  in  ovar¬ 
ian  cancer  with  the  absence  of  HSulf-1  indicate  a  major 
role  of  HSulf-1  in  FA  synthesis.  This  enhanced  produc¬ 
tion  of  FAs  might  meet  the  structural  needs  of  highly 
proliferative  ovarian  cancer  cells'  requirement,  when  the 
tumor  suppressor  HSulf-1  is  lost. 

Sphingolipids  and  lysolipid  or  lysophospholipids  are 
bioactive  lipid  molecules  which  play  pivotal  role  in  can¬ 
cer  pathogenesis.  The  most  important  sphingolipids  are 


ceramide  and  sphingosine-1 -phosphate  and  the  balance 
between  these  two  define  the  cell's  fate.  In  cancer, 
sphingolipid  metabolism  is  altered  and  includes  changes 
in  the  levels  of  sphingolipids  and  the  enzymes  involved 
in  their  metabolism  [76].  Lysolipids  with  their  various 
FA  side  chains  are  also  bioactive  lipids  which  mainly 
function  as  growth-stimulating  factor  and  induce  cell 
proliferation,  differentiation,  and  cell  migration  [77]. 
Herein,  we  report  an  increase  in  sphingolipids  including 
sphinganine,  sphingosine,  palmitoylsphingomyelin  along 
with  altered  choline/inositol  and  lysolipid  content.  Interest¬ 
ingly,  according  to  the  I<EGG  database,  metabolite  mapping 
and  glycerophospholipid,  arachidonic,  and  sphingolipid 
metabolisms  were  identified  as  key  junctions.  These  re¬ 
sults  clearly  indicate  that  with  the  loss  of  HSulf-1,  the  cell 
remodeled  its  own  lipid  synthesis  in  a  manner  to  supply 
both  structural  and  signaling  lipids  to  malignant  cells. 

In  cancer  cells,  the  increased  lipid  synthesis  is  due  to  the 
higher  expression  and  activity  of  lipogenic  enzymes 
[78,79].  Changes  in  the  expression  and  activity  of  enzymes 
involved  in  lipid  metabolism  are  regulated  by  metabolic 
and  oncogenic  signalling  pathways  [80,81].  We  identified 
enhanced  expression  of  FASN,  SREBF1,  SPHK1,  and 
PPARy  related  to  lipid  biogenesis.  Upregulation  of  these 
key  lipogenic  enzymes  with  the  loss  of  HSulf-1  were 
further  corroborated  by  their  mRNA  and  protein  ex¬ 
pressions.  Moreover  higher  enzymatic  activity  of  FASN 
strongly  supports  the  enhanced  biogenesis  of  FAs  with 
HSulf-1  loss.  Additionally,  phospholipase-related  genes 
PLA2G3,  PLA2G4A,  PNPLA2,  PLA2G4C,  and  PNPLA6 
involved  in  lipolysis  were  also  found  to  be  higher  in  Shi 
and  Sh2.  Interestingly,  PLA2G3  is  an  important  enzyme  in 
lipid  metabolism  and  implicated  in  LD  biogenesis  [50,52]. 
Recently,  it  was  reported  that  cancer  cells  contain  in¬ 
creased  numbers  of  lipid  droplets  compared  with  normal 
tissue  [82]  which  are  storage  sites  for  triglycerides  and 
cholesterol  to  be  used  as  energy  source.  Loss  of  HSulf-1 
induced  an  increase  in  PLA2G3  expression  that  may  lead 
to  the  accumulation  of  LD  in  ovarian  cancer.  Furthermore, 
reduction  in  the  number  of  LDs  with  the  rescue  of  HSulf- 
1  expression  in  Shi  cells  establishes  a  direct  connection  to 
the  loss  of  a  putative  tumor  suppressor  HSulf-1  and  en¬ 
hanced  LD  biogenesis  in  ovarian  cancer. 

In  accordance  with  these  results,  we  also  demonstrated 
lower  level  of  carnitine  along  with  higher  level  of  carni¬ 
tine  derivatives  suggestive  of  enhanced  transport  of  FAs 
through  the  mitochondrial  membrane.  More  import¬ 
antly,  the  augmented  CPT1A  level  in  Shi  and  Sh2  cells 
strongly  support  the  possibility  of  enhanced  (3-oxidation. 
Increased  and  altered  FA  synthesis  in  cancer  is  a  well- 
accepted  phenomenon;  however,  enhanced  expression  and 
activity  of  lipolytic  enzymes  in  tumor  cells  is  a  recent  ob¬ 
servation.  In  prostate  cancer,  malignant  cells  solely  rely  on 
FAO  as  their  energy  resource  [83].  Here,  we  report  that 
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loss  of  HSulf-1  also  promotes  enhanced  FAO  in  Shi  and 
Sh2  cells  compared  to  NTC.  Consistent  with  increased 
lipolysis,  our  results  indicated  a  significant  increase  in  the 
expression  of  important  lipolytic  enzymes  including 
MAGL,  DAGLA,  HSL,  and  ASCL1  in  Shi  and  Sh2  cells 
compared  to  NTC  cells.  Higher  MAGL  expression  has 
been  shown  to  stimulate  pro-tumorigenic  signals  and  pro¬ 
motes  survival,  tumor  growth,  and  migration  [84].  In  this 
regard,  the  relationship  between  synthesis,  storage,  and 
utilization  of  free  fatty  acids  (FFA)  through  beta  oxidation 
for  energy  production  is  very  poorly  understood.  While 
the  prevailing  belief  that  cells  that  have  increased  FA  syn¬ 
thesis  usually  do  not  undergo  beta  oxidation,  it  has  also 
been  suggested  that  newly  synthesized  FFA  are  immedi¬ 
ately  converted  into  neutral  lipids  and  stored  in  lipid  drop¬ 
lets  [84].  The  FFAs  are  released  by  the  action  of  lipases 
such  as  mono-acyl  glycerol  (MAGL)  which  is  then  used 
for  new  membrane  synthesis,  lipid  signaling  and  beta  oxi¬ 
dation  for  energy  production  for  the  anabolic  reactions. 
Our  data  seems  to  support  this  hypothesis  since  the  Shi/ 
Sh2  cells  have  both  increased  the  rate  of  FA  synthesis  and 
beta  oxidation  compared  to  NTC  cells.  The  FFA  may  be 
immediately  converted  to  neutral  lipids  and  stored  in  lipid 
droplets  to  be  released  under  stress  conditions. 

Increasing  evidence  indicates  that  metabolic  alterations 
induced  by  loss  of  tumor  suppressor  genes  are  common 
in  cancer  [16].  These  alterations  are  critical  for  growth 
and  survival  of  cancer  cells.  This  is  the  first  comprehen¬ 
sive  report  on  metabolic  alterations  induced  by  loss  of 
HSulf-1,  a  putative  tumor  suppressor  gene  in  ovarian  and 
breast  cancer.  Our  results  clearly  indicate  that  loss  of 
HSulf-1  remodels  lipid  metabolism  in  ovarian  cancer.  It  is 
well  accepted  that  HSulf-1  is  a  major  regulator  of  growth 
factor-mediated  signaling  and  altered  the  tumor  micro¬ 
environment  in  ovarian  cancer;  however,  how  HSulf-1  dir¬ 
ectly  reconstructs  the  lipid  metabolism  in  terms  of  FA 
synthesis,  lipolysis,  and  enhanced  LD  formation  is  yet  to 
be  unraveled.  A  better  understanding  of  the  function  of 
HSulf-1,  its  role  in  energy  metabolism  (glycolysis,  TCA  cy¬ 
cles,  PPP)  and  its  regulation  on  target  genes  involved  in  al¬ 
tered  cellular  metabolism  could  lead  to  a  better 
understanding  of  tumorigenesis  as  well  as  development  of 
new  targeted  therapy. 

Conclusions 

HSulf-1  is  reported  to  be  an  important  tumor-suppressor 
gene  and  its  expression  is  lost  in  a  majority  of  ovarian  tu¬ 
mors.  This  study  demonstrates  a  significant  alteration  of 
cellular  metabolism  upon  loss  of  HSulf-1  in  OV202  cells. 
Microarray  analysis  and  metabolite  profiling  is  performed 
to  ascertain  the  impact  of  HSulf-1  on  the  overall  meta¬ 
bolic  changes  including  alterations  in  lipid  and  amino  acid 
pathways.  The  major  finding  from  our  study  shows  the 
metabolic  reprogramming  of  cells  toward  an  enhanced 


lipid  metabolism  upon  absence  of  HSulf-1.  Cells  adopt  a 
lipogenic  phenotype  which  is  manifested  with  an  excess 
fatty  acid  synthesis  and  an  upregulated  beta-oxidation. 
Furthermore,  HSulf-l-deficient  cells  accumulate  a  huge 
amount  of  cytoplasmic  lipid  droplets  to  accommodate  the 
excessive  fatty  acid  syntheses.  The  present  findings  are 
supported  by  the  increasing  evidences  of  enhanced  lipid 
metabolism  in  cancer  cells  facilitating  cell  survival,  prolif¬ 
eration,  and  signaling.  Our  results  indicate  that,  loss  of 
HSulf-1  is  enabling  the  cells  to  synthesize  more  lipids  to 
expedite  the  high  proliferation  rate  and  survival. 

Additional  files 


Additional  file  1:  Figure  SI.  Loss  of  HSulf-1  induced  enhanced 
proliferation  in  OV202  cells.  Equal  number  of  cells  of  NTC,  Shi,  and  Sh2 
were  plated  in  triplicates  (1  x  1 0s)  and  counted  after  24, 48,  and  72  h. 

The  increase  in  cell  count  in  both  Shi  and  Sh2  were  statistically  significant 
in  48  h  (*p  <  0.05)  and  72  hr  (**p  <  0.001).  These  experiments  were  repeated 
twice. 

Additional  file  2:  Table  SI.  Microarray  analyses  of  OV202  cells  upon 
HSulf-1  loss. 

Additional  file  3:  Figure  S2.  Multivariate  regression  results  between 
major  metabolite  classes.  The  eight  major  metabolite  classes  are  amino  acid, 
carbohydrate,  cofactors  and  vitamins,  energy,  lipid,  nucleotide,  peptide,  and 
xenobiotics.  'Estimate'  and  'Std.  Error'  are  estimate  of  regression  coefficient 
and  estimation  of  standard-deviation  error;  't  value'  is  the  f  statistics  for  each 
coefficient  estimate,  and  'Pr(>|t|)'  is  the  corresponding  p  value  for  each 
coefficient  estimate. 

Additional  file  4:  Table  S2.  Changes  in  major  metabolic  pathways 
including  amino  acids,  lipids  and  nucleotides.  List  of  metabolites  identified 
through  mass  spectrometry  and  the  super-pathway  and  sub-pathway  for 
them  are  shown.  The  green-  and  red-shaded  boxes  are  metabolites  that  are 
downregulated  and/or  upregulated  in  Shi  and  Sh2  cells  compared  to  NTC 
cells  respectively.  The  p  values  and  the  platform  used  to  identify  these 
metabolites  are  also  shown. 

Additional  file  5:  Table  S3.  Global  metabolic  changes  upon  HSulf-1  loss. 
List  of  significantly  altered  metabolites  in  Shi  and  Sh2  cells  compared  to 
NTC  cells  in  different  pathways  (XLS).  Heat  map  of  statistically  significant 
biochemical  profiles  in  this  study.  By  paired  comparisons,  shaded  cells 
indicate  p  <  0.05  (red  indicates  that  the  mean  values  are  significantly  higher 
compared  untreated  control;  green  values  significantly  lower).  Blue-bolded 
text  indicates  0.05  <  p  <  0.1 0.  All  data  were  normalized  using  Bradford 
protein  concentration  (red — upregulated,  green — downregulated,  and 
blue — approaching  significance). 

Additional  file  6:  Table  S4.  Alterations  in  lipid  metabolites.  Lipid 
metabolic  pathways  in  Shi  and  Sh2  cells  compared  to  NTC  cells.  Heat 
map  of  statistically  significant  (p  <  0.05,  FDR  <  0.05)  biochemicals  are 
shown.  Red  indicates  that  the  mean  values  are  significantly  higher  by 
comparison  with  untreated  control;  green  values  significantly  lower. 

All  data  were  normalized  using  Bradford  protein  concentration 
(red — upregulated,  green — downregulated. 

Additional  file  7:  Figure  S3.  Metaboanalyst  pathway  analysis.  (A) 
Statistics  for  pathways  with  major  change  based  on  high  impact  (linoleic 
acid  metabolism)  or  p  value  (pathways  glycerophospholipid,  arachidonic, 
and  sphingolipid  metabolic  pathways).  (B)  Of  the  12  highly  significant  KEGG 
pathways  plotted  according  to  global  test  p  value  (intensity  of  color  in  the 
vertical  axis)  and  impact  factor  (size  of  the  circles  in  the  horizontal  axis),  all  4 
belong  to  the  lipid  pathways. 

Additional  file  8:  Figure  S4.  Effect  of  etomoxir  on  cellular  growth  with 
increasing  concentration  of  Etomoxir  treatment  (0  to  100  pM)  in  NTC, 

Shi,  and  Sh2  cells  (n  =  2).  At  60  pM  and  onwards,  the  cell  growth 
inhibition  was  statistically  significant  (p  <  0.05)  in  both  Shi  and  Sh2  cells 
compared  to  NTC. 
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ACSL1:  long-chain  acyl-CoA  synthetase;  CPT1A:  carnitine  palmitoyltransferase 
1;  DAGLA:  diacylglycerol  lipase  alpha;  FDR:  false  discovery  rate;  FA:  fatty  acid; 
FAO:  fatty  acid  oxidation;  FASN:  fatty  acid  synthase;  GC/MS:  LC/MS/MS: 

PCA:  principal  component  analysis;  HSL:  hormone  sensitive  lipase; 

KEGG;  Kyoto  Encyclopedia  of  Genes  and  Genomics;  LD:  lipid  droplet; 

MAGL:  monoacylglycerol  lipase;  NTC:  OV202:  non-targeted  control  cells; 
PPARy:  peroxisome  proliferator-activated  receptor  gamma;  SAM:  significance 
analysis  of  microarray;  Shi:  stably  knockdown  HSulf-1  clonal  line  1  in  OV202; 
Sh2:  stably  knockdown  HSulf-1  clonal  line  2  in  OV202;  SPHK1 :  sphingosine 
kinase  1;  SREBF1:  sterol  regulatory  element-binding  transcription  factor  1. 
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ABSTRACT 

Warburg  effect  has  emerged  as  a  potential  hallmark  of  many  cancers.  However, 
the  molecular  mechanisms  that  led  to  this  metabolic  state  of  aerobic  glycolysis, 
particularly  in  ovarian  cancer  (OVCA)  have  not  been  completely  elucidated.  HSulf-1 
predominantly  functions  by  limiting  the  bioavailability  of  heparan  binding  growth 
factors  and  hence  their  downstream  signaling.  Here  we  report  that  HSulf-1,  a  known 
putative  tumor  suppressor,  is  a  negative  regulator  of  glycolysis.  Silencing  of  HSulf-1 
expression  in  OV202  cell  line  increased  glucose  uptake  and  lactate  production  by 
upregulating  glycolytic  genes  such  as  Glutl,  HKII,  LDHA,  as  well  as  metabolites. 
Conversely,  HSulf-1  overexpression  in  TOV21G  cells  resulted  in  the  down  regulation 
of  glycolytic  enzymes  and  reduced  glycolytic  phenotype,  supporting  the  role  of 
HSulf-1  loss  in  enhanced  aerobic  glycolysis.  HSulf-1  deficiency  mediated  glycolytic 
enhancement  also  resulted  in  increased  inhibitory  phosphorylation  of  pyruvate 
dehydrogenase  (PDH)  thus  blocking  the  entry  of  glucose  flux  into  TCA  cycle. 
Consistent  with  this,  metabolomic  and  isotope  tracer  analysis  showed  reduced 
glucose  flux  into  TCA  cycle.  Moreover,  HSulf-1  loss  is  associated  with  lower  oxygen 
consumption  rate  (OCR)  and  impaired  mitochondrial  function.  Mechanistically,  lack 
of  HSulf-1  promotes  c-Myc  induction  through  HB-EGF-mediated  p-ERK  activation. 
Pharmacological  inhibition  of  c-Myc  reduced  HB-EGF  induced  glycolytic  enzymes 
implicating  a  major  role  of  c-Myc  in  loss  of  HSulf-1  mediated  altered  glycolytic 
pathway  in  OVCA.  Similarly,  PG545  treatment,  an  agent  that  binds  to  heparan  binding 
growth  factors  and  sequesters  growth  factors  away  from  their  ligand  also  blocked 
HB-EGF  signaling  and  reduced  glucose  uptake  in  vivo  in  HSulf-1  deficient  cells. 
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INTRODUCTION 

Cancer  cell  survival  and  growth  vastly  depends  on 
the  deregulated  overlapping  signaling  cascades  which 
are  initiated  at  the  extracellular  matrix  by  a  range  of 
growth  factors  [1].  Constitutive  activation  of  these  signals 
instructs  the  cells  for  continual  biomass  production  to 
accumulate  metabolic  intermediates  as  the  source  of 
building  blocks.  Mounting  evidence  indicates  that  this 
altered  cell  metabolism  can  be  actively  regulated  by 
oncogenes  and  tumor  suppressors.  Among  the  tumor 
suppressors,  the  roles  of  P53  and  PTEN  leading  to 
metabolic  reprogramming  of  cancer  are  well  documented 
[2,  3].  However,  it  is  unclear  whether  other  putative  tumor 
suppressors  also  play  a  critical  role  in  altering  the  Warburg 
effect.  Thus,  a  comprehensive  understanding  of  the 
molecular  alterations  by  other  putative  tumor  suppressors 
or  oncogenes  that  lead  to  altered  metabolism  may  not  only 
shed  new  light  on  tumorigenesis  but  also  provide  novel 
strategies  for  therapeutically  targeting  these  alterations. 

Majority  of  the  growth  factors  like  FGF-1,-2  [4], 
VEGF  [5],  HB-EGF  [6]  and  cytokines  including  IE-6 
[7]  and  IL-8  [8]  depend  on  heparan  sulfate  proteoglycans 
(HSPGs)  for  binding  and  signaling  through  their  cognate 
receptor.  Heparan  sulfate  (HS)  is  a  glycosaminoglycan 
consisting  of  repeating  unbranched  negatively  charged 
disaccharide  units  variably  sulfated  at  the  3-0,  6-0,  or 
N  sites  on  glucosamine  and  the  6-0  site  on  glucuronic/ 
iduronic  acid  [9].  Growth  factors  and  cytokines  form  the 
ternary  complexes  with  their  cognate  receptors  and  HS 
resulting  in  ligand-mediated  activation.  We  had  previously 
reported  that  HSulf-1  (also  known  as  Sulfatase  1,  Sulf- 
1,  KIAA1077  and  Extracellular  Sulfatase  Sulf-1  [6]  is 
downregulated  in  a  majority  of  ovarian  cancer  cell  lines 
[6]  (Fig. SI)  and  tumors  including  serous,  endometrioid 
and  clear  cell  tumors  of  the  ovary  [10].  Also,  we  have 
demonstrated  that  loss  of  HSulf-1  modulates  the  signaling 
of  HS  binding  growth  factors  such  as  FGF-2,  VEGF,  HGF, 
and  HB-EGF  in  ovarian  [11],  head  and  neck  squamous 
carcinoma  [11]  and  metastatic  breast  carcinomas  [12] 
respectively  and  plays  an  important  role  in  tumor 
progression,  metastasis  and  angiogenesis  [10,  13,  14]. 
Further,  we  showed  that  HIF-la,  a  major  transcription 
factor  that  also  promotes  altered  metabolic  signature, 
negatively  regulates  HSulf- 1  expression  in  breast  cancer 
[15].  Moreover,  HSulf-1  silencing  increases  the  ability 
to  form  anchorage-independent  colonies  in  vitro  and 
enhanced  tumorigenicity  in  vivo  [16].  Other  studies  also 
demonstrated  that  HSulf-1  blocks  cell  proliferation, 
migration  and  growth  in  vitro  and  in  vivo  in  hepatocellular 
carcinoma  [17,  18]  and  suppresses  the  malignant  growth 
in  lung  and  gastric  cancer  by  inhibiting  ERK,  AKT  and 
hedgehog  signaling  respectively  [19,  20].  Based  on 
these  findings,  we  hypothesized  that  HSulf-1,  due  to  its 
regulation  of  growth  factor  mediated  signaling,  might 
play  a  critical  role  in  altering  cellular  metabolism.  Indeed, 


our  recent  findings  demonstrate  that  loss  of  HSulf-1 
potentially  contributes  to  the  metabolic  alterations  in  the 
lipogenic  phenotype  of  ovarian  cancer  cells  [21].  Here, 
we  investigated  whether  HSulf-1  deficiency  would  also 
affect  other  metabolic  pathways  such  as  glycolysis  and 
TCA  cycle.  By  combining  bioenergetics  and  metabolic 
studies,  our  results  show  for  the  first  time  that  growth 
factor  signaling  modulated  by  HSulf-1  loss  increases 
glucose  uptake  and  lactate  production  and  alters  the 
energy  metabolism  and  subsequently  promoting  c-Myc 
activation.  In  this  study  we  utilized  PG545,  a  novel 
synthetic  agent  currently  in  Phase  1  clinical  trials  (Clinical 
Trials  gov.identifier,  NCT02042781)  and  essentially 
mimics  HSulf-1  function.  PG545  functions  as  HS 
mimetic  by  simultaneously  blocking  HS-mediated  growth 
factor  signaling  leading  to  inhibition  of  angiogenesis  and 
carcinogenesis  [22-24]  including  in  ovarian  cancer  [22]. 
However,  whether  PG545  also  modulates  the  glycolytic 
phenotype  has  never  been  explored  before.  We  show  for 
the  first  time  that  PG545  treatment  resulted  in  significant 
reduction  in  glycolytic  phenotype  induced  by  loss  of 
HSulf- 1  and  downregulated  c-Myc  and  expression  of  key 
glycolytic  enzymes  in  vitro.  More  importantly,  PG545 
decreased  in  vivo  glucose  uptake,  lactate  production  and 
markedly  inhibited  tumor  progression. 

RESULTS 


HSulf-1  reprograms  the  glycolytic  pathway 

To  understand  the  impact  of  HSulf-1  on  glycolytic 
metabolism  in  OVCA  cells,  we  analyzed  the  levels  of 
glycolytic  genes  of  OV202  non-targeted  control  (NTC) 
and  HSulf- 1-ShRNA  downregulated  clones,  Shi  and 
Sh2  cells  [16].  The  microarray  analysis  (accession  no- 
GSE67086)  revealed  aberrant  glycolytic  gene  expression 
in  Shi  and  Sh2  compared  to  NTC  cells  (Figure  1A). 
Q-PCR  validation  showed  that  glycolytic  genes  including 
hexokinase  II  (HKII),  6-phosphofructo-2-kinase/fructose- 
2, 6-biphosphatase  3  (PFKFB3),  6-phosphofructokinase, 
liver  type  (PFKL),  aldolase  C,  fructose-bisphosphate 
(ALDOC)  and  related  genes  including  glucose  transporter 
1  (Glutl),  and  moncarboxylase  transporter  4  (MCT4) 
were  upregulated  in  Shi  and  Sh2  cells  (Figure  IB). 
Similar  results  were  also  observed  at  the  protein  level  by 
immunoblot  including  PGAM  and  PKM2.  Importantly, 
most  of  the  protein  levels  were  rescued  after  re-expression 
of  HSulf-1  in  Shi  (C17)  cells  (Figure  1C). 

To  determine  if  ectopic  expression  of  HSulf- 1  will 
reduce  glycolytic  enzyme  levels,  we  generated  a  HSulf- 
1 -overexpressing  stable  clone  in  TOV21G  (Figure  ID). 
Results  showed  that  enhanced  HSulf-1  expression  in 
TOV21G  clonal  line  (Cl  1 )  resulted  in  the  downregulation 
of  all  of  these  enzymes  compared  to  vector  control  (Figure 
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ID).  To  further  substantiate  loss  of  HSulf-1  induced 
alterations  in  glycolysis,  we  checked  the  enzyme  levels 
in  immortalized  Wild-Type  (WT)  and  HSulf-1  Knock-Out 
(KO)  MEFs  and  observed  higher  Glut  1,  F1KI1  and  LDHA 
levels  in  MEF  KO  cells  (Figure  IE). 

To  assess  the  role  of  putative  tumor  suppressor 
HSulf-1  in  other  cancer,  we  transiently  overexpressed 
HSulf-1  in  prostate  cancer  (PC3),  lung  cancer  (H358)  and 
breast  cancer  (MDA-231)  cells.  Results  showed  a  marked 
reduction  in  Glutl,  HK1I  and  LDHA  (Figures  1F-1H)  in 
HSulf- 1  expressing  cells  compared  to  vector  transfected 
controls.  These  results  suggested  that  HSulf-1  might  play 
a  significant  role  in  promoting  glycolysis  in  cancer. 


Loss  of  HSulf-1  induced  enhanced  glycolytic 
phenotype  in  ovarian  cancer  in  vitro  and  in  vivo 

To  evaluate  if  upregulation  of  glycolytic  genes  and 
enzyme  levels  upon  HSulf-1  loss  resulted  in  increased 
glucose  uptake,  we  stained  the  live  cells  using  2-NBDG 
a  fluorescent  glucose  analogue  which  is  used  to  monitor 
glucose  uptake  in  live  cells.  Data  shows  that  glucose 
uptake  was  enhanced  in  Shi  and  Sh2  cells  by  2.8  and  1.8 
fold  respectively  as  compared  to  NTC  cells  (Figures  2A- 
2B).  Interestingly,  re-expression  of  HSulf-1  resulted  in 
significant  reduction  in  glucose  uptake  in  C17compared 
to  Shi  cells  (Figures  2A-2B).  Moreover,  in  vivo  imaging 
using18FDG-PET  showed  a  marked  increase  in  glucose 
uptake  in  2  week  old  Shi -xenograft  but  not  in  the  NTC- 
xenograft  model  (Figures  2C-2D,  S2).  Our  results  also 
demonstrated  a  significant  increase  in  lactate  secretion 
and  higher  levels  of  ATP  in  Shi  and  Sh2  compared  to 


B 


</)(/)  O 


HSulf-1 

Glutl 

HKII 

PFKFB3 

ALDOC 

PGAM 

PKM2 

Tubulin 


&  1-5 

ut 

c 

»  1 
o 

2  0.5 
o 

“■  o 


II 


NTC  Shi  Sh2  CI7 
a>  4  Glutl 


o  2 

2 

o 

0 


a 


U 


LL 


NTC  Shi  Sh2  CI7 
>  4  HKII 


j,1-5 

\  i 

j  0.5 
) 

-  0 


NTC  Shi  Sh2  CI7 


M 


NTC  Shi  Sh2  CI7 


.c  o  - 

o  £- 

2 

o 

0  ■ 

o  1.5 
a ) 

I  1 

1 0-5 
if  o 

>4 
T> 

} 

I  2 

l 
> 

'  0 


a 


n  n 


NTC  Shi  Sh2  CI7 
PGAM 


_  3.0 

5) 

c 

<Q 

1  2-0 

o 

u. 

* 

2  i.o 
(0 

E 

o 

z 

o.o4- 

5.0 

a> 

O) 

jS  4.0 

O 

■o 

O  3.0 

LL 

T3 

jj  20 

re 

i  1.0 

o 

z 

o.o 


NTC  Shi  Sh2  CI7 
PKM2 


nil 


SLC2A1 

(Glutl) 


L  I 

NTC  Shi  Sh2 

PFKFB3 


2.0 


HK2 


PFKL 


o.o 


NTC  Shi  Sh2 


ALDOC 


1 


I 


NTC  Shi  Sh2 


nj  5.0 
£ 
o 

2  4.0 

O 

S  3.0 
0) 

—  2.0 
(0 

E 

o  1.0 
z 

0.0 


NTC  Shi  Sh2 


O 

2 

o 

u-  1.0 


0.0 

3.0i 

i 

i) 

2.5 

| 

j  2.0 

J 

■  1.5 

I 

I  J 

[ 

j  0.5 
0.0 


H358 


NTC  Shi  Sh2 

SLC16A3 
(MCT4) 


i 


NTC  Shi  Sh2 


hi 


HSulf-1 

Glutl 

HKII 

LDHA 

GAPDH 


MDA231 


NTC  Shi  Sh2  CI7 


Figure  1:  Absence  of  HSulf-1  augmented  glycolytic  key  enzymes.  A.  Microarray  analysis  of  glycolytic  genes  in  NTC,  Shland 
Sh2  cells  in  triplicates.  B.  QRT-PCR  analysis  of  relative  mRNA  levels.  C.  Immunoblot  analysis  of  glycolytic  enzymes  in  NTC,  Shi,  Sh2 
and  C17  cells  (N=  2).  Fold  change  as  determined  by  densitometric  analysis  of  these  enzymes  in  Shi  and  2  cells  compared  to  NTC  cells  are 
shown  as  bar  graph.  Immunoblot  analysis  of  glycolytic  enzymes  in  TOV21GVec  and  Clone  1  cells  D.,  in  WT  and  HSulf-1  KO  MEF  cells 
E..  Protein  levels  of  glycolytic  enzymes  after  transient  over-expression  of  HSulf-1  in  PC3  F.,  H358  G.  and  MDA231  H.  cells. 
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Figure  2:  HSulf-1  loss  induces  enhanced  glycolytic  phenotype  in  ovarian  cancer.  A.  Live  cell  imaging  showing  glucose  uptake 
in  NTC,  Shi,  Sh2  and  C17  using  2-NBDG.  B.  Quantification  of  fluorescence  intensity  of  the  glucose  uptake  in  OV202NTC,  Shi,  Sh2 
and  C17  cells.  C.  In  vivo  glucose  uptake  using  1SFDG  in  OV202NTC  and  Shi  xenografts.  D.  In  vivo  glucose  uptake  was  calculated  using 
PMOD  Biomedical  Image  Quantification  software,  *p  <  0.001.  E.  Lactate  secretion  in  media  of  OV202NTC,  Shi,  Sh2  and  C17  cells,  F. 
Intracellular  ATP  levels  in  OV202NTC,  Shi,  Sh2  and  07  cells,  (*p  and  **p  <0.05  compared  to  NTC,  #p<  0.05  compared  to  Shlj.Glucose 
uptake  G., glucose  uptake  intensity  H.,  lactate  secretion  I., and  intracellular  ATP  levels  J.  in  TOV21G  Vec  and  Cl  cells  where  *p  <  0.05. 
Glucose  uptake  K.,  its  intensity  L.,  lactate  secretion  M.,  and  intracellular  ATP  levels  N.  in  WT  and  KO  MEFs,  where  *p  <  0.05.  O.  Fold 
increase/decrease  of  glycolytic  metabolite  were  evaluated  using  GC/MS  and  LC/MS/MS  platforms  and  fold  change  were  calculated  by  the 
average  metabolite  level  ofShl/NTC,  Sh2/NTC  and  Shl-C17/Shl.  *p  <0.05. 
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NTC  cells,  while  this  phenotype  is  reversed  in  Shl-C17 
(Figures  2E-2F).  Conversely,  forced  expression  of  HSulf-1 
in  TOV21G  (Cl  1 )  showed  a  significant  reduction  (35.4%) 
in  glucose  uptake,  lactate  production  (21.5%)  and  ATP 
levels  (50%)  (Figures  2G-2J).  Similarly,  HSulf-1  KO 
MEFs  showed  84%  increase  in  glucose  uptake  when 
compared  to  WT-MEFs  (Figures  2M-2N),  and  exhibited 
22%  increase  in  lactate  secretion  and  37%  increase  in  ATP 
content  (Figures  2K-2N),  further  substantiating  that  loss 
of  HSulf-1 -induced  aerobic  glycolysis.  Also,  transient 
overexpression  of  HSulf-1  in  PC3,  H358  and  MDA231 
cancer  cell  lines  reduced  lactate  secretion  (Figure  S3). 
Additionally,  our  metabolite  analysis  demonstrated  that 
HSulf-1  deficiency  was  closely  associated  with  increased 
levels  of  glycolytic  metabolites  including  glucose 
6-phosphate  [fold  change  (FC)  compared  to  NTC  (FC 
=10.68  p  =<  0.001)],  fructose  1,6  bisphosphate  (FC=2.89, 
p=0.03),  3-phosphoglycerate  (FC  =  6.02,  p  =  0.01)  and 


intracellular  lactate  (FC  =  1.5,  p  =  0.002)  (Figure  20). 
Altogether,  these  results  suggest  a  key  role  of  HSulf-1  in 
modulating  the  glycolytic  phenotype. 

HSulf-1  decouples  the  energy  flow  from  glycolysis 
to  TCA  cycle 

To  elucidate  the  role  of  HSulf-1  in  the  next  level 
of  glucose  oxidation,  we  examined  the  inhibitory 
phosphorylation  of  PDH  (S-293)  along  with  PDK1 
expression  as  PDH  and  PDK  regulates  the  entry  of  acetyl 
CoA  into  TCA  cycle  [23].  Our  results  showed  increased 
p-PDH  in  El  complex  in  Shi  and  Sh2  compared  to 
NTC  cells  (Figure  3A).  Similar  results  were  observed  in 
TOV21G  Vec/Cll  (Figure  3B)  and  in  MEFs  (Figure  3C). 
It  is  well  established  that  increased  levels  of  pyruvate 
that  accumulate  due  to  higher  aerobic  glycolysis  is  either 
converted  to  acetyl  CoA  or  lactate  [24].  Since  the  level  of 
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Figure  3:  Disruption  of  energy  flow  from  glycolysis  to  TCA  cycles.  Immunoblot  analysis  of  p-PDH,  PDH,  and  LDHA  in  OV202 
clonal  cells  A.,  TOV21G  vec  and  Cll  B.  and  HSulf-1  WT  and  KO  MEF  cells  C..  LDH  activity  in  OV202  clonal  lines  D„  TOV21G  vec 
and  Cll  E.  and  in  MEFs  F.  where*/><  0.05  and  **p<  0.05  compared  to  Shi  cells  and  /or  compared  to  vector  transfected  and  WT-  MEF 
cells  respectively.  G.  Fold  increase/  decrease  of  TCA  metabolites  evaluated  on  metabolon  platform  and  fold  change  were  calculated  by  the 
average  metabolite  level  of  Shl/NTC,  Sh2/NTC  and  Shl-C17/Shl,  *p  <  0.05  and  #p  <  0.9.  H.  Isotope  tracer  analysis  using  13C6  labeled 
glucose  to  estimate  glucose’s  contribution  towards  TCA  cycle  metabolites  level.  Data  is  normalized  by  cell  number  for  OV202  NTC,  Shi 
and  Sh2  cells  ( n  =  7). 
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p-PDH  were  increased  with  loss  of  HSulf-1,  we  further 
monitored  the  activities  and  levels  of  LDH,  a  key  enzyme 
involved  in  conversion  of  pyruvate  to  lactate  in  order  to 
maintain  high  glycolytic  rate.  HSulf-1  downregulated 
Shi  and  Sh2  clones  showed  a  marked  increase  in 
LDHA  expression  (Figure  3A)  and  activity  (Figure  3D). 
Moreover,  we  found  that  with  the  HSulf- 1  re-expression, 
LDHA  expression  and  activity  were  reduced  in  C17  cells 
(Figure  3A  and  3D).  Likewise,  we  found  reduced  LDHA 
level  (Figure  3B)  and  activity  (Figure  3E)  in  the  HSulf-1 
over  expressing  TOV21GC1  compared  to  vector.  Similar 
results  were  also  observed  in  HSulf-1  KO  MEF  compared 
to  MEF  WT  cells  (Figure  3C  and  3F). 

We  next  investigated  whether  HSulf-1  deficiency 
reduced  glucose’s  contribution  into  TCA  cycle.  To  this 
end,  we  performed  unbiased  metabolic  profiling  of  TCA 
metabolites  using  the  Metabolon  platform  (Metabolon 
Inc,  NC)  in  NTC,  Shi,  Sh2  and  C17  cells.  Results  showed 
a  significant  reduction  of  most  of  the  TCA  metabolites 
including  citrate,  fumarate,  malate  and  succinate  in  cells 
with  loss  of  HSulf-1  (Figure  3G).  To  further  demonstrate 
that  loss  of  HSulf-1  reduces  TCA  cycle  fluxes  in  ovarian 
cancer  cells,  we  performed  13-C  GC-MS  based  isotope 
tracer  analysis  using  labeled  U-13C6  glucose  to  reveal  the 
distribution  of  metabolites  in  this  metabolic  pathway. 
As  seen  in  Figure  3H,  M6  (6  carbon  labeled)  glucose 
is  converted  into  M3  pyruvate  and  then  M2  acetyl- 
CoA  through  pyruvate  dehydrogenase.  M2  acetyl-CoA, 
combining  with  unlabeled  oxaloacetate  (MO  OAA) 
derived  from  other  nutrient  sources,  is  then  condensed  into 
M2  citrate.  Citrate  is  decarboxylated  into  a-ketoglutarate, 
and  then  into  additional  TCA  cycle  metabolites.  The 
abundance  of  each  isotopomer  metabolite  was  used  to 
illustrate  the  glucose  contribution  to  TCA  cycle  fluxes  in 
ovarian  cancer  cells.  Figure  3H  showed  that  complete  loss 
of  HSulf- 1  resulted  in  decreased  glucose  contribution  to 
the  majority  of  TCA  cycle  metabolites  (citrate,  fumarate, 
malate, aspartate)  (Figure  S4).  The  M2  citrate  level, 
directly  reflecting  glucose’s  entry  flux  into  TCA  cycle, 
decreased  from  0.58nmole  per  million  of  control  cells 
to  around  0.3nmole  per  million  for  both  Shi  and  Sh2 
cells  (Figure  3H).  Similar  trends  were  also  observed 
with  other  metabolites  including  fumarate,  malate,  and 
aspartate.  However,  the  I3C -labeled  metabolites  levels 
of  a-ketoglutarate,  2-hydroxyglutarate,  and  glutamate 
were  unchanged.  Collectively,  these  results  suggest  that 
HSulf-1  loss  reduced  the  glucose  contribution  to  TCA  in 
mitochondria. 

Loss  of  HSulf-1  is  associated  with  reduced 
oxidative  phosphorylation  (OXPHOS)  and 
impaired  mitochondria 

In  order  to  investigate  if  the  increased  glycolytic 
flux  and  altered  state  of  TCA  cycle  caused  by  the  loss 


of  HSulf-1  would  have  any  impact  on  mitochondrial 
function,  we  measured  the  oxygen  consumption  rate 
(OCR).  High  resolution  respirometry  measurements  in 
intact  cells  showed  lower  basal  OCR  in  cells  lacking 
HSulf-1  (Figure  4 A).  Shi  and  Sh2  cells  showed  reduction 
in  OCR  in  all  respiratory  states  compared  to  the  controls 
(Figures  4B-4D).  Moreover,  these  cells  showed  a  lower 
coupled  respiration  rate  compared  to  cells  expressing 
HSulf-1,  suggesting  a  regulatory  role  of  HSulf-1  in 
mitochondrial  ATP  production  (Figure  4D).  Re-expression 
of  HSulf-1  rescues  the  OCR  to  a  level  closer  to  NTC  cells, 
suggesting  that  HSulf-1  might  regulate  mitochondrial 
oxidative  phosphorylation.  Moreover,  Transmission 
Electron  Microscopy  (TEM)  images  clearly  showed 
morphological  alterations  (Figure  4E)  in  mitochondria, 
reduced  mitochondrial  perimeter  (Figure  4F),  reduced 
area  (Figure  4G)  and  less  functional  mitochondria  as 
evaluated  by  mitotracker  staining  with  loss  of  HSulf-1  in 
Shi  (Figure  4H). 

Conversely,  enhanced  expression  of  HSulf-1 
in  TOV21G  Cll  clearly  showed  higher  OCR  in  all 
respirometry  states  compared  to  its  vector  (Figure  41). 
The  TOV21G  Cll  cells  also  recovered  the  morphological 
features  (Figures  4J-4L)  and  functionality  of  mitochondria 
(as  evidenced  by  mitotracker  staining  in  TOV21G  vector 
transfected  compared  to  HSulf-1  transfected  cells  (Figure 
4M).  Consistent  with  these  findings,  HSulf-1  KO  MEFs 
had  lower  OCR  compared  to  WT  MEF  cells  (Figure  S5), 
suggesting  a  direct  relationship  between  HSulf-1  loss  and 
lower  OXPHOS  and  mitochondrial  alterations. 

HSulf-1  deficiency  altered  glucose  metabolism 
through  c-Myc  induction 

Our  gene  microarray  and  ingenuity  pathway 
analysis  in  NTC  and  Shi  cells  [6]  demonstrated  that 
the  genes  in  Ras-Raf-MAPK/ERK  pathway  were 
expressed  at  a  much  higher  levels  in  the  Shi  and-2  cells, 
suggesting  a  major  involvement  of  Ras-Raf-MAPK/ERK 
pathway  in  the  survival/  growth  of  cells  in  the  absence 
of  HSulf-1  (Figure  5A).  Validation  of  microarray  results 
with  immunoblot  analysis  revealed  a  robust  increase  in 
ERK  phosphorylation  and  c-Myc  expression  in  Shi  cells 
(Figure  5B).  More  importantly,  HS-mimetic  PG545, 
EGFR  inhibitor  lapatinib  and  MEK  inhibitor  U0126 
attenuated  ERK  phosphorylation  and  c-Myc  expression 
(Figure  5C),  indicating  an  involvement  of  Raf-MAPK/ 
ERK  mediated  c-Myc  activation.  Furthermore,  c-Myc 
mRNA  expression  was  2.5  fold  higher  with  enhanced 
cytoplasmic  and  nuclear  localization  of  c-Myc  in  Shi  cells 
compared  to  the  NTC  cells  respectively  (Figures  5D-5E). 
To  determine  if  c-Myc  could  potentially  play  a  role  in  this 
metabolic  alteration,  we  used  a  specific  c-Myc  inhibitor 
(10058-F)  and  observed  a  dose  dependent  decrease  in 
Glut  1,  HKI1  and  LDHA  expression  in  Shi  cells  (Figure 
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Figure  4:  HSulf-1  is  associated  with  reduced  OXPHOS  and  impaired  mitochondria.  A.  Basal  oxygen  consumption  in  OV202 
clonal  cells.  Oligomycin,  FCCP  and  Rotenone  were  added  subsequently  to  measure  leak  B.,  maximal  C.  and  coupled  respirations  D.  E. 
Cultured  NTC  and  Shi  cells  were  photographed  following  TEM  at  5000X  and  in  inset  25000X  magnification  to  show  altered  mitochondrial 
morphology.  Mitochondrial  perimeter  F.  and  area  G.  were  calculated  using  Image  J  software  from  TEM  micrographs  in  10  cells  from  three 
different  micrographs,  *p  <  0.05.  H.  Mitotracker  staining  of  OV202  clonal  cells.  I.  Basal,  leak,  maximal  and  coupled  respiration  rates 
in  TOV21G  Vec  and  Cll  cells.  J.  Electron  micrographs  of  TOV21G  Vec  and  Cll  cells  with  5000X  and  in  inset  25000X  magnification. 
Mitochondrial  perimeter  K.  and  area  L.  were  calculated  in  TOV21G  Vec  and  Cl  using  Image  J  software  from  TEM  micrographs  in  10  cells 
from  three  different  micrographs,  *p  <  0.05.  M.  Mitotracker  staining  of  TOV21G  Vec  and  Cl  cells. 
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5F,  panel  2).  However,  Glutl  and  HKII  levels  were 
downregulated  only  at  the  highest  concentration  of  the 
c-Myc  inhibitor  (lOOpM)  in  the  NTC  cells  with  minimal 
changes  in  LDHA  levels.  To  correlate  whether  HB-EGF 
mediated  enhanced  growth  factor  signaling  in  Shi  cells 
results  in  higher  c-Myc  levels,  we  exogenously  added 
HB-EGF  and  showed  a  time  dependent  increase  of  c-Myc 


expression  along  with  augmented  Glut  1,  HKII  and  LDHA 
expression  (Figure  5G).  Pre-treatment  of  10058-F  reduced 
the  HB-EGF  mediated  increment  of  Glut  1,  HKII  and 
LDHA  levels  (Figure  5H).  These  data  suggests  that  loss 
of  HSulf-1  can  has  a  more  pronounced  effect  in  activating 
c-Myc  through  HB-EGF  mediated  ERK  activation  and 
subsequently  reprograms  the  energy  metabolism. 
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Figure  5:  c-Myc  is  important  for  driving  glycolysis  in  HSulf-l-deficient  cells.  A.  Normalized  unsupervised  hierarchical 
clustering  of  Wnt,  ERK  and  EGF  pathway  genes  for  NTC  and  Shi,  each  class  representing  three  biological  replicates.  B.  Immunoblot 
analysis  of  p-ERK,  ERK  and  c-Myc  in  NTC  and  Shi.  C.  Cell  lysate  from  24  hr  treatment  with  PG545  (20pM),  Lapatinib  (20nM)  and 
U1026  (20pM)  were  immunoblotted  against  p-ERK.  ERK  and  c-Myc.  D.  QRT-PCR  analysis  of  c-Myc  in  NTC  and  Shi  cells,  *P  <  0.001 . 
E.  NTC  and  Sh  1  cells  were  immuno-stained  with  anti-c-Myc  antibody  (green)  and  mounted  with  D  API  and  photographed  using  Zeiss  LSM 
510  confocal  microscope.  OV202NTC  and  Shi  cells  were  treated  with  increasing  concentration  of  c-Myc  inhibitor  for  24  hr  and  F.  with 
HB-EGF  for  indicated  time  and  immunoblotted  G.  for  c-Myc,  Glut  1,  HKII  and  LDHA  (N=  2).  H.  Shi  cells  were  serum  starved  for  2  hr, 
treated  with  c-Myc  inhibitor  (10058-F4)  for  another  2  hr  and  then  treated  with  HB-EGF  for  indicated  time  points.  Protein  lysates  from  each 
condition  were  probed  against  Glut  1,  HKII  and  LDHA,  where  GAPDH  served  as  loading  control. 
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Tumor  development  and  metabolic  alteration 
by  loss  of  HSulf-1  is  attenuated  by  HS  mimetic 
PG545 

To  better  understand  the  role  of  the  putative  tumor 
suppressor  HSulf-1  in  HS-dependent  growth  factor 
signaling,  and  the  subsequent  metabolic  reprogramming 
and  pro-tumorigenicity,  we  treated  Shi  cells  with  PG545, 
which  mimics  the  action  of  HSulf-1.  PG545  treatment 
significantly  reduced  glucose  uptake  (Figures  6A-6B)  and 
lactate  production  (Figure  6C)  in  these  cells.  Moreover, 
PG545  dose-dependently  reduced  the  phosphorylation 
of  ERK  and  the  levels  of  Glut  1 ,  HKII,  and  LDH A  more 
in  NTC  cells  expressing  HSulf-1  compared  to  Shi  cells 
(Figures  6D  and  6E  respectively).  Densitometric  analysis 
of  the  fold  change  in  the  levels  of  cMyc,  Glutl  and  HKII 
in  NTC  and  Shi  cells  is  shown  as  bar  graph  following 
PG545  treatment  in  Figures  6F  and  6G.  We  postulated  that 
PG545  might  attenuate  glucose  uptake  in  vivo,  as  PG545 
treatment  decreased  c-Myc  activation  through  ERK 
inhibition  reduced  glucose  uptake  and  Glut  1  expression 
in  vitro.  18FDG-PET  imaging  in  control  and  PG545  treated 
Shi -xenografts  showed  a  50%  reduction  in  18FDG  uptake 
in  PG545  treated  mice  (Figures  6H  and  61).  These  findings 
demonstrated  that  PG545  selectively  suppresses  glucose 
uptake  in  vivo. 

In  a  parallel  experiment,  PG545  treatment  markedly 
decreased  the  size  and  number  of  metastatic  nodules 
compared  to  control  mice  (Shown  as  inset  in  Figure  6J). 
The  mean  excised  tumor  weight  was  reduced  by  65.04% 
in  PG545  treated  mice  compared  to  control  (Figure  6J). 
Moreover,  the  abdominal  circumference  which  reflects 
the  tumor  burden  of  the  peritoneum  was  significantly 
less  in  the  treated  mice  compared  with  untreated  mice 
(Figure  6K).  To  test  whether  PG545  abrogated  the  HB- 
EGF  mediated  signaling  and  subsequently  enhanced 
glycolysis,  we  determined  the  expression  of  the  glycolytic 
enzymes,  p-ERK  and  cMyc  in  the  xenografts  of  control 
untreated  and  PG545  treated  mice.  Western  blot  analysis 
showed  that  PG545  significantly  reduced  the  expression  of 
p-ERK,  c-Myc  and  the  levels  of  Glut  1,  HKII  and  LDHA 
in  the  tumor  lysates  (Figure  6L).  Additionally,  increase  in 
cleaved  PARR  and  cleaved  capase3  in  PG545  xenografts 
indicate  that  the  reduction  in  the  tumor  growth  could  be 
the  result  of  increased  apoptosis  in  these  cells  (Figure  6M). 

DISCUSSION 

Increasing  evidence  suggest  that  the  Warburg 
effect  [25,  26]  is  common  in  cancer  and  is  driven  by 
oncogene  addiction/tumor  suppressor  loss  or  mutations 
in  the  mitochondrial  enzymes  [27].  In  this  study,  we 
report  that  HSulf- 1  loss  promotes  glycolysis  and  impairs 
mitochondrial  function  leading  to  a  significant  reduction 
in  OXPHOS.  This  phenotype  is  coupled  to  increased 


growth  factor  mediated  ERK  activation  resulting  in  the 
upregulation  of  c-Myc  activation.  The  frequent  loss  of 
HSulf-1  in  OVCA  may  thus  represent  a  genetic  alteration 
promoting  metabolic  phenotype  of  OVCA. 

Using  HSulf-1  ShRNA  downregulated  OV202 
cells  and  a  clonal  line  of  TOV21G  with  enhanced 
HSulf- 1  expression,  we  clearly  showed  modulation 
of  the  glycolytic  pathway  both  in  vitro  and  in  vivo. 
Additionally,  transient  overexpression  of  HSulf-1  in 
prostate,  lung  and  breast  cancer  cell  lines  reverses  the 
glycolytic  phenotype,  suggesting  a  critical  role  of  HSulf- 1 
in  metabolic  rearrangement  in  cancer.  Glycolysis  and 
mitochondrial  respiration  are  tightly  coupled  processes. 
A  key  branch  point  in  the  glycolytic  pathway  is  the 
production  of  pyruvate,  which  under  anaerobic  conditions 
is  metabolized  to  lactate  by  LDH  [28]  and  in  normoxia 
by  PDH  to  form  acetyl-CoA  [23].  The  conversion  of 
pyruvate  to  acetyl  CoA  by  PDH  is  a  crucial  event  which 
directs  the  energy  flow  from  glycolysis  to  the  TCA  cycle 
[24].  Phosphorylation  in  the  El  complex  of  PDH  by  PDK 
inhibits  its  activity  and  energy  flow  to  the  mitochondria. 
In  most  cancers,  LDHA  is  highly  expressed  and  diverts  the 
energy  flow  from  mitochondria  to  from  lactate.  Indeed,  our 
results  demonstrated  that  silencing  of  HSulf-1  markedly 
enhanced  LDHA  level  and  activity  in  Shi  and  HSulf-1  KO 
MEF  cells,  whereas  re -expression  of  HSulf-1  in  Shi  cells 
C17  resulted  in  the  reduction  of  LDHA  activity.  Moreover, 
direct  evidence  from  isotope  tracer  analysis  confirmed  less 
glucose  entry  into  the  TCA  cycle.  Together,  these  results 
demonstrate  that  HSulf-1  loss  hindered  the  energy  flow 
from  glycolysis  into  TCA  cycle. 

Of  note,  it  was  originally  hypothesized  that  these 
metabolic  changes  reflect  a  non-functional  state  of 
mitochondrial  oxidative  phosphorylation  OXPHOS 
and  that  the  high  glycolytic  rate  in  the  tumors  is  due  to 
impaired  mitochondrial  respiration  [29].  Interestingly, 
our  results  showed  that  stable  knockdown  of  HSulf-1 
exhibited  lower  oxidative  phosphorylation  compared 
to  NTC,  while  overexpression  of  HSulf-1  reversed  this 
phenotype  in  TOV21G  Cll  cells.  Because  tumors  arise 
through  a  Darwinian  evolutionary  process,  the  continuous 
absence  of  OXPHOS  in  tumor  cells  could  dictate  one 
or  multiple  reversible  and/or  irreversible  adaptations 
that  confer  a  growth  advantage  [30].  It  is  likely  that  the 
reduction  of  OXPHOS  capacity  of  these  tumor  cell  lines 
had  been  irreversibly  diminished  and  could  reflect  on 
mitochondrial  biogenesis  and  morphology  as  previously 
described  [30].  Indeed,  our  results  demonstrated  that 
there  is  lower  mitochondrial  activity  in  HSulf-1 -deficient 
cells  which  can  be  rescued  by  re-expression  of  HSulf-1. 
Interestingly,  TEM  analysis  demonstrated  a  reduction 
in  mitochondrial  number  and  altered  mitochondrial 
morphology  (extended  to  round  shape)  in  HSulf- 1  deleted 
Shi  cells.  Alternatively,  HSulf-1  overexpression  in 
TOV21G  cells  resulted  in  a  decrease  in  the  mitochondrial 
activity  and  rescue  of  normal  mitochondrial  structure. 
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Figure  6:  PG545  diminished  HSulf-l-deficiency  induced  glycolysis  and  tumor  growth.  A.  Glucose  uptake  in  live  Shlcells 
after  PG545  (15  pM)  treatment  for  1  hr.  B.  Fluorescent  intensity  of  glucose  uptake  in  control  and  PG545  treated  Shi  cells.  C.  Lactate 
secretion  after  12  and  24  hr  PG545  treatment  in  Shi  cells.  D.  and  E.  Immunoblot  analysis  of  OV202  NTC  and  Shi  cells  treated  with  PG545 
(0-20pM)  for  24  hr.  F.  and  G.  Augmented  downregulation  of  c-Myc,  Glutl  and  HKII  in  OV202  NTC  cells  compared  to  Shi  cells  following 
PG545  treatment  shown  as  fold  change  respectively.  H.  Micro-PET  imaging  of  18FDG  uptake  in  control  and  PG545  (20mg/kg,  i.p.)  treated 
mice.  A  representative  image  is  provided  in  the  transverse,  coronal,  and  sagittal  regions  of  interest  in  the  tumor.  I.18FDG  uptake  intensities 
were  calculated  in  at  least  20  planes  and  analyzed  using  PMOD  Software,  *P  <  0.01.  G.  A  separate  experiment  with  randomized  tumor¬ 
bearing  mice  ( n  =  10),  were  treated  with  water  or  PG545  (20mg/kg)  twice  a  week  and  after  28  days  mice  were  euthanized  and  imaged.  J. 
Excised  Tumor  weight  from  control  and  PG545  treated  mice,  *P  <  0.00 1 .  K.  Abdominal  circumference  from  untreated  and  PG545  treated 
mice  ( n  =  10)  measured  at  the  day  of  sacrifice. *P  <  0.001.  L.-M.  Tumor  lysates  from  control  (n  =  5)  and  PG545  treatment  (n  =  5)  were 
selected  randomly  and  immunoblotted. 
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These  results  demonstrated  a  fine  interplay  between 
glycolysis  and  mitochondrial  metabolism  regulated  by 
HSulf-1.  Collectively,  our  findings  clearly  established 
a  strong  connection  between  loss  of  HSulf- 1  and 
reprogramming  of  metabolism,  by  inducing  glycolysis, 
disrupting  the  energy  flow  to  mitochondria  and  affecting 
mitochondrial  morphology  and  activity. 

These  metabolic  changes  appear  to  be  strongly 
associated  with  HSulf- 1-HB-EGF  mediated  signaling 
resulting  in  c-Myc  activation.  Studies  suggests  that 
the  metabolically  active  pathways  are  regulated  by  an 
upstream  cascade  of  signaling  mechanisms  involving 
growth  factor  receptor  tyrosine  kinases  (RTK)  and 
this  involvement  of  RTK  is  implicated  in  glucose 
homeostasis  [31].  Previous  reports  from  our  lab  and 
others  unequivocally  established  the  role  of  HSulf- 1  in 
modulating  heparin-binding  growth  factor  signaling 
leading  to  enhanced  ERK  activation  [6].  Consistent  with 
these  observations  our  microarray  analysis  and  in  vitro 
results  showed  activation  of  EGFR-Ras-Raf-MAPK  axis 
which  culminates  in  ERK-mediated  c-Myc  activation. 
c-Myc  has  been  established  as  a  master  regulator  of 
metabolism  and  directly  regulates  the  glycolytic  pathway 
by  activating  Glut  1,  PFK,  PKM2,  LDHA  and  enolase 
[32].  Additionally,  c-Myc  has  been  found  to  regulate 
mitochondrial  function,  biogenesis  and  morphology  [33]. 
Indeed,  our  results  showed  that  HB-EGF  stimulation 
resulted  in  a  time-dependent  induction  of  c-Myc  both 
in  NTC  and  Shi  cells.  However,  HB-EGF  treatment 
upregulated  the  glycolytic  enzymes  more  in  the  Shi  cells 
compared  to  the  NTC  cells.  Importantly,  c-Myc  inhibition 
markedly  reduced  the  basal  and  HB-EGF -induced  Glut 
1,  HK-II  and  LDHA  levels  in  Shi  cells.  These  results 
provide  a  plausible  explanation  of  how  loss  of  HSulf- 1 
mediated  enhanced  HB-EGF-growth  factor  signaling 
can  activate  c-Myc  and  as  a  consequence,  directly  alter 
glucose  metabolism.  However,  we  noticed  that  HB-EGF 
enhanced  the  cMyc  target  genes  at  early  time  points  and 
tapered  at  90  mins,  while  sustaining  c-Myc  expression. 
This  possibly  could  be  due  to  two  reasons.  It  has  been 
reported  that  growth  factor  stimulation  results  in  biphasic 
effects  leading  to  sudden  increase  in  the  target  gene 
expression  followed  by  a  tapering  but  sustained  effect  [34, 
35].  Second,  it  is  possible  that  c-myc  acts  in  conjunction 
with  other  transcriptional  complex  factor  in  regulating 
these  genes.  Thus  availability  of  all  the  factors  is  essential 
for  expression  of  Glut  1,  HKII  and  LDHA.  Limited 
availability  of  other  factors  and/or  shorter  half-  life  of  the 
factors  involved  may  be  a  reason  for  decreased  expression 
of  target  genes  at  later  time  points. 

Metabolic  targeting  for  cancer  therapy  is  currently 
under  investigation  and  small  molecules  specifically 
inhibiting  key  glycolytic  steps  have  been  evaluated  as 
an  anti-cancer  approach.  However,  therapeutic  agents 
that  target  glycolysis  do  not  always  provide  satisfying 
outcomes  in  clinical  trials  [36].  We  report  here  that 


PG545  attenuated  glucose  uptake  and  lactate  production 
through  the  inhibition  of  p-ERK,  c-Myc  and  the  key 
glycolytic  enzymes  in  vitro  and  in  vivo.  Moreover,  PG545 
decreased  tumor  growth  possibly  by  inducing  apoptosis 
as  evidenced  by  increased  cleaved  Caspase  3  and  PARP. 
More  intriguing  is  the  possibility  that  PG545,that  mimics 
the  action  of  HSulf- 1  by  competing  for  growth  factor 
binding  and  thus  inhibiting  HB-EGF  mediated  signaling 
has  additional  effects  in  vivo  by  targeting  the  tumor 
microenvironment  [6].  It  is  noteworthy  that  PG545 
treatment  of  the  NTC  cells  expressing  HSulf- 1  had  a  more 
pronounced  effect  in  downregulating  the  expression  of 
the  Glutl,  HKII  and  LDHA,  thus  supporting  the  role  of 
HSulf-1  as  a  modulator  of  glycolysis.  Our  study,  therefore, 
has  important  implications  for  the  development  of  ovarian 
cancer  therapies  targeting  the  altered  metabolism,  where 
tumor  suppressor  HSulf-1  is  lost. 

MATERIALS  AND  METHODS 


Materials 

PG545  supplied  by  Progen  Pharmaceuticals  Ltd 
(Brisbane,  Australia). Antibodies  used:  HSulf-1,  PFKFB3 
(Abeam,  MA,  USA),  Glutl  (Santa  Cruz  Biotechnology, 
CA,  USA),  HKII,  ALDOC,  PGAM,  PKM2,  LDHA, 
PDH,  PDK1,  Tubulin,  GAPDH  (GeneTex,  CA,  USA),c- 
Myc,  p-ERK,  T-ERK  (Cell  signal  technology,  MA,  USA). 
Secondary  antibodies  anti-HRP-Rabbit,  anti-HRP-Mouse 
antibodies  were  from  BD  Pharmingen  (CA,  USA).  ATP 
and  LDH  assay  kits  were  from  Biovision  (CA,  USA)  and 
Sigma-Aldrich  (MO,  USA). 

Cell  lines  and  culture  condition 

H358  is  a  non-small  cell  lung  carcinoma  cell 
line  derived  from  metastatic  sites  (http://www.phe 
culturecollections.org.uk/products/celllines/generalcell/ 
detail  .jsp?refld=95111733&collection=ecacc_gc ),  highly 
metastatic  prostate  cancer  cell  line  PC3  [37]and  HSulf-1 
deficient,MDA231,  a  triple  negative  breast  cancer  cell 
line  [15]  were  cultured  in  10%  RPMI-1640  and  DMEM 
respectively.  HSulf- 1+,+ and  HSulf- l'1' cells  were  isolated 
from  HSulf-l+/+  and  HSulf- 1  mice  and  immortalized 
as  previously  described  [38].  MEFs  were  maintained  in 
DMEM  supplemented  with  10%  FBS  and  650  pg/ml  of 
G418.  H358,  PC3  and  MDA231  are  HSulf-1  deficient 
cells. 

Transfection,  plasmids  and  shRNA 

Transient  Transfection:  pcDNA-HSulf- 1  plasmid 
was  transiently  overexpressed  in  H358,  PC3  and  MDA231 
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cells  using  Lipofectamine  Plus  (Invitrogen,  Grand  Island, 
NY)  according  to  the  manufacturer’s  instruction. 

Generation  of  HSulf-1  overexpression  stable  clones 

Upon  60%  confluency,  HSulf-1  deficient  ovarian 
clear  cell  cancer  line  TOV21G  cells  [10]  were  transfected 
with  pcDNA3.1-WT-HSulf-l  using  lipofectamine-LTX- 
plus.  After  24hr  of  transfection,  cells  were  selected  on 
1  OOug/m  I  G418  and  individual  clones  were  picked  using 
propagated  as  HSulf-1  overexpressing  TOV21G  stable 
clones. 

Generation  of  HSulf-1  downregulated  Stable  clone 

OV202  stable  clones  OV202Shl  and  OV202Sh2 
were  generated  as  previously  reported  [16]. 

Rescue  of  HSulf-1  in  Shi  cells 

As  HSulf-1  ShRNA  targets  the  3’UTR  of  HSulf-1 
[16],  HSulf-1  expression  in  Shi  cells  was  rescued  with 
CMV-driven  WT-HSulf-1  construct  in  OV202Shlcells 
We  used  OV202Shl-clone-7  as  our  rescue  model  and 
designated  as  C17  throughout  the  study.  OV202  clones 
were  cultured  in  minimal  essential  medium  supplemented 
with  20%  FBS  with  proper  antibiotic  selection. 

Glucose  uptake 

Glucose  uptake  of  the  live  cells  was  measured  using 
2-NBDG  (Cayman  Chemicals,  Michigan,  USA)  according 
to  manufacturer  instruction.  Fluorescent  intensities  were 
calculated  using  Image  J  software. 

Immunoblot  analysis 

lmmunoblot  analysis  was  carried  out  as  previously 
described  [39]. 

ATP,  LDH  assay 

Cellular  ATP  content  and  LDH  activity  were 
measured  using  kits  from  Bio  Vision  and  calculated  as  per 
manufacturer’s  instruction. 

Quantitative  real  time  PCR 

Quantitative  real-time  PCR  was  carried  out  with 
SYBR-Green  PCR  Master  Mix  (Applied  Biosystems, 
NY,  USA),  using  specific  primers  with  GAPDH  or  18S 
ribosomal  subunit  as  an  internal  control  as  previously 
described  [21]. 


Mitotracker  staining  and  confocal  microscopy 

Cells  (5xl04/500pl)  were  seeded  in  a  4-well 
chamber  slide,  incubated  overnight  and  stained  with 
lOOnM  mitotracker  red  for  30  minutes.  Cells  were  then 
washed,  fixed  and  mounted  with  DAPI  and  analyzed  under 
Zeiss  LSM  510  confocal  microscope.  Intensity  of  red 
florescence  was  measured  using  Image  J  software. 

Transmission  electron  microscopy 

Cultured  cells  were  washed  and  fixed  in 
Trumps  fixative  containing  4%  formaldehyde  and  1% 
glutaraldehyde  in  a  phosphate  buffer  pH  -7.3,  post- 
fixed  in  1 .0%  Os04,  dehydrated  with  ethanol  gradation, 
and  transitioned  into  propylene  oxide  for  infiltration  and 
embedding  into  spurr  epoxy  resin.  Ultrathin  sections  were 
cut  onto  grids,  stained  with  uranyl  acetate  and  lead  citrate, 
and  examined  with  JEM- 1400  (JEOL  USA)  transmission 
electron  microscope  and  digitally  photographed. 

Microarray  expression  data  analysis 

OV202NTC,  Shi  and  Sh2  cells  in  triplicates  were 
profiled  using  Illumina  Whole  Genome  DASL  assay  as 
previously  described  [21]. 

Liquid  chromatography/Mass  Spectrometry  (LC/ 
MS,  LC/MS2)  and  Gas  chromatography/Mass 
Spectrometry  (GC/MS) 

The  samples  destined  for  LC/MS,  LC/MS2  and  GC/ 
MS  analysis  were  prepared  and  analyzed  as  previously 
described  [21]. 

Metabolic  tracing  using  GC-MS 

Metabolic  extraction 

At  -80%  confluence,  OV202  cells  (NTC,  Shi, 
Sh2)  were  cultured  with  U-13C6  glucose,  then  extracted, 
derivatized  and  analyzed  by  GC/MS  as  previously 
described  [40,  41]. 

Animal  study 

Both  the  animal  experiments  were  performed  with 
the  approved  protocol  of  Mayo  Clinic  Institutional  Animal 
Care  and  Use  Committee  (IACUC). 

FDG-PET  imaging 

Female  athymic  mice  (5-6  weeks)  were  injected 
intraperitoneally  (i.p.)  with  5X106  OV202  NTC  ( n  =  5) 
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and  OV202  Shi  ( n  =  5)  cells.  After  2  weeks,  mice  were 
anesthetized  by  isoflurane  and  500  pCi  of  18FDG  was 
administered  intravenously  (i.v.)  via  the  tail  vein.  After 
45  minutes  of  18FDG  injection,  micro-PET  imaging  was 
acquired  using  Siemens  Inveon  PET/CT  scanner  (Siemens 
Medical  Solutions,  Molecular  Imaging,  TN,  USA). 
Regions  of  interest  in  the  tumor  were  quantified  using 
the  PMOD  Biomedical  Image  Quantification  and  Kinetic 
Modeling  Software  (PMOD  Technologies,  Switzerland). 

Tumor  xenograft 

Female  athymic  mice  (5-6  weeks)  were  randomized 
in  two  groups  {n  =  10)  and  OV202  Shi  cells  (4X106) 
were  injected  i.p.  PG545  treatment  commenced  5  days 
after  tumor  inoculation  at  20mg/kg,  i.p.,  twice  a  week. 
Treatment  continued  for  the  next  4  weeks.  However 
control  mice  were  sacrificed  during  week  3  as  the  tumor 
burden  exceeded  10%  of  their  body  weight.  Mice  were 
sacrificed  according  to  standard  1ACUC  procedure  and  the 
tumors  and  tissues  were  excised  and  preserved  either  in 
formalin  or  -80°C. 

Statistical  analysis 

All  the  results  were  expressed  as  the  mean  ±  S.D. 
data  obtained  from  three  separate  experiments.  All 
statistical  analysis  was  evaluated  using  Graph  pad  Prism 
software  (San  Diego).  Data  were  analyzed  by  the  paired  t 
test,  and  P  values  less  than  0.05  or  mentioned  otherwise, 
was  considered  statistically  significant. 
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Abstract 


As  defective  autophagy  and  deranged  metabolic  pathways  are  frequently  observed  in 
cancer,  simultaneous  pharmacological  targeting  of  these  two  pathways  could  provide  a 
viable  therapeutic  option.  These  pathways  have  the  ability  to  sense  nutrient  deprivation 
and  thereby  help  to  conserve  and  maintain  the  supply  of  cellular  energy  rich  molecules. 
However,  it  is  unknown  how  these  pathways  are  regulated  by  limited  availability  of 
extracellular  growth  factors.  Here  we  show  that  an  endosulfatase,  HSulf-1,  a  known 
tumor  suppressor  which  attenuates  heparan  sulfate  binding  growth  factor  signaling  also 
regulates  the  interplay  between  autophagy  and  lipid  biogenesis.  Silencing  of  HSulf-1  in 
ovarian  cancer  cell  lines  OV202  and  TOV2223  resulted  in  suppressed  autophagic  flux, 
LC3B  puncta  and  autophagosome  formation  and  increased  lipid  droplet  formation.  In 
contrast  HSulf-1  proficient  cells  exhibited  significantly  more  autophagic  vacuoles  (AVs) 
and  reduced  number  of  lipid  droplets  (LDs).  Increased  accumulation  of  LDs  in  HSulf-1 
depleted  cells  was  associated  with  increased  ERK  mediated  cPLA2S505 
phosphorylation,  a  well-known  marker  for  lipid  biogenesis.  Conversely,  ectopic 
expression  of  HSulf-1  in  SKOV3  cells  reduced  the  number  of  LDs  and  increased  the 
numbers  of  AVs  compared  to  vector  transfected  controls.  Furthermore,  pharmacological 
(AACOCF3)  and  shRNA-mediated  downregulation  of  cPLA2  resulted  in  reduced  LDs 
confirming  the  role  of  cPLA2  in  LD  biogenesis  and  increased  autophagy.  Using  mouse 
model  of  OV202  Shi  derived  xenograft,  we  show  that  AACOCF3  effectively  attenuated 
tumor  growth,  LD  biogenesis  and  reduced  tumorigenesis  in  vivo.  Collectively,  these 
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results  show  a  reciprocal  regulation  of  autophagy  and  lipid  biogenesis  by  HSulf-1  in 
ovarian  cancer. 

Introduction 

Previous  reports  have  shown  that  downregulation  of  HSulf-1  is  common  in 
ovarian  cancer  (OvCa)  and  regulates  heparan  sulfate  binding  growth  factor  signaling 
which  subsequently  promotes  tumorigenesis.1  We  recently  reported  that  loss  of  HSulf-1 
promotes  a  “lipogenic  phenotype”  as  evidenced  by  an  increase  in  lipid  related 
metabolites,  fatty  acid  synthesis  and  beta-oxidation  indicating  an  important  role  of 
HSulf-1  in  metabolic  regulation.2 

Although  adipocytes  were  described  as  the  primary  site  for  LD  biogenesis,  3’  4 
recent  findings  suggest  lipid  droplets  as  an  important  source  of  energy  in  cancer  cells.5'7 
Enhanced  LD  biogenesis  in  cancer  cells  signifies  its  role  in  cell  signaling,  membrane 
trafficking  and  lipid  metabolism,  all  associated  with  increased  survival  and  growth  of 
cancer  cells.8  LDs  are  considered  as  cellular  hallmarks  of  many  different  diseases  such 
as  diabetes,  atherosclerosis  and  cancer.8'12  Though,  the  presence  of  LDs  is  associated 
with  disease  progression,  the  functional  significance  in  promoting  inflammation  and 
tumorigenesis  is  not  well  understood.  More  importantly,  the  molecular  alterations  that 
promote  LD  accumulation  in  cancer  cells  have  not  been  defined. 

Primarily,  LDs  are  storage  organelles  for  neutral  lipids  and  cholesterol  esters.13 
Stress  induced  release  of  fatty  acids  from  the  stored  LDs  in  adipocytes  provide  energy 
which  subsequently  promote  tumor  growth,  metastasis  and  cell  survival  of  OvCa.14 
Several  of  the  LD  associated  proteins  SREBP1,  PLINs,  PLA2G4A,  PLA2G3,  ATGL, 
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HSL,  MAGL,  PPARy  involved  in  LD  biogenesis  and  release  of  fatty  acids  are 
upregulated  upon  HSulf-1  loss.3, 15-20  Many  of  these  genes  are  overexpressed  in  cancer 
cells  and  implicated  as  potential  contributors  towards  tumorigenesis.3, 21-24 

Autophagy  is  a  lysosomal  pathway  by  which  long  lived  proteins  and  damaged 
organelles  are  degraded  to  provide  energy  to  the  cell  and  maintain  cellular 
homeostasis.25-27  Although  it  is  generally  accepted  that  LDs  are  substrates  for  lipases, 
Singh  et  al  28,  29  identified  lipid  droplets  as  substrates  for  macroautophagy  and  coined 
the  term  lipophagy  to  describe  the  engulfment  of  LDs  by  autophagosomes  that 
ultimately  fuse  with  lysosomes  for  the  breakdown  of  LD  components  under  stress 
conditions.  Following  this  seminal  finding,  there  are  several  reports  on  the  regulation  of 
LDs  by  autophagy.30-32  However,  the  interplay  and  the  molecular  players  involved  in 
these  two  metabolic  pathways  are  not  known  in  the  context  of  OvCa.  In  the  present 
work,  we  identify  HSulf-1,  a  major  regulator  of  growth  factor  signaling  as  a  missing  link 
between  autophagy  and  lipid  droplets  in  OvCa  acting  through  cPLA2a  as  one  of  the 
molecular  pathways  to  promote  LDs  and  inhibit  autophagy. 

Results 

HSulf-1  loss  promotes  lipid  droplet  biogenesis  and  defective  autophagy 

We  previously  reported  that  stable  genetic  downregulation  of  HSulf-1  in  OV202 
cells  achieved  by  two  different  shRNAs  targeting  HSulf-1  (OV202Sh1  and  Sh2  cells) 
displayed  increased  LDs  compared  to  non-targeted  control  transduced  cells  (NTC)2.  We 
further  confirmed  these  findings  using  OV202Sh1  to  OV202Sh1  Clone  Cl 7  cells 
(generated  by  re-expression  of  CMV  driven  HSulf-1  construct  in  Shi  cells.)  HSulf-1 
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levels  were  verified  in  these  clones  by  western  blot  analysis  as  shown  in  Fig  1A. 
Subsequently,  lipid  droplets  were  monitored  by  Bodipy  staining  in  these  cells.  Rescue  of 
HSulf-1  in  OV202Sh1CI7  cells  reduced  the  number  of  LDs  compared  to  OV202Sh1  and 
Sh2  cells  (Fig. IB,  Insets).  Mean  fluorescence  intensity  of  Bodipy  staining  in  these  cells 
is  shown  in  Fig.  1C. 

We  further  performed  quantification  of  the  changes  in  the  LDs  by  Transmission 
Electron  Microscopy  (TEM).  Data  showed  significantly  higher  numbers  of  LDs  in 
OV202Sh1  and  Sh2  cells  compared  to  NTC  cells  (Figs.  ID  and  E).  The  increased  lipid 
droplets  indicate  the  nutrient  rich  or  energy  proficient  state  of  the  cells  and  therefore  we 
hypothesized  that  cells  will  not  activate  catabolic  process  such  as  autophagy. 
Therefore,  to  verify,  we  determined  the  extent  of  autophagy  via  TEM  analysis  in  these 
cells  and  observed  that  HSulf-1  depleted  OV202Sh1  and  Sh2  cells  exhibit  less 
autophagic  vesicles  (AV)  when  compared  to  OV202NTC  cells.  Also,  rescue  of  HSulf-1 
in  Sh1CI7  cells  resulted  in  increased  AV  compared  OV202Sh1  cells.  To  rule  out  the 
possibility  of  cell  specific  effects,  downregulation  of  HSulf-1  in  TOV2223  cells  with 
ShRNA  (Fig. IF)  also  showed  increased  numbers  of  LDs  in  TOV2223Sh1  cells 
compared  to  TOV2223NTC  cells  by  both  Bodipy  staining  (Fig.lG)  and  TEM  analysis 
(Fig.  1H)  indicating  that  these  metabolic  alterations  were  not  unique  to  the  OV202  cell 
line.  Similarly,  a  high  number  of  AVs  was  observed  in  both  OV202NTC  cells  and 
TOV2223NTC  cells  as  compared  to  corresponding  HSulf-1  depleted  cells.  The 
quantification  of  LDs  and  AVs  in  TOV2223  cells  is  shown  in  Fig.  II.  Furthermore,  HSulf- 
1  knockout  (KO)  mouse  embryonic  fibroblasts  (MEFs)  also  displayed  increased 
numbers  of  LDs  compared  to  wild-type  MEFs  (Fig. IK).  Similarly  upon  ectopic  expression 
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of  HSulf-1  in  SKOV3  cells,  there  was  a  significant  decrease  in  the  numbers  of  LDs 
compared  to  vector  transfected  controls  (Fig. SI).  In  contrast,  vector  transfected  SKOV3 
cells  exhibited  higher  degree  of  AVs  as  quantified  through  TEM  analysis,  compared  to 
HSulf-1  transfected  SKOV3  cells  (Figs  SIB  and  D  respectively).  Quantification  of  LDs 
and  AVs  in  25  cells  is  shown  in  Fig. SIC. 

cPLA2  is  activated/phosphorylated  in  HSulf-1  depleted  ovarian  cancer  cells 

Since  our  initial  observation  indicated  increased  lipid  droplet  accumulation  upon 
HSulf-1  depletion  in  ovarian  cancer  cells,  we  next  wanted  to  determine  the  underlying 
mechanisms.  LD  often  accumulates  from  anabolic  process  known  as  lipid  biosynthesis. 
Previously  we  showed  that  lipid  droplet  associated  proteins  were  upregulated  in  HSulf-1 
depleted  cells.  Among  them,  one  of  the  key  proteins  involved  in  lipid  droplet  formation  is 
cPLA2a  previously  shown  to  be  essential  for  LD  biogenesis.33  cPLA2a  is  activated  and 
phosphorylation  at  Ser505  by  p-ERK  has  been  shown  to  be  critical  for  cPLA2  mediated 
LD  biogenesis.34  Therefore,  we  first  determined  the  activated/  phosphorylated  levels  of 
cPLA2a  by  Immunoblot  analysis.  Our  data  show  that  p-cPLA2ser505  was  clearly 
increased  in  OV202Sh1  and  to  a  lesser  extent  in  OV202Sh2  cells  (Fig.  2A).  Re¬ 
expression  of  HSulf-1  also  reduced  the  cPLA2  phosphorylation,  indicating  that 
phosphorylation  was  reversed  by  increasing  levels  of  HSulf-1. 

Pharmacological  Inhibition  of  cPLA2  attenuates  LD  biogenesis  in  OV202SM  cells. 

To  further  investigate  the  role  of  cPLA2  activation  and  lipid  droplet 
accumulation/biosynthesis,  we  treated  OV202Sh1  and  OV2008  cells  with  a  cPLA2 
specific  inhibitor,  AACOCF3  (10  and  20|iM).  Western  blot  analysis  show  decreased 
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levels  of  p-cPLA2  in  inhibitor  treated  cells  (Figs.2B  and  C).  Bodipy  staining  upon 
inhibition  of  cPLA2  activity  with  lOpM  of  AACOCF3  and  MAFP  (another  cPLA2  inhibitor) 
in  OV202Sh1  and  -2  cells  (Fig.  2D)  showed  almost  complete  inhibition  of  LD  biogenesis 
in  these  cells  compared  to  untreated  control  cells.  Consistent  with  this  data,  TEM 
analysis  of  OV202Sh1  cells  with  AACCOF3  and  MAFP  also  revealed  significantly  lower 
number  of  LDs  (Fig.2E,  with  quantitation  of  LDs  in  25  cellsshown  next  to  the 
representative  TEM  images).  Furthermore,  transient  downregulation  of  cPLA2 
expression  with  two  different  ShRNA  (ShcP-1  and  ShcP-2)  against  cPLA2  (Fig.  2F)  in 
OV202Sh1  cells  resulted  in  the  decrease  of  LDs  as  revealed  by  Bodipy  staining  (Fig. 
2G).  Similarly,  stable  downregulation  of  cPLA2  in  OV2008  cells  (Fig.2H)  also  resulted  in 
reduced  numbers  of  LDs  (Fig. 21).  To  determine  whether  cPLA2  inhibitors  were  able  to 
inhibit  cPLA2  activity,  we  performed  arachidonic  acid  (AA)  release  assay  in  OV202Sh1 
and  Sh2  cells.  Cells  treated  with  lOpM  of  AACOCF3  showed  significant  reduction  of  AA 
release  indicating  that  activity  of  cPLA2  was  inhibited  (Figs.  2J  and  K). 

Inhibition  of  growth  factor  mediated  signaling  reduced  p-cPLA2  and  LD 
biogenesis 

cPLA2  is  activated  by  phosphorylation  at  Ser  residues  505,  515,  and  727. 35  Within 
these  sites,  Ser-505  phosphorylation  leads  to  increased  catalytic  activity,  36  which 
depends  on  MAPK,  ERK  and  p38. 36-38  We  previously  showed  that  downregulation  of 
HSulf-1  in  OV202  cells  results  in  increased  ERK  activity.39  Conversely,  re-expression  of 
HSulf-1  in  HSulf-1  deficient  cells  results  in  attenuation  of  growth  factor  mediated  ERK 
activity.40"42  To  determine  if  inhibition  of  ERK  activity  will  alter  p-cPLA2  levels, 
OV202NTC,  Shi  and  Sh2  cells  were  treated  with  lOpM  U0126  (MEK  inhibitor)  for  24 
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hrs.  As  shown  in  figure  3A  (Panels  2  and  4),  treatment  with  U0126  resulted  in  the 
attenuation  of  p-ERK  and  p-cPLA2ser505  levels  in  OV202Sh1  and  Sh2  cells  compared  to 
NTC  cells.  Bodipy  staining  of  OV202  Shl-and  -Sh2  cells  show  figure  3B  shows 
reduced  LD  accumulation  as  a  result  of  U0126  mediated  inactivation  of  ERK  and  p- 
cPLA2  (Fig.3B). 

To  further  elucidate  the  role  of  growth  factor  signaling,  we  treated  OV202Sh1  cells  with 
PG545,  a  clinically-relevant  HS  mimetic  that  mimics  the  action  of  HSulf-1,43  and 
evaluated  the  expression  levels  of  total  and  p-cPLA2  (Fig  3C).  PG545  functions  as  a 
FIS  mimetic  and  blocks  heparan  sulfate-mediated  growth  factor  signaling  (similar  to 
HSulf-1),  inhibits  angiogenesis  and  carcinogenesis  44  including  OvCa.45  Our  recent 
report  indicated  that  PG545  effectively  inhibited  OV202Sh1  induced  tumorigenesis  and 
also  modulated  levels  of  glycolytic  enzymes  both  in  vitro  and  in  vi vo.43  This  is  the  first 
report  to  demonstrate  that  PG545  (20pM)  treatment  reduces  the  level  of  p-cPLA2, 
indicating  its  effect  on  the  regulation  of  lipid  metabolism.  Consistent  with  this  finding, 
PG545  treatment  reduced  the  number  of  LDs  in  OV202Sh1  cells  (Fig. 3D).  To  evaluate 
the  effect  of  PG545  on  lipid  metabolism  in  vivo  we  determined  the  level  of  p-cPLA2  in 
PG545  treated  OV202Sh1  xenografts  43  by  immunoblot  analysis.  As  shown  in  Fig.3E, 
PG545  reduced  the  phosphorylation  levels  of  p-cPLA2  in  the  treated  xenografts 
compared  to  untreated  controls.  Collectively,  these  results  clearly  indicate  that 
extracellular  growth  factor  signaling  is  directly  associated  with  LD  biogenesis. 

Defective  autophagy  is  associated  with  loss  of  HSulf-1 
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Previous  reports  suggest  that  enhanced  lipid  metabolism  may  lead  to  autophagic 
impairment.32,  46  To  understand  if  HSulf-1  loss  mediated  lipogenic  phenotype  of 
OV202Sh1/Sh2  cells  may  be  associated  with  an  altered  autophagic  response,  we 
analyzed  the  lipidated  form  of  LC3BII  levels  by  immunoblot  analysis.  As  shown  in  figure 
4A,  while  the  OV202NTC  and  the  CI7  cells  (rescue  clone)  with  HSulf-1  expression 
expressed  LC3BII,  HSulf-1  deficient  OV202Sh1  and  -2  cells  had  negligible  levels  of 
LC3BII.  Consistent  with  these  findings,  immunofluorescence  (IFC)  analysis  also 
showed  significantly  less  LC3B  puncta  and  less  acidic  vesicles  as  evidenced  by  acridine 
orange  staining  in  the  Shi  cells  compared  to  NTC  cells  (Fig.4B,  panels  2).  We  next 
analyzed  the  differences  in  autophagic  flux  in  NTC  and  Shi  cells.  Cells  were  grown  in 
the  presence  of  EBSS  medium  (amino  acid  depleted  medium)  to  induce  autophagy  for 
24  hrs  followed  by  1  hr  incubation  with  lOOnM  bafilomycin  to  inhibit  autophagic 
degradation  in  the  lysosomes.  Immunoblot  analysis  showed  no  LC3BII  at  the  basal  level 
in  the  Shi  cells  (Fig.  4C,  lane  1  in  Shi  cells).  However,  induction  of  autophagy  by 
EBSS  treatment  induced  the  expression  of  LC3BII  in  the  Shi  cells  as  evidenced  by 
immunoblot  analysis  (Fig.  4C)  although,  the  level  of  LC3BII  was  much  higher  in  the  NTC 
cells  exposed  to  EBSS  (Fig.4C,  lane  2  compared  to  Lane  2  in  the  Shi  cells).  Inhibition 
of  autophagosome/lysosome  fusion  by  bafilomycin  treatment  resulted  in  higher  levels  of 
LC3BII  both  in  NTC  and  Shi  cells,  although  it  was  more  pronounced  in  the  NTC  cells 
(Lanes  3  in  NTC  and  Shi).  Consistent  with  these  results,  Cyto-ID  staining  that  detects 
autophagosomes  and  autolysosomes  also  indicated  higher  levels  of  autophagy  in  NTC 
cells  compared  to  Shi  cells  (Fig.  4D,  control  compared  to  Fig. 4E, control).  Exposure  of 
both  OV202NTC  and  Shi  cells  to  EBSS  induced  autophagy  as  evidenced  by  increased 
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Cyto-ID  staining  (Figs.  4D  and  E,  panels  2  respectively),  which  was  even  more  evident 
upon  bafilomycin  treatment  (Figs.4D  and  E,  panels  3  respectively),  although  it  was  more 
pronounced  in  the  NTC  cells.  These  data  indicate  that  diminished  autophagic  flux  in 
HSulf-1  depleted  cells  is  associated  with  increased  accumulation  of  lipid  droplets. 

Autophagy  induction  abrogates  LD  biogenesis  in  Shi  cells 

Previous  reports  clearly  showed  that  inhibition  of  AV  in  rat  hepatocytes  resulted  in 
increased  accumulation  of  triglycerides  in  the  form  of  LDs  and  that  activation  of 
autophagy  decreased  LD  levels.29  Based  on  these  and  other  studies,  we  next 
examined  the  consequence  of  inducing  autophagy  in  Shi  cells  on  lipid  droplets.  EBSS 
treatment  that  induces  autophagy  in  Shi  cells  resulted  in  the  reduction  of  LDs 
accumulation  (Fig.4F,  lower  left  panel).  Furthermore,  inhibition  of  autophagy  by 
bafilomycin  rescued  LDs  signifying  autophagy  mediated  lipophagy  may  be  responsible 
for  the  reduction  in  LD  levels  in  cells  deficient  in  HSulf-1  (Fig.4F,  lower  right  panel). 

Previous  reports  by  other  groups33,  47  and  our  data  shown  in  figure  2D  indicated  that 
inhibiting  p-cPLA2  activity  with  AACOCF3  abrogated  LD  biogenesis  in  Shi  cells.  Based 
on  the  notion  that  lipid  accumulation  beyond  a  certain  physiological  level  impaired 
autophagic  function  28,  we  surmised  that  AACOCF3  inhibition  of  p-cPLA2  activity 
resulting  in  the  inhibition  of  LDs,  may  promote  autophagy.  IFC  and  immunoblot  analysis 
of  OV202Sh1  cells  treated  with  AACOCF3  showed  induction  of  autophagy  by  Cyto-ID 
staining  (Fig.4G)  and  increase  in  LC3BII  levels  on  immunoblot  analysis  (Fig. 41,  top 
panel).  Since  it  is  also  well  established  that  increased  growth  factor  signaling  seems  to 
be  associated  with  increased  lipogenesis48'51,  and  our  results  showing  PG545  treatment 
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decreases  the  number  of  LDs  in  these  cells  (Fig. 3D),  we  determined  the  level  of 
autophagy  by  Cyto-ID  staining  following  PG545  treatment  of  OV202Sh1  cells.  There 
was  a  dose  dependent  increase  in  the  Cyto-ID  staining  following  PG545  treatment 
(Fig.4H)  and  increase  in  LC3BII  levels  by  immunoblot  analysis  (Fig.4J).  Notably,  similar 
to  OV202NTC  and  Shi  cells,  HSulf-1  knockout  (KO)  MEFs  showed  negligible 
autophagy  compared  to  wild-type  (WT)  MEFs  (Fig.S2A,  Panel  1).  More  importantly, 
EBSS  treatment  induced  more  autophagy  which  was  even  more  pronounced  upon 
bafilomycin  treatment  in  the  WT  MEFs  compared  to  HSulf-1  KO  MEFs  as  evidenced  by 
Cyto-ID  staining  (Fig.S2A,  panels  2  and  3  respectively).  Also,  AACOCF3  and  PG545 
treatment  of  HSulf-1  KO  MEFs  resulted  in  the  induction  of  autophagy  and  reduction  of 
LDs  in  these  cells  (Fig.S2B). 

Consistent  with  the  notion  that  lipid  accumulation  beyond  a  certain  physiological  level 
impaired  autophagic  function,  addition  of  exogenous  lipids  in  the  form  of  BSA 
conjugated  palmitate  in  OV202NTC  enhances  LD  accumulation  in  the  cytoplasm 
(Fig. S3,  top  panel)  and  inhibits  autophagy  (Fig. S3,  lower  panel).  Additionally,  we  have 
demonstrated  greater  LD  accumulation  in  the  ATG5  knockdown  MEFs  cells  compared 
to  the  wild  type  MEFs,  whereas  ATG5  knockdown  cells  showed  a  lower  level  of  Cyto-ID 
staining  compared  to  the  wild  type  cells  (Fig  S4). These  results  clearly  show  a 
relationship  between  induction  of  autophagy  and  abrogation  of  LDs. 

Although  increase  in  LC3BII  is  indicative  of  autophagic  flux,  a  more  reliable  assay  is  the 
use  of  fluorescent  tagged  Cherry-GFP-LC3B  construct  to  monitor  autophagic  flux  52, 53 
following  autophagy  induction.  As  the  GFP  signal  is  pH  sensitive  and  is  quenched  in  the 
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lysosomes,  autophagosomes  light  up  as  both  red  and  green  while  autolysosmes  are 
distinguished  mostly  by  red  signals  (seen  as  orange  to  yellow  in  merged  images)  only. 
Following  transient  transfection  (48  hrs)  of  Cherry-GFP-LC3B  in  OV202Sh1  cells,  the 
cells  were  treated  with  lOpM  of  AACOCF3  and/or  PG545.  Induction  of  autophagy  by 
these  agents  caused  an  increase  in  Cherry+/GFP+  signals  (Fig.4K,  panels  2  and  3)  as 
well  as  Cherry+  signals  (Panels  4  in  the  treated  cells)  compared  to  un-treated  control 
cells  (Panel  1 )  clearly  indicating  the  activation  of  autophagic  flux. 

AACOCF3  alone  and  in  combination  with  carboplatin  suppresses  tumor  growth, 
and  inhibits  lipid  droplet  biogenesis  in  vivo 

To  determine  whether  blocking  cPLA2  activities  via  AACOCF3  confer  any  advantage  in 
attenuating  tumor  growth  in  combination  with  chemotherapeutic  agent  such  as 
carboplatin,  we  performed  testing  of  these  agents  in  vivo  using  a  mouse  model.  The 
effect  of  AACOCF3  alone  and  in  combination  with  carboplatin  (CBP)  on  primary  tumor 
growth  was  evaluated  in  OV202Sh1  cells  bearing  nude  mice  i.p.  (intraperitoneally). 
5X1 06  Sh  clones  expressing  luciferase  cells  in  serum-free  RPMI  1640  were  injected 
intraperitoneally  into  female  athymic  nu/nu  mice  at  4  to  5  wk  of  age  (National  Cancer 
Institute-Frederick  Animal  Production  Area  (Frederick,  MD).  Once  i.p  implants  were 
visible  via  non-invasive  imaging  (4  days  after  inoculation),  mice  were  randomized  (10 
mice/group)  and  treated  i.p.  with  the  lOmg/kg  body  weight  (b.w.)  of  cPLA2  inhibitor 
AACOCF3  (referred  to  as  F3  in  the  figures)  every  other  day  till  the  end  of  the  study, 
51  mg/kg  b.w.  carboplatin  (CBP)  every  5  days  till  the  end  of  the  study  and  a  combination 
of  CBP+F3  every  5  days,  with  vehicle  treated  group  as  control.  Luciferase  imaging  of 
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the  representative  mice  from  four  groups  is  shown  in  Fig.5A.  Higher  luciferase  intensity 
in  the  control  group  and  CBP  group  indicates  increased  tumor  volume,  progression  and 
metastasis.  Representative  excised  tumor  from  one  of  the  mouse  in  each  group  is 
shown  in  Fig.5B.  At  necropsy,  the  measured  tumor  weight  and  circumference  of  mice 
in  each  group  showed  that  the  combination  treatment  was  more  effective  at  reducing 
cancer  progression  compared  to  all  of  the  other  groups  (Figs.5C-D).  There  was  no 
significant  body  weight  loss  in  F3  or  combination  treated  group  compared  to  control 
group  suggesting  both  F3  and  CBP  treatment  were  well  tolerated  without  apparent 
toxicity  to  the  animal  (Fig.5E).  More  importantly,  Bodipy  staining  of  frozen  sections  of 
the  xenografts  clearly  showed  high  levels  of  lipid  droplets  (Fig.5G)  in  the  control  and 
the  CBP  treated  group  and  significantly  less  LDs  in  the  F3  and  combination  treated 
group  consistent  with  the  in  vitro  data  shown  in  figure  2D,  top  panel.  In  contrast, 
TUNEL  staining  showed  increased  TUNEL  staining  in  the  F3  and  the  combination 
treated  groups  compared  to  untreated  control  or  CBP  treated  groups  alone  (Fig.5H) 
Anti-tumor  effects  of  AACOCF3  monotherapy  and  combination  with  CBP  correlated  with 
significant  reductions  in  tumor  cell  proliferation  marker,  Ki67  and  phospho-cPLA2  as 
determined  by  IHC  staining  using  the  specific  antibodies  to  these  proteins  (Figs.  51  and 
K). 

Discussion: 

Lipid  droplet  biogenesis  is  a  mechanism  by  which  cells  store  excess  free  fatty  acids 
(FFA)  in  the  form  of  triacylglycerides  and  cholesterylesters.54  Under  stress  conditions, 
FFAs  are  released  by  lipolysis  from  the  LDs  and  undergo  beta  oxidation  in  the 
mitochondria  to  provide  energy  to  the  dying  cell.3  During  oncogenic  transformation, 
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many  cancer  cells  manifest  a  high  rate  of  de  novo  lipid  synthesis  which  results  in  the 
accumulation  of  cytoplasmic  lipid  droplets  (LDs).  40,  48, 55  We  recently  reported  that  loss 
of  HSulf-1  promotes  a  “lipogenic  phenotype”  by  upregulating  lipid  related  metabolites, 
fatty  acid  synthesis  and  beta-oxidation  indicating  an  important  role  of  HSulf-1  in  the 
metabolic  regulation.2  In  this  study  we  have  identified  HSulf-1  as  a  molecular  link  -  in 
regulating  LDs  accumulation  by  modulating  autophagy.  Using  HSulf-1  ShRNA 
downregulated  OV202  and  TOV2223  ovarian  cancer  cells,  we  observed  that  HSulf-1 
depleted  cells  accumulated  LDs  and  exhibited  diminished  autophagic  flux  as  evidenced 
by  reduced  autophagic  vesicles  (AVs).  Autophagy  maintains  the  cellular  homeostasis  by 
recycling  cellular  components,  clearance  of  lipids  and  proteins  26.  Previous  studies  have 
shown  that  increased  growth  factor  signaling  may  lead  to  defective  autophagy  56  which 
is  predominant  in  cells  with  loss  of  HSulf-1.  Our  data  indicate  that  HSulf-1  expressing 
cells  attenuate  growth  factor  signaling  and  ERK-1/2  activation  thereby  limiting  the  direct 
activation  of  cPLA2,  a  key  protein  involved  in  lipid  biogenesis.  Limited  growth 
promoting  signaling  and  nutrient  depletion  are  two  positive  regulators  of  autophagy.  Our 
results  show  that  by  genetic  depletion  of  HSulf-1,  autophagic  flux  was  significantly 
decreased  and  lipid  droplets  were  significantly  increased.  Additional  evidence  indicates 
that  autophagic  clearance  of  lipid  droplets  maintains  the  cellular  lipid  homeostasis  in  a 
steady  state;  and  impaired  autophagy  leads  to  accumulation  of  lipids  in  the  form  of  LDs 
28,  48, 57  w|-|jc|-|  subsequently  leads  to  growth  and  progression  of  the  tumor.  Consistent 

with  these  reports  our  results  suggest  that  HSulf-1  loss  leads  to  increased  LD 
biogenesis,  defective  autophagic  clearance  of  cytoplasmic  LDs,  thus  leading  to  a 
lipogenic  phenotype  as  we  previously  reported  2  providing  favorable  conditions  for  cell 
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growth  and  tumorigenesis.  We  previously  reported  that  loss  of  HSulf-1  also  promotes 
enhanced  fatty  acid  oxidation  in  Shi  and  Sh2  cells  compared  to  NTC  cells.  Our  data 
also  suggest  that  the  free  fatty  acids  (FFA)  released  by  the  action  of  lipases  may  be 
immediately  converted  to  neutral  lipids  and  stored  in  lipid  droplets  to  be  released  under 
stress  conditions.  Additionally,  we  have  also  demonstrated  that  supplementation  of  BSA 
conjugated  palmitic  acid  in  the  OV202  NTC  cells  promotes  cytoplasmic  LD 
accumulation  but  inhibits  autophagy,  which  indicates  an  inverse  association  of  LD 
biogenesis  and  autophagy. 

This  is  the  first  report  linking  loss  of  HSulf-1  with  increased  cPLA2  activity  that  triggered 
LD  biogenesis.  Consistent  with  LD  accumulation,  OV202  Sh  cells  showed  an  increased 
AA  release  compared  to  NTC  cells.  Released  AA  is  critical  for  inducing  inflammatory 
responses  in  cancer  cells  though  cyclooxygenase  and  lipoxygenase  pathways.58 
Pharmacological  inhibition  of  cPLA2  phosphorylation  reduces  AA  release  and  LD 
biogenesis,  suggesting  that  loss  of  HSulf-1  mediated  LD  biogenesis  and  AA  release  is 
dependent  on  cPLA2  activity. 

Additionally,  treatment  with  PG545,  a  HS  mimetic  that  mimics  the  action  of  HSulf-1  in 
inhibiting  heparan  sulfate  binding  growth  factor  signaling  43  reduced  the  number  of  LDs, 
enhanced  autophagy,  reduced  tumor  growth  43  and  reduced  p-cPLA2  levels. 
Significantly,  EBSS  treatment  promoted  autophagy  and  reduction  in  LDs  in  HSulf-1  KO 
MEFs  suggesting  the  role  of  HSulf-1  in  regulating  these  two  interconnected  metabolic 
pathways  in  vivo.  Although  a  recent  report  by  Shteingauz  et  al.,  59  using  heparanase- 
overexpressing  (Hepa)  U87  glioma  cells  showed  that  PG545  inhibited  heparanase 


15 


induced  autophagy,  our  results  clearly  shows  PG545  as  a  promoter  of  autophagy  in 
OvCa  cells  further  affirming  the  role  of  role  of  autophagy  as  a  promoter  of  cell  death 
(and/or  a  survival  mechanism)  is  context  and  cell  type  dependent  phenomenon. 

Our  in  vivo  data  using  a  specific  cPLA2  inhibitor,  AACOCF3  showed  less  tumor  burden 
compared  to  the  untreated  control  indicating  that  a  higher  level  of  cPLA2  itself  can  play 
a  crucial  role  in  ovarian  tumorigenesis.  Interestingly,  AACOCF3  and  carboplatin 
combination  therapy  was  more  effective  to  reduce  tumor  growth  than  carboplatin  alone. 
This  is  a  first  report  to  show  that  HSulf-1  loss  associated  ovarian  tumor  formation  is 
regulated  by  higher  activity  of  cPLA2.  More  studies  are  required  to  test  the  synergistic 
effect,  absorption  and  bioavailability  of  AACOCF3  alone  versus  the  combination  with 
carboplatin  in  the  mouse  model.  Our  findings  indicate  that  cPLA2  inhibition  is  a  potent 
therapeutic  approach  to  synergize  with  conventional  chemotherapies  to  reduce  ovarian 
tumor  growth. 


Materials  and  Methods: 

Cell  culture 

Human  ovarian  cancer  cell  line  OV202  NTC,  OV202  Shi,  OV202  Sh2  and  OV202  CI7 
were  grown  in  minimum  essential  medium  alpha  IX  (Cellgro,  Cat#  15-012-CV) 
supplemented  with  20%  fetal  bovine  serum  (Biowest,  Cat  #  SI  620)  and  1%  pen-strep. 
OV202  clones  and  NTC  are  grown  in  the  presence  of  1pg/ml  puromycin  as  a  selection 
marker  for  the  HSulf-1  shRNA.  Growth  media  for  CI7  cells  are  additionally 
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supplemented  with  400pg/ml  of  G418  as  a  selection  marker  for  overexpression  plasmid 
of  HSulf-1. 

MEF  (WT  and  HSulf-1  knock  out)  cells  are  obtained  from  Germany.  They  are  grown  in 
the  DMEM  IX  medium  in  the  presence  of  10%  FBS  and  1  %  Pen  strep  with  addition  of 
650  pg/ml  of  G418.  OV2008  cells  were  grown  in  RPMI  1640  medium  (Cellgro,  Cat#  10- 
040-CV)  supplemented  with  10%  fetal  bovine  serum  and  1%  pen-strep. 

Inhibitors 

AACOCF3,  an  inhibitor  of  phospho-cPLA2  (Cayman  chemicals,  Cat  #  62120)  was  used 
to  inhibit  the  cPLA2  activity  and  lipid  droplet  biogenesis.  MEK  inhibitor  U0126  (Sigma- 
Aldrich,  Cat  #11120)  was  used  to  inhibit  ERK  activity.  PG545,  a  heparanase  inhibitor 
used  to  block  the  effect  of  extracellular  growth  factor  signaling  (Progen 
pharmaceuticals). 

Cell  proliferation 

Cell  proliferation  was  quantified  by  MTT  assay  (Promega,  Cat  #G4000)  as  previously 
described  60. 

Generation  of  cPLA2  downregulated  Stable  clone 

OV202  Shi  cell  line  and  OV2008  cells  were  transfected  with  Sh  cPLA2  (ShCPLA2  1 : 

T  GAT  ACAAAT  GT  AG  G  GAT  AT  A  and  ShCPLA2  2:  CCTTGTATTCTCACCCTGATT) 
using  standard  transfection  protocol  and  reagents.  Stable  closes  were  further  selected 
by  puromycin.  OV202  cells  were  cultured  in  5%  C02-95%  air  humidified  atmosphere  at 
37°C  with  minimal  essential  medium  supplemented  with  20%  fetal  bovine  serum  with 
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their  specific  antibiotic  selection.  All  cell  lines  were  tested  using  a  PCR-based  assay 
and  found  to  be  free  of  Mycoplasma  contamination. 


Western  blot  analysis 

Western  blot  analysis  was  performed  as  described  previously8.  Whole  cell  lysates  were 
analyzed  with  the  following  antibodies:  HSulf-1  (Abeam,  Cat  #  ab31960),  phospho- 
cPLA2  (Sigma,  Cat  #SAB4503812),  total-cPLA2  (SantaCruz,  Cat  #SC454),  phospho- 
ERK  (Cell  Signaling,  Cat  #9106),  total-ERK  (Cell  Signaling,  Cat  #9102),  LC3B  (Cell 
Signaling,  Cat  #3868),  GAPDH  (Cell  Signaling,  Cat  #3683),  (3-actin  (  Gene  Tex,  Cat  # 
GTX629630)  and  (3-tubulin  (GeneTex,  Cat  #GTX1 09639  ). 

Detection  of  autophagic  vesicles  by  Cyto-ID  staining 

OV202  NTC  and  Shi  cells  were  seeded  in  the  4-well  tissue  culture  chambered  slides 
and  treated  with  indicated  concentrations  of  AACOCF3  or  starved  with  EBSS  in  the 
presence  or  absence  of  bafilomycin.  Then,  autophagic  vesicles  were  detected  with  a 
Cyto-ID  autophagy  detection  kit  (Enzo  Life  Sciences,  Cat  #ENZ-51031).  Images  were 
captured  using  a  Zeiss  LSM  510  microscope. 

Detection  of  acidic  vesicular  organelles  using  Acridine  orange  staining 

OV202  NTC  and  Shi  cells  were  seeded  in  4-well  tissue  culture  chambered  slides  and 
grown  overnight.  Next  day  cells  were  stained  with  acridine  orange  (Ipg/ml;  Molecular 
Probes,  Cat  #A1301)  in  a  serum  free  medium  for  15  minutes  at  37°  C.  Acridine  orange 
was  removed  and  images  were  captured  in  Zeiss  LSM  microscope.  Red  fluorescence 
indicates  the  acidic  vesicles  as  described  earlier61. 
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Transmission  Electron  microscopy 

Cultured  cells  were  washed  twice  with  PBS  and  then  fixed  in  Trumps  fixative  containing 
4%  formaldehyde  and  1%  glutaraldehyde  in  a  phosphate  buffer  pH  —7.3,  post-fixed  in 
1.0%  0s04,  dehydrated  with  ethanol  gradation,  and  transitioned  into  propylene  oxide 
for  infiltration  and  embedding  into  spurr  epoxy  resin.  Ultrathin  sections  were  cut  onto 
grids,  stained  with  uranyl  acetate  and  lead  citrate,  and  examined  with  aJEM-1400 
(JEOL  USA)  transmission  electron  microscope,  and  digitally  photographed. 
Morphologically  we  assessed  the  cytoplasmic  lipid  droplets  and  autophagosomes  and 
the  numbers  were  counted.  A  total  of  25  cells  per  sample  were  analyzed. 

Bodipy  staining 

Cells  (50,000)  were  seeded  on  a  coverslip  in  a  24-well  plate  and  were  grown  for  24  hrs 
in  the  presence  of  complete  growth  medium.  Cells  were  washed  and  fixed  in  4% 
paraformaldehyde  for  ten  minutes  at  room  temperature  before  staining  with  Ipg/ml 
Bodipy  (493/503;  Sigma,  Cat  #790389)  in  PBS  for  ten  minutes  at  room  temperature. 
Coverslips  were  washed  with  PBS  and  mounted  in  a  slide  with  Prolong  Gold  Antifade 
Reagent  (Invitrogen,  Cat  #P36935).  Bodipy  stained  cells  were  examined  under  Zeiss- 
LSM  510  fluorescence  microscope. 

Immunoflouroscence 

Cells  were  grown  on  coverslips  overnight,  and  then  fixed  with  100%  methanol  followed 
by  blocking  with  1%  BSA  in  PBS.  After  blocking,  the  cells  were  incubated  with 
respective  primary  antibodies  at  room  temperature  for  1  hr,  then  washed  three  times 
with  IX  PBS  and  then  incubated  in  dark  with  alexa  fluor  goat  anti-rabbit  (488  nm) 
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or/with  alexa  fluor  rabbit  anti-mouse  (593  nm)  in  1%  BSA  in  PBS.  Coverslips  were 
washed  three  times  before  mounting  with  Prolong  Gold  Antifade  reagent  (Invitrogen). 
Stained  samples  were  visualized  using  a  Zeiss-LSM  510  fluorescence  microscope. 


Arachidonic  acid  release  assay 

Arachidonic  acid  (AA)  release  was  evaluated  in  OV202  NTC,  Shi,  Sh2,  CI7  cells  and 
OV2008  vector  and  clones.  In  a  6-well  plate  0.5x1 06  cells  were  plated  and  grown 
overnight  in  complete  media.  On  the  next  day,  growth  medium  was  removed  and  the 
cells  were  incubated  with  3H-AA  under  serum  free  medium  for  24  hrs.  The  media  was 
removed  and  cells  were  washed  three  times  with  IX  PBS  and  complete  growth  medium 
was  added  with  or  without  AACOCF3  (lOpM)  for  another  24  hrs.  Aliquots  of  growth 
medium  from  the  cells  were  measured  for  radioactivity  at  the  end  of  treatment. 
Radioactivity  is  showed  as  counts  per  minute  (cpm/ml). 

Nude  mouse  tumor  model 

5x1 06  OV202  Shi  cells  expressing  luciferase  were  injected  intraperitoneally  (i.p)  to  the 
nude  mice.  4  days  following  i.p  inoculation  the  mice  were  treated  with  (lOmg/kg) 
AACOCF3  every  3rd  day  until  the  end  of  study,  or  carboplatin  (51  mg/kg)  in  every  5th 
day  for  4  cycles,  or  combination  of  carboplatin  (51  mg/kg)  and  AACOCF3  (lOmg/kg) 
every  5th  day  for  4  cycles.  All  the  animals  were  sacrificed  by  the  28th  day  of  the 
experiment.  Mayo  clinic  animal  care  and  use  committee  approved  all  procedures. 

In  vivo  bioluminescence  measurements 
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Weekly  bioluminescent  reporter  imaging  was  performed  to  monitor  the  seeding  of 
OV202Sh1-luc  cells,  using  the  Xenogen  I  VIS  200  System  (Caliper  life  sciences)  as 
previously  described45. 

Immunohistochemistry:  Immunohistochemical  (IHC)  staining  was  done  on  tumors 
from  control,  F3,  CBP  and  F3+CBP  treated  mice.  Paraffin  embedded  tissue  sections 
from  the  tumors  were  stained  for  Ki67  and  phospho-cPLA2.  The  standard  IHC 
procedure  was  performed  at  the  Pathology  Research  Core  (Mayo  Clinic,  Rochester, 
MN)  as  previously  described45. 

TUNEL  Staining:  The  terminal  dUPT  nick  end-labeling  (TUNEL)  assay  of  frozen  tissue 
sections  of  tumors  from  control,  F3,  CBP  and  F3+CBP  treated  mice  were  performed 
using  ApopTag  Fluorescein  Direct  In  Situ  Apoptosis  Detection  Kit  (Millipore,  Cat  # 
S7110)  according  to  the  manufacturer’s  instructions.  Green  fluorescence  indicates 
TUNEL  staining  and  red  fluorescence  indicates  propidium  iodide  (PI)  staining. 

Statistical  analyses 

All  results  were  expressed  as  the  mean  ±  S.D.  data  obtained  from  three  separate 
experiments.  All  statistical  analysis  was  evaluated  using  Graph  pad  Prism  software 
(San  Diego,  CA).  Data  were  analyzed  by  the  paired  t  test,  and  P  values  less  than  0.05 
or  mentioned  otherwise,  was  considered  statistically  significant. 
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Figure  Legends 

Figurel:  HSulf-1  loss  promotes  lipid  droplet  biogenesis  and  defective  autophagy 
in  OV202,  TOV2223  and  in  HSulf-1  knockout  mouse  embryonic  fibroblast  cells 

A.  Protein  expression  of  HSulf-1  was  assessed  by  western  blot  in  OV202  clones 
(NTC,  Shi,  Sh2  and  CI7  cells).  Actin  was  probed  as  a  loading  control.  B.  OV202-NTC, 
Shi,  Sh2  and  CI7  cells  were  labeled  with  Bodipy  493/503  (green)  and  DAPI  (Blue)  to 
determine  the  cytoplasmic  lipid  droplets  (LD)  and  nuclei  respectively.  C.  Mean 
fluorescence  of  Bodipy  in  NTC,  Shi,  Sh2  and  CI7  is  measured  using  Carl  Zeiss  Zen 
(2009)  software  and  the  data  is  plotted  as  a  bar  diagram.  D.  Representative  trans¬ 
electron  microscopic  (TEM)  images  of  NTC,  Shi,  Sh2  and  CI7  cells  are  shown; 
cytoplasmic  LDs  are  marked  with  green  arrows,  whereas  red  arrows  indicate 
autophagic  vesicles  (AV).  E.  LDs  and  AVs  quantified  in  NTC,  Shi  and  Sh2  cells  (total 
25  cells  counted  in  each  category)  demonstrated  in  bar  graph.  F.  Protein  expression  of 
HSulf-1  in  TOV2223  NTC  and  Shi  cells  is  detected  by  immunoblot  analysis;  (3-actin  is 
probed  as  a  loading  control.  G.  Bodipy  staining  performed  in  TOV2223  NTC  and  Shi 
cells  to  detect  LDs.  H.  Representative  TEM  images  of  TOV2223  NTC  and  Shi  cells 
were  shown.  I.  AVs  and  LDs  in  TOV2223  cells  are  quantified  J.  Protein  expression  of 
HSulf-1  in  Mouse  Embryonic  Fibroblast  (MEF)  wild  type  (WT)  and  HSulf-1  knock  out 
(Sulf-1  _/')  cells  is  detected  by  western  blotting.  (3-Tubulin  has  been  used  as  a  loading 
control.  K.  MEF  cells  were  labeled  with  Bodipy  493/503  (green)  to  determine  the 
cytoplasmic  lipid  droplets  (LD). 

Figure  2:  Inhibition  of  cPLA2  attenuates  LD  biogenesis  in  OV202SM  cells. 
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A.  Immunoblot  analysis  shows  phospho-cPLA2  and  total-cPLA2  levels  in  OV202  NTC, 
Shi,  Sh2  and  CI7  cells.  B.  &  C.  Immunoblot  analysis  shows  phospho-cPLA2  and  total- 
cPLA2  levels  in  OV202  Shi  and  OV2008  cells  after  treatment  with  lOpM  and  20pM  of 
cPLA2  inhibitor,  AACOCF3.  D.  Shi  and  SH2  cells  were  treated  with  either  10[iM 
AACOCF3  or  MAFP  for  24  hrs  followed  by  Bodipy  (493/503)  staining  to  show  the  LDs. 
E.  Representative  TEM  images  of  OV202  Shi  cells  untreated  and  treated  with  10pM 
AACOCF3  and  MAFP  for  24  hrs.  LDs  are  marked  with  green  arrows,  whereas  red 
arrows  indicate  autophagic  vesicles  (AV).  LDs  are  quantified  from  25  cells  from 
untreated  and  treated  cells  from  TEM  images  and  demonstrated  in  a  bar  diagram.  F.  & 
H.  Total  cPLA2  expression  was  transiently  downregulated  with  two  different  ShRNAs 
targeting  the  open  reading  frame  of  cPLA2  in  OV202Sh1  and  OV2008  cells. 
Immunoblot  analysis  shows  total  cPLA2  expression  in  non-targeted  control  transduced 
(NTC),  ShcP-1  and  ShcP-2  cells.  G.  &  I.  Bodipy  staining  of  LDs  in  OV202  Shi  and 
OV2008  cells  after  transiently  transfecting  them  with  either  empty  vector  or  sh-cPLA2. 
J.  Arachidonic  acid  (AA)  release  is  evaluated  in  OV202  NTC,  Shi,  Sh2  and  CI7  cells 
treated  or  untreated  with  lOpM  AACOCF3.  Cells  were  incubated  with  H3-AA  under 
serum  starved  condition  for  24  hrs.  Fresh  medium  was  added  to  the  cells  after  washing 
and  aliquots  of  growth  medium  were  measured  for  radioactivity  after  24  hrs. 
Radioactivity  is  shown  as  counts  per  minute  (cpm)/ml.  K.  AA  release  is  evaluated  in 
OV202  Shi  and  CI7  cells. 

Figure  3:  Inhibition  of  growth  factor  mediated  signaling  reduced  p-cPLA2  and  LD 
biogenesis 
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A.  OV202  NTC,  Shi  and  Sh2  cells  were  treated  with  20pm  of  U0126,  an  MEK  inhibitor 
for  24  hrs  followed  by  western  blot  analysis  to  detect  the  protein  expressions  of 
phospho-ERK,  total-ERK,  phospho-cPLA2  and  total-cPLA2.  B.  Bodipy  staining  for  LDs 
in  OV202Sh1  and  Sh2  cells  after  treatment  with  20  pM  U0126.  C.  Immunoblot  analysis 
of  phospho-cPLA2  and  total-cPLA2  in  OV202Sh1  cell  is  shown  after  treatment  with  10 
and  20  pM  PG545.  D.  Bodipy  staining  of  LDs  in  OV202  Shi  cells  after  treating  with  10 
and  20pm  of  PG545.  E.  Immunoblot  analysis  of  phospho-cPLA2  and  total-cPLA2  in 
lysates  from  untreated  and  PG545  treated  OV202Sh1  Xenografts. 

Figure  4:  Dysfunctional  autophagy  is  associated  with  loss  of  HSulf-1  and 
activation  of  cPLA2. 

A.  Immunoblot  analysis  with  anti-LC3B  antibody  in  OV202NTC,  Shi,  Sh2  and  Sh1CI7 
cells.  B.  Immunoflourescence  analysis  of  OV202  NTC  and  Shi  cells  with  anti-LC3B 
antibody  (green,  Top  panels),  Acridine  orange  (Red,  Lower  panels)  and  DAPI  (blue). 

C.  Immunoblot  analysis  of  LC3BI  and  II  in  OV202NTC  and  Shi  cells  treated  with  EBSS 
with  or  without  bafilomycin  A1  (lOOnm)  for  1  hour  with  GAPDH  as  loading  control. 
Cyto-ID  staining  of  D.  OV202NTC  E.  OV202  Shi  cells  in  control,  EBSS  and  EBSS+Baf 
treated  cells.  Cells  were  grown  in  either  regular  media  or  EBSS  to  induce  autophagy.  F. 
Bodipy  staining  of  LDs  in  OV202Sh1  cells  treated  with  EBSS  with  or  without  bafilomycin 
A1  (lOOnm)  for  1  hour.  Cells  grown  in  complete  media  is  used  as  controls  (-EBSS). 
Immunoflourescence  analysis  of  Cyto-ID  staining  of  OV202Sh1  cells  treated  with  10 
and  20pM  G.  AACOCF3  and  H.  PG545  showing  increase  in  autophagy.  Immunoblot 
analysis  of  increase  in  LC3BII  levels  in  OV202Sh1  cells  treated  with  increasing 
concentrations  of  I.  AACOCF3  and  J.  PG545.  K.  Autophagic  flux  in  AACOCF3  and 


28 


PG545  treated  OV202Sh1  cells  following  transient  expression  of  Cherry-GFP-LC3B  (48 
hrs).  The  cells  were  treated  for  24  hrs.  Confocal  microscopy  analysis  shows  changes 
in  the  autolysosome  formation.  Green  (GFP),  and  red  fluorescence  (Cherry)  signals 
and  the  merged  images  show  orange  signals  in  the  autolysosomes  in  the  treated  cells. 
Representative  pictures  from  two  independent  experiments  are  shown. 

Figure  5:  AACOCF3  alone  and  in  combination  with  carboplatin  suppresses  tumor 
growth,  and  inhibits  lipid  droplet  biogenesis  in  vivo 

A.  Randomized  OV202Sh1  tumor-bearing  mice  (n=10),  were  treated  with  water  or 
AACOCF3  (F3)  (lOmg/kg),  or  carboplatin  (CBP)  (51  mg/kg)  or  with  a  combination  of 
AACOCF3  and  CBP  for  4  weeks  and  after  28  days  mice  were  euthanized. 
Representative  images  of  the  mice  prior  to  sacrifice  from  control,  F3  treated,  CBP 
treated  and  CBP+F3  combination  treatment  groups  using  I  VIS  luminescence  imaging 
system  series  2000.  Color  bar  shows  photon  intensity.  B.  Representative  images  of 
excised  tumors  from  OV202  Shi  xenografts.  C.  Excised  tumor  weights  from  control,  F3, 
CBP  and  CBP+F3  combination  treatment  groups.  D.  Abdominal  circumference  from 
different  groups  of  mice  measured  on  the  day  of  sacrifice.  E.  Total  body  weight  of 
untreated  control  and  treated  with  F3,  CBP,  CBP+F3.  F.  Immunoblot  analysis  of 
LC3BI/II  in  lysates  from  untreated,  F3,  CBP  and  F3+CBP  treated  Xenografts  with 
GAPDH  as  loading  control.  G.  Tumor  xenografts  were  stained  with  Bodipy  to  detect  lipid 
droplets;  DAPI  was  used  to  stain  the  nuclei.  H.  Representative  images  of  TUNEL 
staining  of  frozen  section  of  xenografts  from  different  treatment  groups  and  untreated 
control.  Green  fluorescence  indicates  TUNEL  and  red  fluorescence  indicates  propidium 
iodide  (PI)  staining.  I.  Representative  images  of  immunohistochemical  staining  of  Ki67, 
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J.  Bar  graph  showing  the  Ki67  proliferation  index  (*p=<0.05  and  **p=<0.01). 
Proliferation  index  (PI)  of  Ki-67  staining  was  measured  using  ImmunoRatio,  from  public 
domain  image  analysis  software.  K.  phospho-cPLA2  in  FFPE  sections  of  treated  and 
untreated  xenografts. 

Supplementary  figure  legends 
Fig  SI 

A.  Protein  expression  of  HSulf-1  was  determined  by  western  blot  in  SKV03  vector  and 
cl7  cells.  GAPDH  was  probed  as  loading  control  B.  Representative  electron  micrograph 
images  of  SKOV3  vector  and  cl7  cells  are  shown;  LDs  and  AVs  are  indicated  by  green 
and  red  arrows  respectively.  C.  Quantification  of  LDs  and  AVs  from  electron  micrograph 
images  are  determined  in  SKOV3  vector  and  cl7  cells  and  are  demonstrated  in  a  bar 
diagram.  D.  SKOV3  vector  and  cl7  cells  were  labeled  with  Bodipy  493/503  (green)  and 
DAPI  (blue)  to  determine  the  cytoplasmic  LDs  and  nuclei  respectively. 


Fig  S2 

A.  Mouse  Embryonic  Fibroblast  (MEF)  wild  type  (WT)  and  HSulf-1  knock  out  (Sulf-1  _/') 
cells  were  either  grown  in  complete  media  (CM)  or  starved  with  EBSS  in  the  presence 
and  absence  of  bafilomycin  A1  followed  by  Cyto-ID  staining  to  demonstrate  the 
autophagic  vesicles.  B.  MEFs  were  treated  with  10pm  of  AACOCF3  and  PG545 
respectively,  and  then  stained  with  either  Bodipy  or  Cyto-ID  to  detect  LDs  and  AVs. 

Fig  S3 
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OV202NTC  cells  were  treated  with  25  and  50|iM  of  BSA  conjugated  sodium  palmitate 


for  24  hrs.  After  treatment,  cells  were  labeled  with  Bodipy  and  Cyto-ID  to  detect  lipid 
droplets  and  autophagic  vesicles  respectively. 

Fig  S4 

Wild  type  (WT)  and  ATG  5  knockdown  mouse  embryonic  fibroblast  (MEF)  cells  were 
stained  with  Bodipy  and/or  Cyto-ID  to  identify  lipid  droplets  and  autophagic  vesicles. 
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